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SECTION  I 
INTRODUCTION 

This  report  describes  an  analysis  of  the  sliding 
phenomena  that  occur  during  operation  of  a hypervelocity  guide- 
rail  system.  As  such  systems  come  into  wide  use  in  the  hyper- 
velocity experimental  field  it  becomes  increasingly  important 
to  understand ‘ the  behavior  of  the  sliding  contact  between  the 
projectile  and  rails.  Very  large  guide-rail  systems  (both  in 
bore  diameter  and  length)  are  currently  being  contemplated. 

The  technical  feasibility  of  such  systems  will  depend,  among 
other  things,  on  the  magnitude  of  the  frictional  forces  and 
associated  projectile  wear,  and  other  behavior  peculiar  to  hyper- 
velocity sliding  contacts  such  as  gouging.  This  study  was  initiated 
to  investigate  friction,  wear,  and  gouging  in  hypervelocity 
guide-rail  systems.  The  objective  of  the  investigation  was  to 
develop  an  analysis  of  hypervelocity  sliding  which  could  be 
used  to  evaluate  frictional  forces,  projectile  wear  and  the 
probability  of  gouging  behavior  for  both  existing  and  proposed 
guide-rail  systems.  To  maintain  the  required  generality  a funda- 
mental approach  was  employed.  The  basic  phenomena  were  identified, 
analyzed  and  combined  into  an  overall  interface  model.  This 
model  was  then  combined  with  the  projectile/rail  dynamics  to 
provide  the  overall  friction,  wear,  and  gouging  evaluation  tool. 

This  section  begins  with  a general  description  of  guide - 
rail  systems  and  a detailed  description  of  the  Arnold  Engineering 
Development  Center  (AEDC)  Hypervelocity  Tracks  K and  G,  which 
were  the  specific  systems  used  to  exercise  the  model.  It  must 
be  emphasized  that  our  analysis  is  not  restricted  in  any  sense 
to  the  track  K or  track  G parameters  and  can  be  applied  to  any 
guide- rail  configuration  now  in  use,  under  construction,  or  being 
contemplated.  An  introduction  to  the  projectile/rail  sliding 
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phenomena  which  will  be  treated  is  then  presented.  Finally, 
the  effects  of  projectile  dynamics  are  outlined. 

The  remainder  of  the  report  deals  in  detail  with  the 
sliding  phenomena,  the  interface  model,  the  projectile  dynamics, 
a computer  program  designed  to  describe  the  performance  of 
hypervelocity  guide^rail  systems,  and  comparisons  between  results 
computed  with  this  program  and  data  from  experimental  firings. 

1.1  GUIDE  RAIL  SYSTEMS 

A guide-rail  system  is  a device  which  basically  takes  a 
projectile  as  delivered  from  the  muzzle  of  a gun  and  guides  it 
via  a fixed  rail  system  through  a prescribed  atmospheric  en- 
vironment. The  application  under  consideration  is  a hypervelocity 
guide-rail  system  employed  in  a reentry  vehicle  simulation  and 
test  facility.  A model  (the  projectile)  is  launched  through  a 
typical  reentry  environment  and  undergoes  ablation  and  erosion 
similar  to  that  encountered  during  an  actual  reentry. 

The  principal  advantage  of  a guide-rail  system  is  that 
it  provides  precise  control  of  projectile  orientation  and  location 
during  flight.  This  control  permits  satisfactory  ballistic  ex- 
perimentation with  projectiles  that  have  at  least  some  level  of 
aerodyncunic  instability  and  lift.  Free  flight  ranges  are  limited 
to  testing  with  aerodynamically  stable  and  zero  lift  projectiles 
in  order  to  keep  the  projectile  on  a reasonably  straight  and  pre- 
dictable trajectory  such  that  its  behavior  can  be  closely  ob- 
served. This  limitation  becomes  critically  important  for  tasks 
such  as  erosion  testing  where  non-symmetrical  nosetip  mass  loss 
can  generate  lift  and  subsequent  projectile  fly-off. 

Guide  rail  systems  do,  however,  perturb  the  aerodynamic 
flow  around  projectiles  in  flight.  These  perturbations  reduce 
the  area  of  the  forebody  that  can  be  used  for  testing  and  cause 
increased  drag  on  the  projectile.  The  increased  drag  can  result 
in  unacceptable  velocity  drop  over  the  flight  range  and  reduce 
the  effective  test  time  of  the  range.  In  addition,  friction 
and  flow  interaction  phenomena  arising  in  the  projectile/rail 
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interface  region  can  augment  the  drag  (and  thus  the  velocity  drop) 
and  can  threaten  projectile  integrity  by  inducing  large  amounts 
of  projectile  wear. 

There  are  currently  four  hypervelocity  guide  rail  systems 
in  operation;  Tracks  K and  G at  AEDC/  the  McDonnell  Douglas  range 
in  El  Segundo,  California,  and  the  Air  Force  Materials  Laboratory/ 
University  of  Dayton  Research  Institute  (AFML/UDRI)  light  gas 
gun  range  at  Wright-Patterson  Air  Force  Base.  Larger  systems 
are  in  the  preliminary  planning  stages  at  AEDC.  This  report  is 
principally  concerned  with  the  AEDC  ranges;  however,  the  analysis 
has  been  made  sufficiently  general  to  handle  the  scaling  between 
the  various  ranges  currently  in  operation  or  being  planned.  As 
Tracks  K and  G are  currently  operating  and  some  performance  data 
are  available,  typical  parameters  for  these  two  facilities  will 
be  used  for  example  problems. 

1.1.1  Track  K Parameters 

Track  K is  typical  of  current  guide  rail  designs  and 

is  described  in  some  detail  by  Teng,  et  al.^  A drawing  of  the 

system  is  shown  in  Figure  1.  There  are  two  sets  of  opposing 

rails  7.92  mm  wide.  The  inside  of  the  guide  tube  is  accurately 

honed  to  96.52  ram  '*’^9  The  rail  height  is  then  precision  ground 

to  27.43ram  ^ and  the  rails  are  located  on  the  inside  of  the 

- 0 pm 

tube  to  provide  an  inter-rail  gap  of  41.33  mm  ±30  ym.  The  guide 
tubes  are  constructed  in  3.05  m sections  and  ten  of  them  are 
connected  together  to  form  a systCTi  30.5  m long.  Care  was  taken 
to  ensure  proper  inter-section  alignment. 

The  Range  K light  gas  gun  has  a 41.28  mm  diameter  launch 
tube  and  typically  launches  cylinders  with  various  forebodies. 

The  gap  between  the  guide-rails  and  the  projectile  is  nominally 
50  ym*.  The  projectiles  typically  have  a cylindrical  section 


*Recent  in-flight  measurements  of  projectiles  launched  on  the 
Track  K facility  show  projectiles  having  post  launch  diameters 
of  40.64  mm  resulting  in  an  increased  gap  height  of  the  order  of 
350  ym. 
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2 

41. 2B  min  long  giving  a nominal  rail  contact  area  of  327  mm  per 
rail.  The  half-angles  of  the  conical  forebodies  are  typically 
in  the  range  of  45®  to  10®.  For  simulation  of  typical  current 
and  projected  reentry  vehicles,  half-angles  of  10"  are  most 
appropriate.  Therefore,  principal  attention  in  analysis  develop- 
ment was  directed  to  cone  half-angles  between  10°  and  30®. 

Projectiles  are  typically  constructed  from  polycarbonate 
plastic  material.  For  reentry  simulation,  the  cone  section 
surface  must  be  a more  realistic  reentry  vehicle  nose  tip  material 
such  as  graphite.  Typically,  the  cones  are  either  covered  with 
or  fabricated  from  materials  of  interest,  but  the  cylindrical 
portion  which  contacts  the  rails  will  probably  remain  poly- 
carbonate or  some  other  impact  resistant  plastic. 

The  Track  K launcher  typically  launches  these  projectiles 
to  velocities  between  2.5  km/s  and  4.5  km/s  although  reentry 
testing  can  be  extended  up  to  velocities  approaching  6 km/s. 

The  range  atmosphere  is  usually  air  at  densities  between  0.1  and 
1.0  atm  depending  on  the  stage  of  reentry  to  be  simulated. 

The  following  is  a list  of  the  parameters  of  interest 
for  the  Track  K at  AEDC  with  typical  values: 

1.  Projectile: 

Velocity  - 2.5  km/s  to  4.5  km/s 

Material  - Polycarbonate 

Geometry  - 41.28  mm  diameter  cylinder 

10°  to  30°  half-angle  conical  forebody 

2.  Rails: 

Material  - Steel 

Inter-rail  gap  - 41.33  mm 

Length  - 30  m 

3.  Projectile  - Rail  Interface: 

Gap  height  - <0.35  mm  nominal 

2 

Contact  area  ^<300  mm  nominal 

4.  Range  Atmosphere: 

Gas  - air 

Pressure  - 0.1  atm  to  1.0  atm 
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1.1.2  Track  G Parameters 

Track  G is  similar  to  Track  K,  but  considerably  larger. 
The  Range  consists  of  a light  gas  gun  with  a 63.5  mm  diameter 
launcher  tube  and  a nominal  300  ra  long  free  flight,  range  tank. 
The  track  system  is  installed  in  the  range  tank  such  that  it  can 
be  swung  off  trajectory  to  permit  unguided  projectile  testing. 
The  track  and  range  system  is  described  in  detail  by  Norfleet 
et  al.“  Track  G parameters  are  summarized  as  follows: 

1.  Projectile; 

Velocity  - 3 km/s  to  6 km/s 

Material  - Composite  with  polycarbonate 

resin  sabot,  aluminum  or  titanium 
internal  strut  and  heat  shield  back 
up  material,  and  carbon  epoxy  heat 
shield  material 

Geometry  - 63.5  mm  (2.5  inch)  diameter  cylinder 
(nominal) . Various  forebody  designs 

2.  Rails: 

Material  - Steel 

Inter-rail  gap  - 63.5  mm 

Length  - 277.4  m (910  ft).  91  sections 
10  ft/section 

3.  Projectile  - Rail  Interface: 

Gap  height  - ^.50  mm  (0.020  inch)  nominal 

2 

Contact  area  - 500  mm  nominal 

4.  Range  Atmosphere: 

Gas  - air  (can  use  nitrogen,  argon,  helium,  etc) 
Pressure  - 0.00026  atm  to  1.0  atm 

1.2  PROJECTILE-RAIL  SLIDING  PHENOMENA 

This  report  is  particularly  addressed  to  describing 
the  phenomena  associated  with  projectile-rail  sliding  and  the 
effects  of  these  phenomena  on  guide-rail  performance.  Three 
aspects  of  hypervelocity  sliding  will  be  treated?  friction,  wear, 
and  gouging. 
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Friction  occurs  when  any  two  solid  surfaces  in  close 
proximity  slide  with  respect  to  one  another.  The  frictional 
force  opposes  the  slider  (projectile)  motion  and  so  adds  to  the 
aerodynamic  drag  force  on  the  projectile.  This  effect  should 
be  small  for  a successful  guide-rail  system. 

Wear  is  closely  associated  with  the  frictional  process 
and  is  defined  as  the  mass  removed  from  the  rail  and/or  pro- 
jectile by  the  sliding  process.  Wear  is  important  in  guide  rail 
systems  for  three  reasons: 

1.  Wear  increases  the  projectile-rail  gap  separation 
and  thus  degrades  the  position  and  orientation  con- 
trol of  the  system. 

2.  Wear  could  remove  sufficient  material  to  appreciably 
lower  the  available  testing  area  of  the  projectile. 

3.  Wear  could  remove  sufficient  material  to  threaten 
projectile  integrity  either  by  parting  the  pro- 
jectile or  by  providing  sufficient  space  for 
lateral  or  rotational  oscillations  of  the  pro- 

• jectile  to  grow  and  ultimately  disintegrate  the 

projectile. 

As  wear,  no  doubt,  increases  with  the  distance  traveled,  all  of 
these  problems  could  become  very  important  on  the  very  long 
ranges  under  consideration.  Indeed,  wear  has  been  found  to  be 
an  important  constraint  in  Track  G.  Wear  on  Track  K models  has, 
on  the  other  hand,  been  minimal. 

Gouging  is  a phenomenon  peculiar  to  hypervelocity  sliding 
contacts.®’**’®  Under  certain  conditions  tear  shaped  craters 
are  formed  on  the  rail  (and  presumably  the  projectile)  surface 
during  sliding.  Gouging,  therefore,  is  a relatively  violent 
process  and  has  two  possible  detrimental  effects: 

1.  The  rails  can  become  so  damaged  that  they 
must  be  refinished  or  replaced  periodically. 

This  will  be  expensive  and  time  consuming 
when  very  large  and  long  systems  are  developed. 

2.  The  violence  of  the  gouging  process  could  re- 
sult in  disintegration  of  the  projectile  and  com- 
plete loss  of  data.  Subsequent  rail  refinishing 
or  replacement  would  again  be  required. 
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This  investigation  is  addressed  to  understanding  friction, 
wear,  and  gouging  phenomena  sufficiently  to  design  long  rail 
systems  and/or  to  scale  the  behavior  of  current  systems  to  ones 
developed  in  the  future. 

1 . 3 PROJECTILE-RAIL  DYNAMICS 

The  behavior  of  the  interface  region  between  the  projectile 
and  the  rail  is  a function  of  the  dynamics  of  the  projectile  in 
the  rail  system.  In  particular,  the  projectile  can  undergo 
relative  lateral  translation  and  rotation  due  to  aerodynamic 
forces  and  moments  and  the  curvature  of  the  rails.  These  pro- 
jectile motions  result  in  changes  in  the  interface  behavior  which, 
in  turn,  affects  system  dynamic  behavior  performance.  It  is, 
therefore,  necessary  to  incorporate  projectile  and  coupled  in- 
terface dynamics  into  any  realistic  description  of  guide-rail 
performance. 

This  report  identifies  the  phenomena  of  importance  in 
hypervelocity  sliding.  These  phenomena  are  then  analyzed  in 
detail  and  modeled  for  the  particular  conditions  of  interest  in 
guide-rail  systems.  The  dynamics  of  rail  guided  projectiles  are 
then  analyzed  and  modeled.  Finally,  the  interface  model  and 
projectile  dynamics  model  are  combined  to  provide  an  overall 
model  of  the  dynamic  behavior  of  hypervelocity  guide-rail 
sliding  performance.  The  results  obtained  from  exercising 
the  combined  progreun  are  compared  with  the  experimental  data 
currently  available. 
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SECTION  II 

BACKGROUND  REVIEW  OF  SLIDING  PHENOMENA 


The  following  paragraphs  outline  background  information 
on  the  sliding  of  solids  and  describe  the  conceptual  framework 
in  which  the  present  problem  will  be  investigated. 

2.1  LOW  AND  MEDIUM  VELOCITY  SLIDING 

Low  to  medium  velocity  friction  has  been  studied  exten- 
sively and  is  described  in  a book  by  Bowden  and  Tabor.  It  is 
instructive  to  review  briefly  these  well-understood  mechanisms 
before  considering  phenomena  occurring  during  hypervelocity 
sliding. 

Static  friction  between  solids  is  attributed  to  the 
cold  welding  of  asperities  on  the  surfaces  as  they  are  brought 
into  contact.  A tangential  force  applied  to  the  interface  must 
shear  these  welded  regions r thus  giving  rise  to  the  frictional 
force.  The  area  of  welded  contact  is  proportional  to  the 
normal  force  and  the  shear,  or  friction  force  is  proportional 
to  the  area  of  welded  contact.  Therefore,  the  coefficient  of 
friction  which  is  defined  as  the  ratio  of  the  maximum  frictional 
force  to  the  normal  force,  is  predicted  to  be  independent  of  the 
normal  force,  as  is  observed  experimentally.  The  process  is  a 
material  shear  strength  dominated  phenomenon. 

When  motion  takes  place,  the  welded  areas  shear,  new 
ones  form,  and  the  behavior  remains  essentially  the  same  as  for 
static  friction. 

As  sliding  velocities  increase,  there  is  insufficient 
time  for  the  heat  dissipated  during  shearing  to  be  conducted  into 
the  bulk  material.  Thermal  eguilibrium  cannot  be  established, 
and  the  shear  region  heats  up.  The  rise  in  temperature  at  the 
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interface  results  in  a lowering  of  the  shear  strength  of  the 
material  in  the  cold  welded  areas  and  friction  forces  decline. 
Typically,  the  coefficient  of  friction  declines  as  shown  in 
Figure  2.  Finally,  the  entire  interface  becomes  covered  with 
a molten  layer  and  the  friction  becomes  essentially  fluid 
dynamic  and  dependent  upon:  (a)  the  viscosity  of  the  molten 
interface  materials;  (b)  sliding  velocity;  and  (c)  interfacial 
thickness.  Therefore,  over  some  range  of  velocity,  the  friction 
process  transitions  from  strength  dominated  to  viscous  dominated. 

The  viscous  force,  F,  which  takes  the  place  of  the 
classic  friction  force  once  surface  melting  occurs,  may  be 
expressed  in  the  form; 


P = llA  g , (1) 


where  y is  the  coefficient  of  viscosity,  A is  the  area  of  the 
interface,  and  dv/dy  is  the  transverse  velocity  gradient  in  the 


Figure  2.  The  Coefficient  of  Friction,  y,  for  Steel  on  Nylon 
as  a Function  of  Sliding  Velocity.  (Taken  from 
Bowden  and  Tabor, ^ Pg.  476). 
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interface  zone  (i.e.r  transverse  to  the  sliding  direction}.  The 
average  velocity  gradient  between  surfaces  sliding  at  velocity  v 
and  separated  by  distance  h is; 


(dv/dy)^y  = v/h. 


(2) 


Note  that  the  separation,  h,  is  the  distance  between  solid  pro- 
jectile and. rail  surface.  Substituting  this  expression  into  the 
previous  equation,  we  arrive  at  an  "average"  viscous  force, 


_ y A V 
^av  " h 


(3) 


As  shown  in  Equation  (3) , the  viscous  force  varies  directly 
with  the  viscosity  and  velocity  and  inversely  with  separation 
distance. 

It  should  be  noted  here  that,  for  a fluid  dynamic 
interface,  the  frictional  force  will  have  a strong  velocity 
dependence.  Such  frictional  behavior  maJces  it  difficult  to 
distinguish  between  sliding  "friction"  and  aerodynamic  drag  in 
projectile  deceleration  data  acquired  from  guide-rail  systems. 

As  the  aerodynamic  drag  is  also  velocity  dependent,  the  func- 
tional dependences  of  the  "friction"  and  the  aerodynamic  drag 
may  well  be  nearly  identical. 

Further  increases  in  sliding  velocity  result  in  diminish' 
ing  thermal  equilibrium  in  the  interface  and  the  thickhess,  h, 
of  the  interface  zone  contracts,  resulting  in  higher  velocity 
gradients  within  the  fluid  layer  and  increased  drag.  The 
available  experimental  data  on  high  speed  contacts  ends  at  about 
700  m/s,  well  below  the  velocities  at  which  this  investigation 
begins . 
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2.2  HYPERVELOCITY  SLIDING 

As  velocity  increases  from  zero  the  friction  mechanism 
quite  clearly  changes  from  a material  strength  dominated 
phenomena  to  a viscous  dominated  process  and  hypervelocity 
friction  is  almost  certainly  a completely  viscous  process. 

A logical  extension  to  the  previous  discussion  appears 
to  be  that,  as  velocity  is  further  increased,  the  energy  dissipa- 
tion rate  in  the  interface  becomes  sufficient  to  vaporize  signi- 
ficant amounts  of  the  interfacial  material.  In  addition,  for  a 
guide-rail  system,  ambient  air  is  entrained  by  the  projectile 
motion  and  injected  into  the  interface  forming  a self-pumped 
gas  bearing. 

The  gas  velocity  varies  across  the  interface  from  zero 
at  the  rail  surface  to  the  sliding  velocity  at  the  projectile 
surface.  The  average  velocity  of  the  interfacial  material  is, 
therefore,  approximately  one-half  the  sliding  velocity. 

Material  originating  at  the  projectile  or  rail  surface  and 
moving  into  the  interfacial  layer  must  convert,  on  the  average, 
3/4  of  its  original  kinetic  energy  to  other  forms  of  energy. 

This  energy  becomes  available  for  chemical  decomposition  of  the 
material,  material  heating,  and  phase  changes  (i.e.,  melting, 
boiling,  or  sublimation) . 

If  we  neglect  energy  described  in  phase  changes  and 
decomposition  and  consider  that  all  the  available  kinetic  energy 
is  converted  to  thermal  energy  we  obtain  for  a monatomic  gas; 


3nmUf, 

CO 

“8 


_ 3nKr 


(4) 


where  n is  the  number  density  of  atoms,  m is  the  molecular  mass, 
K is  Boltzman  constant,  is  the  sliding  velocity  and  T is  the 
gas  temperature.  Solving  for  T we  obtain; 


T 


oo 

TT~  ' 


(5) 
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Solutions  for  Equation  (5)  are  shown  in  Table  1 for  various  atoms 
that  might  be  present  in  the  interface.  Larger  polyatomic 
molecules  have  larger  masses  (and  kinetic  energy)  but  they  also 
have  more  degrees  of  freedom  in  which  to  absorb  energy.  The 
temperatures  for  polyatomic  molecules  are  therefore  probably 
lower  for  single  atom  molecules,  but  the  same  order  of  magnitude 
as  those  displayed  in  Table  1. 

TABLE  1.  GAS  TEMPERATURES  VS.  SLIDING  VELOCITY  FOR  OXYGEN, 
CARBON,  NITROGEN 


GAS  TEMPERATURE,  “K 


U m/s 

0 (16) 

C (12) 

N (14) 

1000 

484 

363 

423 

2000 

1936 

1452 

1694 

3000 

4356 

3267 

3811 

4000 

7744 

5808 

6776 

5000 

12101 

9076 

10588 

6000 

17426 

13069 

15247 

7000 

23718 

17789 

20753 

The  energy  available  to  the  interfacial  zone  varies  as  the 
square  of  the  sliding  velocity.  Therefore,  vaporization  might 
be  expected  at  moderate  velocities,  perhaps  of  the  order  of  a 
few  kilometers  per  second  for  the  materials  currently  used  for 
hypervelocity  projectiles,  and  the  interface  will  display  gas 
dynamic  behavior.  It  should  also  be  noted  that  many  plastics 
such  as  Teflon  sublime  rather  than  melt  when  heated  and  the 
frictional  behavior  of  these  materials  might  not  display  the 
hydrodynamic  phase  but  transfer  directly  to  the  gas  dynamic 
phase  in  a similar  manner  to  a recently  proposed  solid  CO2 
bearing . ’ 

Frictional  forces  might  be  expected  to  be  extremely 
small  for  an  ablative  gas  dynamic  bearing,  but  will  increase 
with  velocity  because  the  interface  thickness  will  contract 


21 


AEDC-TR.78-19 


(leading  to  higher  transverse  velocity  gradients)  due  to 
increasing  thermal  disequilibrium  and  its  temperature  will  rise 
resulting  in  increasing  gas  viscosity.  For  hypervelocity 
sliding,  therefore,  a gaseous  interface  will  be  considered 
except  when  normal  loads  become  sufficient  to  reduce  the  inter- 
face to  zero  thickness  causing  direct  contact  between  the  pro- 
jectile and  the  rail.  Such  direct  contact  results  in  extremely 
high  transverse  velocity  gradients.  The  high  energy  dissipa- 
tion rates  associated  with  these  gradients  could  account  for 
the  crater-like  surface  damage  characteristics  of  gouging.  This 
conceptual  approach  to  gouging  is  supported  by  the  work  of 
Graff  and  Dettloff^  who  studied  high  velocity  gouging  observed 
during  operation  of  rocket  sled  tracks.  They  found  that 
gouging  occurred  only  when  a critical  normal  force  between  a 
projectile  and  a rail  was  exceeded  and  that  gouging  could  be 
inhibited  by  coating  the  rail  with  a variety  of  soft  materials. 
We  feel  that  both  reducing  the  normal  forces  and  introducing  a 
soft  pressure  bearing  but  essentially  fluid  layer  serve  to 
maintain  separation  between  the  solid  elements  and  keep  the 
transverse  velocity  gradients  sufficiently  low  to  prevent 
gouging . 
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SECTION  III 

SLIDING  IN  A HYPERVELOCITY  RAIL  SYSTEM 

3 . 1 INTRODUCTION 

As  was  mentioned  in  Section  II,  sliding  in  a hypervelocity 
rail  system  is  basically  a fluid  dynamic  phenomenon  except  when 
critical  normal  loads  become  sufficient  to  reduce  the  interface 
to  nearly  zero  thickness.  The  minimum  thickness  below  which  the 
interface  can  no  longer  be  considered  to  be  purely  fluid  dynamic 
is  on  the  order  of  the  surface  roughness  of  either  the  rail  or 
the  projectile  surface,  whichever  is  larger.  The  critical  normal 
load  required  to  reduce  the  interface  to  a thickness  on  the  order 
of  the  dimensions  of  surface  asperities  can  be  expected  to  in- 
crease as  the  projectile  velocity  increases  because  the  viscous 
dissipation  will  increase  and  therefore  the  ablation  rate  will 
increase.  The  presence  in  the  interface  of  gases  from  ablated 
material  makes  an  additional  contribution  to  the  pressure  of 
the  gas  in  the  interface  over  and  above  that  due  to  fluid  dy- 
namic effects. 

The  interface  thickness,  which  will  be  called  gap  height, 
is  the  critical  parameter  in  a discussion  of  hypervelocity  rail 
systems.  Various  flow  regimes  can  be  hypothesized  depending  on 
the  gap  height.  A convenient  gauge  for  determining  the  type  of 
flow  regime  which  can  be  expected  is  the  lubrication  Reynolds 
number  Re*  defined  as; 


of  which  p is  the  gas  density,  B is  the  length  of  the  projectile 
cylindrical  surface,  is  the  projectile  velocity,  y is  the  gas 
viscosity,  and  h is  the  gap  height.  If  Re*<<l,  the  interfacial 
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flow  is  completely  dominated  by  viscous  forces  and  the  inertial 
effects  of  the  gas  flow  are  negligible.®  In  this  case,  the  in- 
terface flow  can  be  treated  as  a purely  laminar  flow  in  which 
only  the  viscous  terms  in  the  governing  equations  (Navier  Stokes 
equations)  are  considered.  If  Re*  is  of  order  1,  then  the  flow 
is  still  laminar,  but  inertia  terms  must  be  included.  If  Re*>>l, 
the  flow  is  turbulent. 

Thus,  determining  the  interface  forces  acting  on  a projectile 
in  a hypervelocity  guide-rail  system  is  a fluid  dynamic  problem 
except  when  normal  forces  are  greater  than  the  critical  normal 
force  required  to  produce  solid-solid  contact.  In  this  section 
we  present  and  discuss  the  methods  and  procedures  which  have 
been  used  to  determine  the  fluid  dynamic  forces  on  the  projectile 
and  to  determine  the  wear  of  the  projectile  due  to  ablation  caused 
by  aerodynamic  heating.  The  special  forms  of  the  governing 
equations  used  to  model  the  various  flow  regimes  are  discussed 
along  with  the  boundary  conditions  to  be  applied  to  the  governing 
equations.  The  methods  and  procedures  used  to  describe  the 
material  properties  including  thermodynamic  and  transport  properties 
of  the  gas  are  also  presented. 

3 . 2 GOVERNING  EQUATIONS 

The  Navier  Stokes  (N-S)  equations  of  fluid  dynamics 
are  the  basic  equations  to  be  applied.  The  full  N-S  equations 
will  not  be  written  out  here  since  they  can  be  found  in  any 
standard  fluid  dynamics  textbook  (see  Chapter  12  of  Reference  9 
for  example) . The  special  forms  of  the  N-S  equations  applicable 
to  the  three  different  flow  regimes  listed  in  the  previous 
paragraphs  will,  however,  be  given  since  they  form  the  basic 
elements  of  the  fluid  dynamic  models  developed  under  this  effort. 

3.2.1  Coordinate  System 

Before  discussing  the  governing  equations  themselves, 
the  coordinate  system  in  which  the  fluid  dynamic  equations  are  to 
be  solved  will  first  be  developed.  First,  it  should  be  emphasized 
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that  the  overall  objective  of  this  effort  is  to  determine  the 
wear  history  of  a projectile  as  it  moves  down  the  rail  system. 

The  forces  required  to  determine  the  trajectory  of  the  pro- 
jectile are  basically  transient  in  nature  since  the  basic  fluid 
dynamics  is  transient.  In  order  to  be  able  to  develop  relatively 
simple  fluid  dynamic  models  for  determining  the  forces  on  the 
projectile,  it  was  necessary  to  perform  a quasi-steady  analysis 
of  the  flow  field.  That  is,  it  was  assumed  that  the  flow  adjusts 
rapidly  enough  to  changes  in  angle-of-attack  and  position  of 
center— of-gravity  that  a steady  flow  solution  to  an  instantaneous 
geometry  would  adequately  describe  the  instantaneous  forces  on 
the  projectile.  From  the  point  of  view  of  an  observer  in  a 
coordinate  system  fixed  to  the  rails,  the  flow  is  always  transient. 
However,  an  observer  riding  with  the  projectile  sees  a relatively 
steady  flow  situation  as  long  as  the  aerodynamic  blockage  in  the 
system  is  not  too  great.*  Thus,  the  coordinate  system  used  in 
this  work  is  one  attached  to  the  projectile  and  it  is  assumed 
that  the  blockage  ratio  is  greater  than  3.  In  this  coordinate 
system  an  observer  riding  with  the  projectile  sees  the  range 
gas  and  the  rails  moving  toward  him  with  a velocity  that  is  equal 
and  opposite  to  the  absolute  velocity  of  the  projectile.  The 
basic  flow  situation  is  sketched  in  Figure  3.  For  the  interface 
flow  analysis  the  coordinate  system  origin  is  taken  as  shown  in 
Figure  4 . 

Some  details  of  the  shock  wave  structure  due  to  the 
presence  of  the  rails  are  shown  in  Figure  3.  Note  particularly 
that  boundary  layer  buildup  begins  on  the  rail  at  the  point  where 
the  bow  shock  intersects  the  rail. 


*In  Reference  1,  it  is  suggested  that  if  the  open  area  of  the 
range  tube  divided  by  the  closed  area  is  greater  than  three,  a 
normal  shock  will  not  move  down  the  range  tube  ahead  of  the  pro- 
jectile and  the  flow  relative  to  the  projectile  will  be  steady. 
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Figure  3.  Boundary  Layer  Build-up. 
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3.2.2  The  Reynolds  Equations 

By  assuming  steady  flow  in  the  interface  and  neglecting 
the  inertial  terms,  body  forces,  z-direction  components  and 
pressure  gradient,  the  classical  momentum  equations  for  the  in- 
terface flow  reduce  to; 

= _L 

3x 

and; 

3P  _ 3 
"Sy  3z 

in  which  P=P  (x,y) , u=u  (x,y,z)  and  v=v  (x,y,z).  In  these 
equations,  P is  the  pressure,  u is  local  gas  velocity  in  the  x 
direction,  v is  the  local  gas  velocity  in  the  y direction, 
is  the  viscous  shear  force  acting  in  the  x-direction  on  a plane 
perpendicular  to  the  z-direction  and  is  the  viscous  shear 
force  acting  in  the  y-direction  on  a plane  perpendicular  to  the 
z-direction. 

The  continuity  (conservation  of  mass)  equation  is; 


(•'ll.)  = ^ (s.)  - 
(''It)'  -k  (Vz)  ' 


_3_ 

3x 


(ou)  + (pv)  + (pw) 


= 0, 


(9) 


where  the  terms  u,  v,  w are- velocities  in  the  x,y,z  directions, 
respectively . 

Integrating  with  respect  to  z,  one  obtains; 


z=h 


pw 


z=0 


(pu)  + 


3 

3y 


■ 


(pv)  dz. 


(10) 


By  interchanging  the  order  of  integration  and  differentiation 
the  following  equation  is  obtained; 
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pw 


1=0  - - 0" 


pv  dz. 


(11) 


The  left  hand  side  of  the  above  equation  is  zero  if  there  is 
no  ablation  since  the  w-component  of  velocity  vanishes  at  solid 
surfaces. 

Integration  of  Equation  (7)  yields; 


3P  3u  , . 

3^  ^ °1 


iH  - i li  T - fl^*^ 

3z  ~ y 3x  y ' 


z z 

u =f  — -1^  z dz  - c-  /"  — dz  + U , 

./q  0 ^ 


(12) 


where  Cj^(x)  is  an  integration  constant  and  the  no-slip  condition, 
u=U^  at  z=0  has  been  applied. 

For  the  simplified  case  where  the  viscosity  is  taken 
to  be  a constant  (or  replaced  by  a suitably  defined  average 
value)  the  following  equation  for  u is  obtained; 


(13) 


For  the  case  of  constant  viscosity,  integration  of 
Equation  (8)  yields; 


hz  (3P/3y)  z\ 

^ 2y  ~ h ) ’ 


(14) 


The  velocity  field  for  constant  fluid  viscosity  is  thus  a 
combination  of  very  elementary  viscous  flows.  Equation  (13) 
represents  the  solutions  to  laminar  Couette  flow^  with  pressure 
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gradient  (Petrov-Couette  flow’)  and  Equation  (14)  represents 
plane  Poiseuille’  flow.  The  velocity  field  obtained  by  super- 
position of  these  two  elementary  flows  is  assumed  to  describe 
the  basic  interface  flow  for  the  case  of  very  small  gap  heights 
(Re*<<l) . The  primary  objective  here,  however,  is  not  necessarily 
to  define  precisely  the  velocity  field,  but  rather  to  determine 
reasonable  approximations  to  the  normal  and  shear  force  distri- 
butions on  the  projectile.  The  pressure  distribution  on  the  pro- 
jectile could  be  obtained  by  combining  Ecjuations  (11),  (12),  and 
(13) , yielding  the  well-known  Reynolds  Equation  of  lubrication 
theory; 


+ 12  pW. 


(15) 


The  term  W is  the  projectile  velocity  component  normal  to  the 
rail.  The  solution  to  this  equation  yields  the  pressure  dis- 
tribution. Integration  of  the  pressure  distribution  over  the 
interface  region  then  determines  the  total  normal  force  on  the 
projectile  in  the  interface.  The  viscous  shear  forces  on  the 

projectile  could  then  be  determined  by  solving  for  t (z=h) 

xz 

and  Ty^(z=h)  in  Equations  (7)  and  (8)  respectively. 

Application  of  Equation  (15)  to  determining  the  pressure 
distribution  and  shear  force  distribution  on  the  projectile  over 
the  interface  region  for  hypervelocity'  guide  rail  systems  is  a 
formidable  task.  The  density  and  viscosity  in  Equation  (15)  are 
highly  variable  in  the  z-direction  because  of  the  large  temper- 
ature variation  through  the  interface  gap.  In  order  to  determine 
the  temperature  variations  through  the  interface  gap,  the  energy 
equation  must  be  coupled  to  Equation  (15) . The  energy  equation 
for  the  case  of  Re*«l  is  given  by;® 
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dX 
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(16) 
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The  term  k in  Equation  (16)  is  the  thermal  conductivity  of  the 
gas  in  the  interface. 

Thus,  a highly  non-linear  strongly  coupled  system  of 
partial  differential  equations  would  have  to  be  solved,  in  order 
to  determine  the  pressure  and  shear  force  distribution  over  the 
projectile.  This  could  only  be  done  numerically  and  would 
require  far  more  computer  time,  even  for  a simple  geometry,  than 
would  be  practical  for  the  present  task. 

Therefore,  approximate  solutions  to  Equations  (15)  and 
(16)  available  in  the  literature  were  used  to  establish  order- 
of-magnitude  estimates  of  the  pressure  and  shear  forces  in  the 
interface  region.  These  approximate  solutions  were  primarily 
those  obtained  for  compressible  Couette  flow.  The  rationale 
for  making  extensive  use  of  compressible  Couette  flow  solutions 
(for  both  laminar  and  turbulent  cases)  in  the  present  effort  is 
discussed  later  in  Section  3.5  after  the  governing  equations  for 
the  Re*!*!  and  Re*>>l  situations  have  been  presented  and  after 
the  boundary  conditions  have  been  discussed. 


3.2.3  Governing  Equations  for  Re*=l 

The  Reynolds  Equation  (15)  is  based  on  purely  viscous 
flow.  If  the  gap  height  is  large  enough  for  Re*  to  be  of  order 
one,  the  fluid  inertia  terms  are  important  and  must  be  included. 
The  addition  of  fluid  inertia  terms  to  Equations  (7)  and  (8) 
yields  the  following  equations; 


and; 


^ / 3u  3u  \ 3P  . 3 

'^3y  ) ” 3x  3z 


(17) 


(18) 
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Solutions  to  Equations  (17)  and  (18)  available  in  the  gas  bearing 
lubrication  literature  show  that  inertia  effects  are  generally 
small  in  terms  of  influencing  the  forces  on  gas  bearing  surfaces 
when  Re*  is  of  order  one,  but  become  increasingly  important  as 
Re*>>l.  However,  when  Re*  becomes  much  larger  than  1,  turbulence 
effects  become  the  primary  inertia  effects  and  so  the  solution 
to  Equations  (17)  and  (18)  are  useful  to  the  present  problem 
only  over  a very  limited  range  of  gap  heights.  Thus,  no  effort 
was  made  to  develop  an  interface  flow  model  which  included 
laminar  flow  with  inertia  terms  included. 

3.2.4  Governing  Equations  for  Re*>>l 

The  applicable  equations  for  the  flow  in  the  interface 
when  Re*  is  much  larger  than  one  are  the  turbulent  boundary  layer 
equations  for  compressible  flow.  Because  the  ratio  h/B  is  a 
very  small  number  in  the  hypervelocity  guide  rail  system,  it  is 
appropriate  to  use  a mixing  length  theory  to  relate  the  Reynolds 
stress  to  density  and  mean  velocity  gradient.  The  turbulent 
boundary  layer  equations  for  three  dimensional  boundary  layers 
are; 


^ (P")+  ^ ^ (pw)  = 0' 


(19) 


{ 3u 

/ 3v 
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4.  A 
3x  3z 
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[ (-pv) ' w' ] 


(20) 

(21) 


Quantities  primed  represent  turbulent  fluctuations  of  the 
specified  quantity.  Following  the  work  of  Korkegi,  et  al., 
the  Von  Karman  mixing  length  model  is  used  to  relate  the  Reynolds 
stress  terms  to  mean  velocity  gradients.  The  last  term  in 
Equations  (20)  and  (21)  are  then  respectively; 
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3u{  3u 

3z|  dz  ' 

(22) 

3wj  3w 

dz|  3z  * 

(23) 

These  equations,  like  the  Reynolds  Equations,  are  much  too 
difficult  to  solve  as  a part  of  the  interface  flow  model  and  re- 
course is  made  to  using  approximate  solutions  available  in  the 
literature.  In  particular,  extensive  use  has  been  made  of  the 
work  of  Korkegi  and  Briggs*®’  on  hypersonic  slipper  bearings 

and  turbulent  Couette  flow.  The  hypersonic  slipper  bearing 
analyses  of  Korkegi  and  Briggs  as  well  as  the  general  literature 
on  gas  bearing  analyses  is  not  directly  adaptable  to  the  present 
problem  because  of  differences  in  boundary  conditions. 

3.3  BOUNDARY  CONDITIONS  FOR  INTERFACE  REGION 

In  order  to  completely  define  an  interface  flow  model, 
the  boundary  conditions  to  be  applied  to  the  interface  must  be 
specified. 

3.3.1  Entrance  Conditions 

The  conditions  of  the  flow  just  as  it  enters  the  interface 
are  primarily  determined  by  the  aerodynamics  of  the  nose  cone 
guide-rail  system.  As  is  depicted  in  Figure  5,  the  conical  shock 
attached  to  the  cone  tip  interacts  with  the  rail  to  produce  a 
complicated  shock-expansion  wave  system  imbedded  in  the  cone 
flow  field.  The  velocity  profile  at  the  entrance  to  the  inter- 
face is  basically  three-dimensional  and  is  influenced  by  the 
boundary  layer  buildup  both  on  the  cone  and  on  the  rail  section 
preceding  the  interface.  Some  of  the  details  of  how  the  inter- 
face entrance  conditions  are  defined  are  given  in  the  following 
paragraphs.  The  other  boundary  conditions  are  discussed  in 
succeeding  paragraphs . 
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Figure  5.  Projectile-Rail  Geometry. 


The  initial  effort  centered  on  determining  the  entrance 
conditions  for  30“  and  45“  half -angle  sharp  cones  because  ex- 
perimental data  were  available  to  guide  model  development.  Care 
was  taken  to  assure  the  generality  of  the  approach  so  that 
application  to  lower  cone  angles  of  interest  (near  10“  half-angle) 
can  be  readily  accomplished.  The  laser  shadowgraphs  presented 
in  Reference  1 for  30“  and  45“  half-angle  cones  at  a Mach  num- 
ber of  about  13  show  that  normal  shocks  (or  possibly .strong 
oblique  shocks)  are  set  up  between  the  rail  and  the  cone.  The 
model  developed  herein  predicts  the  conditions  under  which  the 
normal  (or  strong  oblique)  shocks  will  occur.  Downstream  of 
the  normal  or  strong  oblique  shock,  the  flow  is  subsonic;  but 
the  pressure  is  higher  than  that  experienced  by  the  flow  away 
from  the  rail  which  has  undergone  only  compression  by  the 
conical  shock  attached  to  the  cone.  Therefore,  that  part  of  the 
flow  downstream  of  the  strong  shock  will  expand  and  accelerate 
until  the'  pressure  of  the  gas  which  underwent  the  strong  shock 
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reaches  the  pressure  of  the  gas  which  was  compressed  by  the 
weak  conical  shock.  The  condition  of  the  gas  which  first  under- 
went the  strong  shock  and  then  expanded  back  to  the  pressure 
corresponding  to  the  pressure  behind  the  attached  conical  shock, 
is  assumed  to  be  approximately  the  conditions  at  the  entrance 
to  the  rail-projectile  interface. 

In  actuality,  the  gas  flow  which  is  influenced  by  the 
presence  of  the  rail  is  highly  three-dimensional  since  the  in- 
teraction of  the  conical  shock  with  the  rail  is  a local  phenomenon 
imbedded  in  the  conical  flow  field.  In  Figure  6,  some  details  of 


Figure  6.  Shock  System  on  30°  and  45“  Cones  Showing  Definition 
of  Flow  Regions. 

the  conical  shock-rail  interaction  are  shown.  Point  3 refers 
to  the  gas  which  has  been  compressed  to  a pressure  over  and 
above  the  gas  passing  through  the  attached  conical  shock. 
Surrounding  this  high  pressure  region  is  the  lower  pressure 
region  of  gas,  which  is  not  influenced  by  the  presence  of  the 
rail.  Point  2 refers  to  the  state  of  the  gas  which  only 
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experiences  compression  by  the  conical  shock.  The  gas  at  Point 
3,  therefore,  will  tend  to  expand  to  the  lower  pressure  region 
of  Point  2 where  gas  pressure  is  P2.  This  expansion  process  is 
assumed  to  be  an  isentropic  expansion  to  the  lower  pressure. 

The  minimum  pressure  to  which  gas  at  Point  3 can  expand  is  the 
pressure  corresponding  to  Point  2.  Whether  or  not  there  is 
sufficient  distance  between  Point  3 and  the  entrance  to  the  in- 
terface for  the  gas  at  Point  3 to  expand  completely  down  to  P2 
is  not  certain.  For  30®  and  45®  half -angle  cones,  it  is  felt  that 
this  assumption  is  fairly  realistic.  Under  any  circumstance, 
the  pressure  at  Point  2 is  the  minimum  pressure  which  can  be  ex- 
pected at  the  entrance  to  the  interface  region.  Since  the  flow 
field  downstream  of  a conical  shock  is  not  a uniform  flow  field 
(all  properties  are  assumed  constant  along  ray  lines  except  for 
entropy  which  is  constant  throughout  the  conical  flow  field  down- 
stream of  the  shock) , a consideration  must  be  made  of  the  fact 
that  P2  varies  from  one  ray  line  to  another  in  this  conical  flow 
field.  At  the  high  Mach  numbers  considered  here  (M>10)  nearly 
all  of  the  pressure  increase  occurs  across  the  shock  itself  with 
little  extra  pressure  increase  contributed  by  the  turning  of 
the  flow  behind  the  conical  shock.  Thus,  P2  can  be  considered 
to  be  constant  in  the  conical  flow  region  and  can  be  taken  to  be 
either  the  pressure  just  downstream  of  the  conical  shock  or  to 
be  slightly  higher  pressure  adjacent  to  the  cone  surface. 

During  the  early  part  of  this  investigation,  it  was 
thought  that  a system  of  reflected  shocks  might  be  set  up  be- 
tween the  rail  and  the  projectile  much  like  the  system  of  re- 
flected shocks  set  up  in  a two-dimensional  supersonic  diffuser 
formed  by  one  parallel  wall  and  one  sloping  wall  forming  a con- 
vergent 2-D  passage.  If  there  could  be  a second  reflection 
(reflection  at  the  cone  surface  of  the  shock  reflected  from  the 
rail) , then  this  reflected  shock  would  be  highly  curved  since  it 
would  be  influencing  gas  which  had  already  started  to  expand. 
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For  the  30®  half-angle  cone  traveling  between  Mach  10 
and  15,  it  has  been  established  that  the  second  shock  in  the 
system  (the  first  reflected  shock  from  the  rail)  cannot  undergo 
regular  reflection  at  the  cone  surface  as  a weak  oblique  shock, 
but  must  be  a strong  oblique  shock.  Actually,  this  shock  must 
have  some  curvature  since  the  Mach  number  distribution  just  up- 
stream of  the  second  shock  is  non-uniform.  Thus,  for  a 30®  half- 
angle cone,  the  reflected  shock  at  the  rail  is  assumed  to  be  a 
strong  oblique  shock.  For  the  45®  half-angle  cone,  regular  re- 
flection of  the  conical  shock  itself  is  impossible  in  the  Mach 
number  range  of  10  to  15.  In  this  case,  a normal  shock  is 
assumed  for  the  flow  near  the  rail. 

Thus,  for  both  30°  and  45®  half-angle  cones,  a second 
reflected  shock  is  not  possible  and  the  entire  flow  from  Point 
3 to  the  interface  is  expanding.  In  Section  5.3  it  is  shown 
that  the  strong  shock  assumption  for  the  first  (and  only)  re- 
flected shock  is  indeed  a very  good  assumption  for  30®  half- 
angle cones  on  the  basis  of  total  drag  observations  for  30®  half- 
angle cones  fired  in  Track  K. 

For  smaller  half-angle  cones,  regular  reflection  of  the 
rail  reflected  shock  impinging  on  the  cone  is  theoretically 
possible  although  a reflected  shock  from  the  cone  surface  was 
not  observed  in  the  wind  tunnel  tests  on  10®  cones  described  in 
Reference  13.  For  a 10®  sharp  cone,  the  reflected  shock  off 
the  rail  is  undoubtedly  a weak  nearly  plane  shock.  Between 
cone  angles  of  10®  and  30®  a rather  dramatic  change  in  the  char- 
acter of  the  reflected  shock  must  take  place.  Figure  32  of 
Reference  1 shows  a laser  shadowgraph  of  a 30®  cone.  The  authors 
of  Reference  1 indicate  on  this  figure  the  presence  of  a normal 
shock  although  it  is  not  clear  whether  they  are  referring  to  a 
Mach  reflection  at  the  rail  or  that  the  reflected  shock  is 
nearly  normal  (a  strong  oblique  shock)  as  we  have  assumed. 

At  some  value  of  cone  angle  between  10  and  30  degrees, 
the  gas  dynamic  conditions  produced  by  shock  reflections  off  the 
rail  must  change  from  those  described  by  a weak  shock  solution 
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to  those  of  a strong  shock  solution.  The  pressure  at  Point  3 
(and  thus  the  shock  interference  drag  on  the  cone  due  to  the 
presence  of  the  rails)  will  be  strongly  dependent  on  how  this 
transition  in  character  of  the  reflected  shock  occurred.  It 
is  suspected  that  this  transition  in  character  might  critically 
depend  on  the  projectile  Mach  number  as  well  as  cone  angle.  A 
possible  mechanism  for  the  transition  might  be  a Mach  reflection 
of  the  rail-reflected  wave  at  the  cone  surface,  although  neither 
this  mechanism  nor  any  other  has  been  investigated  as  a part 
of  this  study.  It  might  also  be  possible  that  there  is  actually 
a reflected  wave  from  the  cone  surface  for  10“  cones,  but  it 
might  not  have  been  possible  to  observe  the  reflected  shock  in 
the  Schlieren  photographs  of  Reference  13  due  to  shock  spreading 
and  shock-boundary  layer  interaction,  or  the  quality  of  the 
photographic  reproduction . 

A computer  program  was  written  which  determines  the  flow 
conditions  behind  a normal  shock,  a series  of  oblique  shock 
waves,  or  oblique  shocks  followed  by  either  a strong  oblique 
shock  or  a normal  shock.  After  each  shock  in  a system  of  shocks, 
the  stagnation  conditions,  sonic  conditions,  and  the  conditions 
which  would  exist  if  the  flow  were  isentropically  expanded  to  the 
pressure  existing  immediately  behind  the  first  shock  in  the 
system  are  determined.  Equilibrium  air  composition  is  assumed 
throughout  these  processes.  The  details  of  this  computer  pro- 
gram are  presented  in  Section  3.4. 

Some  results  of  calculations  with  this  computer  program 
for  30“  and  45“  half -angle  cones  appear  in  Table  2.  The 
numbering  of  states  in  Table  2 corresponds  to  the  numbering  in 
Figure  6.  For  30“  half-angle  cones  it  is  observed  from  the 
results  in  Table  2 that  the  flow  conditions  at  the  entrance  to 
the  interface  correspond  to  supersonic  flow  when  the  assximption 
is  made  that  ^4=^2  reflected  shock  is  a strong 

oblique  shock.  For  the  45“  cone,  the  Mach  number  M^,  is  just 
slightly  over  one  assuming  a Mach  reflection  at  the  rail  and  sub- 
sequent isentropic  expansion  to  P2* 
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It  should  be  emphasized  that  this  computer  program 
has  sufficient  flexibility  that  a number  of  hypotheses  on 
the  nature  of  the  rail-cone  shock  interaction  system  can  be 
tested. 


3.3.2  Interface  Sides 

The  sides  of  the  interface  region  are  open  to  the  gas 
flow  around  the  cylindrical  portion  of  the  projectile.  The 
pressure  in  this  region  of  the  flow  may  be  derived  from  the 
Prandtl-Meyer  expansion  of  the  gas  around  the  shoulder  of  the 
cone-cylinder.  The  pressure  along  the  sides  of  the  interface 
is,  therefore,  defined  by  the  pressure  increase  across  the 
attached  conical  shock  (or  conical  bow  shock  for  a blunted  cone) 
followed  by  the  Prandtl-Meyer  expansion  through  a turning  angle 
equal  to  the  cone  half-angle.  The  interface  side  pressures,  P^ 
are  most  informatively  compared  to  the  interface  entrance 
pressures  as  tabulated  in  the  last  column  of  Table  2. 

3.3.3  Interface  Exit 

The  interface  has  an  exit  boundary  at  the  trailing  edge 
of  the  projectile.  Determination  of  the  pressure  at  the  inter- 
face exit  is  perhaps  one  of  the  most  challenging  aspects  of  de- 
veloping a model  of  the  interface.  If  the  interface  flow  could 
be  considered  to  be  two-dimensional,  that  is,  if  there  were  no 
circumferential  pressure  relief,  the  method  used  by  Korkegi  and 
Briggs  could  be  used  to  estimate  the  exit  pressure.  They  were 
able  to  determine  the  local  pressure  in  the  interface  up  to 
the  point  where  the  boundary  layers  on  the  rail  and  the  projectile 
cylindrical  surface  merged  to  form  a turbulent  Couette  flow 
region.  After  the  mergence  of  the  boundary  layers,  they  assumed 
that  the  pressure  remained  constant  with  the  flow  at  the  exit, 
then  expanding  to  the  base  pressure  outside  of  the  interface 
flow  region. 
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TABLE  2.  FLOW  CONDITIONS  NEAR  RAIL  AND  AT  ENTRANCE* ** 


a 

cone 

angle 

"l 

"l 

atm 

deg. 

deg. 

“i 

m/sec 

P2/P,  or 
P4/P1 

“2 

m/sec 

P /P 

3'  1 

m/sec 

"4 

04/01 

T4 

“K 

TO 

4 

®K 

p;/Pi 

30“ 

10 

l.O 

33.57 

27.20 

3,470 

36.54 

2909 

510 

2410 

2.42 

3.91 

2796 

4498 

587 

30“ 

10 

.1 

33.57 

27.20 

3,470 

36.54 

2909 

516 

2385 

2.41 

3.90 

2971 

4229 

588 

30“ 

12 

1.0 

33.08 

27.34 

4.164 

51.73 

3507 

821 

2831 

2.50 

4.17 

3653 

5795 

931 

30“ 

12 

.1 

33.  os 

27.35 

4,164 

51.74 

3507 

832 

2778 

2.50 

4.29 

3507 

5303 

935 

30® 

14 

1.0 

32.64 

27.40 

4,858 

69.33 

4108 

1221 

3269 

2.57 

4.45 

4418 

7268 

1377 

30“ 

14 

.1 

32.64 

27.40 

4,858 

69.49 

4108 

1251 

3208 

2.58 

4.61 

4114 

6612 

1399 

30“ 

15 

1.0 

32.46 

27.47 

5.205 

79.17 

4409 

1462 

3498 

2.60 

4.57 

4801 

8015 

1645 

30“ 

15 

.1 

32.46 

27.58 

5,205 

79.46 

4408 

1512 

3437 

2.60 

4.75 

4445 

7236 

1687 

45® 

10 

1.0 

4946 

40.00 

3,470 

70.28 

2286 

124 

1248 

1.06 

5.21 

3935 

4334 

132.7 

45“ 

10 

.1 

49.46 

40.17 

3,470 

70.46 

2286 

125 

1221 

1.06 

5.44 

3695 

4020 

133.1 

45® 

12 

1.0 

4B.8B 

40.57 

4 ,164 

100.70 

2767 

180 

1446 

1.06 

5.70 

4919 

5458 

191.5 

45“ 

12 

.1 

48.88 

40.91 

4,164 

101.3 

2765 

182 

1405 

1.06 

6.0 

4504 

4976 

191.2 

45“ 

14 

1.0 

48.59 

41.07 

4,858 

137.2 

3241 

246 

1657 

1.05 

5.87 

6086 

6781 

260.9 

45“ 

14 

.1 

48.59 

41.50 

4,858 

138.2 

3238 

247 

1612 

1.05 

6.19 

5623 

6201 

261.6 

45“ 

15 

1.0 

48.30 

41.09 

5,205 

156.6 

3490 

283 

1769 

1.06 

5.93 

6707 

7421 

300 

45“ 

15 

.1 

48.30 

41.52 

5,205 

157.8 

3487 

285 

1715 

1.07 

6.28 

6169 

6722 

301 

*The  flow  regions  denoted  by  the  subscripts  correspond  to  the  numbering  of  the  flow  regimes  shown  in  Figure  7. 
l.e.  Pg  refers  to  the  stagnation  pressure  at  condition  4. 

**  6 = angle  between  bow  shock  and  trajectory. 

eee  d = angle  through  which  gas  turns  as  It  traverses  bow  shock. 
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TABLE  2. 


FLOW  CONDITIONS  NEAR  RAIL  AND  AT  ENTRANCE 


(Continued) . 


“l 

atm 

®1 

deg. 

•’3/^1 

m/sec 

"o'"! 

* 

gm/cm^ 

‘’’'3/P3 

k 

‘’"4 

gm/cin^ 

AAA 

's/p. 

30“ 

10 

1.0 

33.57 

blO.8 

604 

5.7b 

34.06 

.6007 

14.99 

4396 

.043 

9.34 

28.3 

1.58 

30“ 

10 

.1 

33.57 

516 

575 

6.03 

35.67 

.0607 

14.46 

4149 

.0043 

9.35 

27.70 

1.58 

30“ 

12 

1.0 

33.08 

821 

665 

6.26 

40.62 

.9651 

20.192 

5670 

.0608 

12.40 

27.11 

2.00 

30“ 

12 

.1 

33.08 

832 

625 

6.66 

43.16 

.0978 

19.27 

5203 

.0061 

12.15 

26.07 

2.00 

30“ 

14 

1.0 

32.64 

1221 

748 

6.50 

45.75 

1.436 

26.68 

7116 

.0815 

15.58 

25.75 

30“ 

14 

.1 

32.64 

1251 

705 

6.89 

49.11 

.1471 

25.46 

6494 

.0082 

15.04 

24.63 

2.43 

30“ 

15 

1.0 

32.46 

1462 

789 

6.59 

48.33 

1.719 

30.23 

7857 

.0931 

17.33 

25.09 

2.46 

30“ 

15 

.1 

32.46 

1512 

739 

7.04 

52.7 

,1778 

28.67 

7120 

.0094 

16.70 

23.97 

4S“ 

10 

1.0 

49.46 

124 

415 

8.36 

8.36 

.1458 

14.83 

4290 

.0827 

12.79 

27.90 

2.66 

45“ 

10 

.1 

49.46 

125 

393 

8.83 

6.83 

.0147 

14.1b 

3985 

.00829 

12.94 

27.19 

2.73 

45“ 

12 

1.0 

4s.ee 

180 

457 

9.10 

9.10 

.2119 

19.79 

5403 

.1184 

16.86 

26.36 

2.91 

45“ 

12 

.1 

48.88 

182 

430.5 

9.67 

9.67 

.0214 

18.81 

4930 

.01191 

16.88 

25.37 

2.98 

45® 

14 

1.0 

48.59 

246 

515 

9.43 

9.43 

.2893 

26.06 

6715 

.1614 

23.3 

24,87 

45“ 

14 

.1 

48.59 

247 

486 

10.00 

10.00 

.0291 

24.68 

6151 

.01625 

22.32 

24.00 

3.95 

45“ 

15 

1.0 

48.30 

283 

540 

9.63 

9.63 

.3328 

29.36 

7357 

.1842 

24.95 

24.19 

4.11 

45“ 

15 

,1 

48.30 

285 

505 

10.32 

10.32 

.0335 

27.71 

6676 

.0186 

25.1 

23.33 

*Siibscrlpt  H denotes  conditions  defined  at  the  temperatupe  of  the  rail  which  Is  assumed  to  be  at  a temperature 
of  300OK. 

**HH  Is  the  molecular  weight  of  the  air  at  state  3. 

***State  5 denotes  conditions  at  the  sides  of  the  Interface  (see  Section  3.3.2). 
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For  the  purpose  of  this  report,  we  have  used  the  Reynolds 
number  based  on  average  gap  height  Re^  to  serve  as  a guide  in 
estimating  the  pressure  at  the  exit.  If  the  interface  flow  is 

3 

laminar,  that  is,  if  Re^  is  less  than  10  , the  pressure  in  the 
interface  is  assumed  to  be  constant  at  the  entrance  value.  The 
procedure  is  justified  on  the  basis  of  information  given  on  pages 
81-83  of  Reference  8 where  it  is  shown  that,  under  conditions  of 
extremely  high  values  of  U^,  the  flow  in  a parallel  surface  slider 
bearing  with  an  imposed  pressure  gradient  in  the  direction  of 
flow  is  such  that  the  pressure  is  constant  at  the  initial  value 
until  at  the  very  end  of  the  bearing.  Beyond  that  point,  the 
pressure  drops  very  rapidly  to  the  surrounding  value. 

If  the  Reynolds  number  based  on  gap  height  lies  between 
10^  and  10^  the  Korkegi  and  Briggs  procedure  is  used.  It  is 
presumed  that  the  boundary  layers  on  the  rail  and  the  projectile 
will  merge  at  some  point  along  the  interface  and  that  the  pressure 
in  the  interface  will  decrease  continuously  until  the  boundary 
layers  merge,  after  which  the  pressure  will  remain  constant  until 
the  end  of  the  interface  is  reached.  Because  the  width  of  the 
interface  zone  in  the  present  problem  (essentially  the  width  of 
the  rail}  is  not  infinite  as  was  presumed  in  the  problem  addressed 
by  Korkegi  and  Briggs,  the  rate  at  which  the  pressure  decreases 
in  the  region  upstream  of  the  merged  boundary  layers  is  taken  to 
vary  inversely  with  the  ratio  of  rail  width  to  length  of  the 
interface  region. 

If  the  Reynolds  number  based  on  gap  height  is  larger  than 
10^,  it  is  presumed  that  the  boundary  layers  have  never  merged 
and  that  the  pressure  in  the  gap  decreases  from  the  entrance  pressure 
to  the  pressure  existing  along  the  sides  of  the  interface  at 
the  very  entrance  to  the  interface  region  so  that  the  pressure 
along  the  interface  is  presumed  to  be  constant  at  the  value 
fixed  by  the  pressure  existing  on  the  cylindrical  portion  of  the 
projectile  which  is  not  influenced  by  the  presence  of  the  rail. 
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3.3.4  Interface  Walls 

The  interface  is  bounded  by  two  solid  walls — the  rail  and 
the  projectile.  Each  of  these  walls  presents  a separate  boundary 
condition  for  the  problem. 

The  projectile  passes  rapidly  over  the  rail  so  that  any 
particular  section  of  the  rail  is  only  in  contact  with  the  in- 
terface for  a short  time.  This  time  is  too  short  for  signifi- 
cant heat  to  flow  into  any  portion  of  the  rail,  or  for  appreciable 
temperature  rise  in  the  rail,  so  the  rail  will  appear  to  be  a 
cool  constant  temperature  boundary. 

The  projectile  surface  remains  in  contact  with  the 
interface  during  the  entire  event.  There  can,  therefore,  be 
appreciable  heat  flow  to  the  projectile  and  the  temperature  will 
only  be  limited  by  the  heat  capacity,  thermal  conductivity,  and 
heat  capacity  associated  with  phase  changes  of  the  projectile 
material.  For  high  heat  fluxes,  phase  change  dominates  and 
the  projectile  surface  temperature  approaches  the  ablation  tem- 
perature of  the  material.  Strictly  speaking,  thermo-chemical 
equilibrium  must  be  established  and  processes  similar  to  those  in 
an  ablating  boundary  layer  occur. 

3.4  MATERIAL  PROPERTIES 

In  this  section  the  methods  and  procedures  used  to  define 
the  thermodynamic  and  transport  properties  of  the  various 
materials  involved  in  the  projectile-rail  guide  system  are  dis- 
cussed. The  material  properties  of  interest  are  the  gas  in  the 
system  (air) , the  solid  and  vapor  states  of  the  projectile 
material,  and  the  rail  itself. 

3.4.1  Air  Properties 

The  thermodynamic  properties  of  high  temperature  air  are 
required  in  a convenient  form  in  order  to  solve  the  normal  and 
oblique  shock  equations  for  the  high  Mach  numbers  and  cone 
angles  associated  with  the  projectiles.  For  a projectile  running 
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into  air  initially  at  BOO^K  with  a velocity  of  4.6  km/s,  the 
stagnation  temperature  reached  by  the  air  will  be  about  7500“K 
just  prior  to  entering  the  interface.  At  temperatures  below 
about  7500 ®K  the  primary  species  in  dry  equilibrium  are  N2» 

0,  N,  NOr  and  Ar  provided  that  the  density  of  the  air  is  no 
less  than  about  one  percent  the  density  of  eiir  at  standard  con- 
ditions and  no  greater  than  about  100  times  the  density  of  air 
at  standard  conditions.  In  addition,  the  individual  species 
comprising  air  still  behave  as  ideal  gases  under  these  limits  on 
temperature  and  density.  Thus,  the  dissociation  of  ©2  and  N2 
and  the  formation  of  nitric  oxide  (NO)  can  be  used  quite 
adequately  to  determine  the  composition  and  thermodynamic 
properties  of  equilibrium  air  within  the  above  limitations  on 
temperature  and  density. 

Existing  methods,  procedures,  and  parts  of  computer 
programs  available  at  UDRI  which  were  originally  developed  for 
combustion  of  hydrocarbon  fuels  with  air  were  used  in  conjunction 
with  the  elementary  gas  dynamic  equations  for  normal  and  oblique 
shocks  to  produce  a computer  program  which  allows  one  to  deter- 
mine the  thermodynamic  and  flow  conditions  downstream  of  a 
normal  shock,  an  oblique  shock,  or  a combination  of  oblique  shocks 
followed  by  a normal  shock.  The  program  is  very  efficient  both 
in  terms  of  running  time  and  central  memory  requirements.  The 
running  time  is  very  short  due  primarily  to  a subroutine  which 
determines  initial  guesses  on  the  composition  of  the  air  to 
within  20  percent  for  each  species. 

An  outline  of  a few  of  the  details  of  the  program 
follows.  The  key  part  of  the  program  is  a subroutine  which  solves 
for  the  equilibrium  composition  and  thermodynamic  properties  of 
air  at  a specified  temperature  and  pressure  (EQAIR) . The  gas 
dynamic  calculations  are  all  carried  out  in  the  Subroutine 
GASDYN.  The  oblique  shock  equations  are  solved  for  a specified 
upstream  Mach  number  temperature,  pressure,  and  wedge  half-angle. 
The  oblique  shock  equations  are  written  with  the  downstream 
temperature  (T) , pressure  (P) , and  shock  angle  (6^)  as  the 
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independent  variables.  (The  normal  shock  equations  are  just  the 
oblique  shock  equations  with  half-angle  equal  to  zero  and  shock 
angle  equal  to  90°. ) Perfect  gas  oblique  shock  and  thermo- 
dynamic relations  are  used  to  provide  initial  estimates  of  T,P, 
and  e^.  In  order  to  carry  out  the  Newton- Raphson  procedures, 
the  thermodynamic  derivatives  (9H/9T)p,  (SH/9P)^,  (9p/9T)p,  and 
(9p/9P)^,  where  p is  the  density  and  ri  is  the  enthalpy,  are 
required.  These  quantities  are  determined  in  a separate  sub- 
routine. After  T,p,  and  6^  have  been  determined,  the  equilibrium 
speed  of  sound,  a,  is  calculated.  The  calculation  of  the  speed 
of  sound  is  done  in  the  subrountine  SONIC  which  calculates  the 
derivative  (9n^/9T)g  for  each  species  where  s is  the  entropy. 

The  speed  of  sound  is  calculated  using  the  formula: 

a^  = {9P/9p)g  (24) 

Knowing  T,P,0^  and  a,  the  remaining  unknowns  such  as  Mach  number, 
density,  velocity,  entropy,  and  enthalpy  are  then  determined. 

The  program  also  contains  a subroutine  (CONSTS)  which 
performs  isentropic  compressions  and  expansions  for  equilibrium 
air.  This  subroutine  is  used  to  determine  stagnation  pressures 
and  temperatures  and  is  called  by  another  subroutine  PMEXP  which 
does  Prandtl-Meyer  expansions.  A subroutine  TRANSP  is  also 
included  which  uses  the  equilibrium  composition  to  determine 
the  viscosity  and  thermal  conductivity  of  the  equilibrium  air. 

This  program  facilitates  the  determination  of  interface 
entrance  conditions  for  any  desired  cone  angle  and  Mach  number 
and  provides  boundary  layer  edge  conditions  at  any  location  where 
heat  transfer  calculations  are  needed. 

3.4.2  Projectile  Material  Properties 

Mass  is  deposited  in  the  interface  by  the  decomposition 
of  the  slider  contact  surface.  The  interface  gases  are  very  hot 
and  the  resulting  heat  flow  to  the  slider  surface  is  undoubtedly 
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large  enough  to  cause  ablation  of  the  material.  The  mass  loss 

rate,  ih^^,  is  determined  by  the  heat  of  ablation/  L,  the  heat  flux 

to  the  wall,  d , and  the  ablating  (or  contact)  area,  A ; 
w 

m = q A /L.  (25) 

a ^w  c 

If  the  length  of  the  interface  is  B,  then  the  ablative 
mass  per  unit  width  of  gap  is; 

*a/”r  " 

It  is  then  necessary  to  evaluate  d„  and  L. 

3. 4. 2.1  Heat  of  Ablation,  L 

The  material  of  concern  in  this  study  is  Lexan^,  which 
is  a polycarbonate  material.  The  heat  of  ablation,  L,  is  defined 
as; 


L 


q /A  . 

^w  a 


(27) 


At  steady  state  ablation; 

L =^pdT  + AL  polymerization  + AL  vaporization,  (28) 


where; 


Cp  = 2.65  kJ/kg“K, 

AL  vap  = 0.68  MJ/kg, 

AL  poly  = 0.98  MJ/kg, 

or 

L = 2.65  X 10^  (T  - 295)  + 1.66  x 10®  J/kg.* 
This  equation  is  plotted  in  Figure  7 . 


®G.E.  Trademark  for  polycarbonate. 

*Private  communication  with  Glen  Ormbrek,  ATML. 


(29) 


(30) 
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Figure  7 . Heat  of  Ablation  versus  Ablation  Temperature  for 
Polycarbonate . 

The  ablation  temperature  of  polycarbonate  was  found  to 
be  between  450®C  and  550®C  from  a Thermal  Gravametric  Analysis 
(TGA)  at  atmospheric  pressure.  This  would  indicate  a heat  of 
ablation  in  the  range  of  2.1  to  2.3  MJ/kg.  Laser  ablation  studies 
indicate  a value  of  3.76  MJ/kg  (and  a temperature  of  about 
1100°K).*  Other  values  are  found  in  Reference  14  (5  MJ/kg)  and 
Reference  15  (4.66  MJ/kg). 

3.5  MODELS  FOR  DETERMINING  FORCES  AND  WEAR  FROM  A FLUID 

DYNAMIC  INTERFACE  BETWEEN  A PROJECTILE  AND  GUIDE-RAIL 

Before  presenting  a description  of  the  model  developed 
for  describing  projectile/rail  interface  regions,  a few  words 
are  in  order  concerning  some  practical  constraints  imposed  upon 
the  computational  approaches.  Ideally,  it  would  be  desired  to 
use  solutions  to  the  Navier-Stokes  equations  to  determine  the 
pressure,  temperature,  and  viscous  stress  distributions  which 


*Private  communication  with  Glen  Ormbrek,  AFML. 
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exist  in  the  interface  fluid  in  order  to  determine  the  forces, 
moments,  and  wear  (ablation)  on  the  projectile.  However,  the 
time  required  to  perform  these  computations  would  not  be  com- 
patible with  the  present  problem  even  if  computer  programs  were 
available  which  could  perform  the  necessary  computations  since 
these  computations  need  to  be  performed  thousands  of  times  during 
analysis  of  a single  passage  of  a projectile  through  the  rail 
system.  The  same  comment  applies  to  the  Reynolds  lubrication 
equation  which  is  very  difficult  to  solve  for  all  but  the  simplest 
cases.  Also,  the  lubrication  Reynolds  number  can  be  expected  to 
be  of  order  one  only  when  the  average  gap  height  is  extremely 
small  because  of  the  extremely  high  projectile  velocities.  Thus, 
the  lubrication  model  can  be  applied  for  only  those  cases  where 
a rail-projectile  collision  is  imminent  or  when  extremely  small 
clearances  exist  between  the  rail  and  projectile.  As  will  be 
shown  later  in  this  section,  a gap  height  small  enough  for  the 

laminar  flow  lubrication  model  to  apply  to  the  present  problem 

* • 

is  smaller  than  can  be  realized  physically  (in  the  range  of 
hundreds  of  nanometers) . 

3.5.1  Boundary  Layer  Considerations  for  the  Interface  Flow 

As  was  mentioned  in  Section  3.2.1,  boundary  layers 
develop  on  both  the  rail  and  the  projectile  surface  upstream 
of  the  interface  (see  Figure  3) . A boundary  layer  begins  to 
develop  on  the  rail  at  the  point  where  the  conical  shock  inter- 
sects the  rail  and  the  boundary  layer  on  the  cone  begins  to 
develop  at  the  tip  of  the  cone.  The  boundary  layer  on  the  conical 
portion  of  the  projectile  also  undergoes  a classical  boundary 
layer- shock  interaction  at  the  point  on  the  cone  where  the  re- 
flected shock  from  the  rail  intersects  the  cone.  At  the  juncture 
of  the  conical  portion  of  the  projectile  with  the  cylindrical 
portion,  the  projectile  boundary  layer  undergoes  another  change 
as  the  flow  on  the  projectile  surface  turns  to  enter  the  interface 
region.  The  character  of  the  rail  and  projectile  boundary  layers 
in  the  region  just  upstream  of  the  interface  can  have  a strong 
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influence  on  the  nature  of  the  flow  in  the  interface  itself. 

For  instance,  if  these  two  boundary  layers  have  not  merged  at 
the  entrance  to  the  interface,  then  an  inviscid  core  exists  be- 
tween the  two  boundary  layers,  and  the  initial  behavior  of  the 
flow  in  the  interface  is  dominated  by  gas  dynamic  considerations 
up  to  the  point  in  the  interface  where  the  separate  boundary 
layers  become  fully  merged.  After  mergence,  the  interface  flow 
is  dominated  by  viscous  effects.  If  the  gap  height  is  sufficiently 
large,  it  is  possible,  of  course,  that  the  two  boundary  layers 
never  merge.  If,  on  the  other  hand,  the  two  boundary  layers 
merge  upstream  of  the  interface,  then  the  flow  throughout  the 
interface  is  dominated  by  viscous  effects. 

Because  the  work  of  Korkegi  and  Briggs  comes  closest 
of  any  available  treatment  to  describing  the  boundary  layer  as- 
pects of  the  problem  at  hand,  it  bears  close  consideration. 

Korkegi  and  Briggs  studied  the  interface  region  of  a hypersonic 
slipper  bearing  on  rocket-boosted  sleds  guided  by  rails. 

Figures  8 through  11  describe  the  basic  features  of  the  problem 
they  addressed  and  the  idealizations  of  their  flow  model. 


Slipper 


Figure  8.  Typical  Rocket  Sled  Configuration  and  Slipper  Details. 
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a)  Probable 


b)  Assumed 


Figure  10.  Bow  Shock  Wave. 
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Figure  11.  Boundary  Layer  Growth  in  Gap. 

In  Reference  16  many  simplifications  were  made  on  the 
sled-slipper  motion  problem  such  as  assuming:  a steady  state 
configuration;  locally  two-dimensional  flow;  and  that  the  flow 
model  consists  of  a laminar  stagnation  region  at  the  leading 
edge  of  the  slipper  followed  by  separate  turbulent  boundary  layers 
on  slipper  and  rail  (shock  tube  type)  eventually  merging  into 
a full  shear  layer  which  tends  toward  a Couette  flow  asymptote. 

The  inlet  shock  structure  was  determined  to  be  dependent  upon: 
the  shape  of  the  leading  edge  of  the  slipper;  the  ratio  of  gap 
height  to  slipper  thickness;  and  the  mass  flow  through  the  gap 
which,  in  turn,  dictates  the  location  of  the  dividing  streamline 
separating  the  flow  entering  the  gap  from  that  going  over  the 
slipper.  The  mass  flow  was  found  to  be  dependent  on  the  pressure 
field  generated  by  the  shock  compression.  Hence,  the  mass  flow 
and  shock  structure  were  considered  coupled.  In  the  sled-slipper- 
track  problem,  the  gap  heights  were  much  smaller  than  the  blunt 
slipper  leading  edge;  thus,  it  was  assumed  that  the  bow  wave  was 
normal  between  the  dividing  streamline  and  the  rail,  so  that  all 
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the  flow  entering  the  gap  was  compressed  through  a single  normal 
shock  as  shown  in  Figures  8 through  11. 

One  of  the  major  differences  between  the  slipper/rail 
problem  and  the  projectile/rail  problem  is  that  slender  nose 
cones  will  create  a Mach  reflection  of  the  bow  wave  on  the  rail 
rather  than  the  single  normal  shock.  If  the  semi-cone  angle 
becomes  large,  i.e.,  >45“  the  assumption  of  a single  normal 
shock  at  the  gap  entrance  would  be  reasonable;  however,  three- 
dimensional  pressure  relief  is  present  along  the  projectile-rail 
interface.  The  single  most  important  difference  between  the 
slipper/rail  problem  and  the  projectile/rail  problem  is  the  fact 
that  in  the  slipper/rail  problem,  two-dimensional  flow  could 
reasonably  be  presumed  to  exist  over  the  interface,  whereas  in 
the  projectile/rail  problem  the  two-dimensional  assumption  is 
never  strictly  valid,  and  grossly  misrepresents  the  flow  situation 
when  the  boundary  layers  in  the  flow  have  not  merged  to  form  a 
viscous  dominated  flow  region. 

While  the  entire  interface  model  of  Korkegi  and  Briggs 
is  not  applicable  to  the  projectile/rail  problem,  some  of  the 
thinking  that  went  into  their  model  is  applicable.  In  addition, 
their  work  on  fully  developed  turbulent  Couette  flow  is  directly 
applicable.  In  principle,  their  basic  approach  to  determine 
both  the  point  in  the  interface  where  the  boundary  layers  merged 
and  the  pressure  distribution  upstream  of  the  mergence  point  can 
be  used  to  advantage  in  the  rail/projectile  problem.  They  cal- 
culated the  displacement  thicknesses  of  the  two  boundary  layers 
which  defined  the  boundaries  of  a two-dimensional  flow  channel 
and  then  used  isentropic  flow  gas  dynamic  considerations  to  de- 
termine the  pressure  distribution  along  the  flow  channel  or 
interface.  Because  they  assumed  two-dimensional  flow,  their 
model  was  relatively  straightforward  to  implement.  In  the  present 
problem  their  procedure  would  have  to  be  considerably  altered 
to  account  for  the  pressure  relief  afforded  by  the  circumferential 
flow  from  the  center  of  the  interface  to  the  sides  of  the  interface. 
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No  detailed  computations  of  the  displacement  thicknesses 
were  attempted  in  the  present  analysis  because  of  the  difficulty 
involved  in  performing  such  calculations  for  three-dimensional 
boundary  layers.  No  simple  formulations  exist  for  determining 
the  displacement  thickness  of  three-dimensional  boundary  layers . 
Therefore,  no  effort  was  made  to  develop  a rigorous  method  for 
determining  the  exact  point  of  mergence  of  the  rail  and  projectile 
boundary  layers.  Instead,  a simple  application  of  the  results 
of  Korkegi  and  Briggs  was  made  tempered  with  the  fact  that 
pressure  distribution  upstream  of  the  mergence  point  is  adjusted 
such  that  the  pressure  in  the  interface  rapidly  achieves  ec[uilibration 
with  the  pressure  existing  along  the  sides  of  the  interface  when- 
ever the  boundary  layers  have  not  already  merged  upstream  of  the 
interface.  The  results  of  computations  of  the  rail  boundary 
layer  are  given  in  the  next  paragraph  in  order  to  provide  some 
insight  into  the  conditions  under  which  mergence  of  the  boundary 
layers  upstream  of  the  interface  region  could  be  expected. 

3.5.2  Rail  Boundary  Layer  Computations 

From  the  data  given  in  Table  2 some  basic  boundary  layer 
properties  were  calculated  for  the  boundary  layer  which  develops 
along  the  rail  having  its  leading  edge  at  the  juncture  of  the 
strong  shock  and  the  rail.  The  purpose  of  these  boundary  layer 
calculations  was  to  determine  whether  this  boundary 
layer  is  a laminar  or  turbulent  boundary  layer  by  the  time  it 
reached  the  entrance  of  the  gap,  and  also  to  determine  the  thick- 
ness of  the  boundary  layer  at  the  entrance  to  the  gap.  Knowledge 
of  the  boundary  layer  thickness  and  Reynolds  number  at  the  en- 
trance to  the  gap  allows  one  to  determine  the  conditions  under 
which  the  flow  at  the  gap  entrance  could  be  expected  to  have  a 
significant  unifoirm  velocity  region  (an  inviscid  core)  and  whether 
or  not  the  shear  flow  region  could  be  expected  to  already  be 
turbulent.  If  the  boundary  layer  thickness  is  greater  than  the  gap 
height,  no  significant  uniform  mean  flow  region  exists  at  the 
gap  entrance,  and  the  flow  in  the  entire  gap  length  can  be 
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considered  viscous  dominated  with  little  pressure  relief  along 
the  interface.  If,  at  the  same  time,  the  shear  flow  region 
at  the  gap  entrance  is  laminar,  then  the  initial  portion  of  the 
interface  flow  region  can  be  expected  to  be  laminar  regardless 
of  the  laminar  or  turbulent  state  achieved  well  into  the  gap. 
Conversely,  if  the  boundary  layer  is  already  turbulent  at  the 
entrance  to  the  gap,  then  the  flow  in  the  interface  region  could 
have  one  of  two  different  basic  characteristics  depending  on  the 
Reynolds  number  of  the  interface  flow,  If  is  such  that 

the  basically  Couette  flow  is  turbulent,  then  the  entire  in- 
terface flow  regime  can  effectively  be  considered  to  be  a tur- 
bulent Couette  flow.  If  Re.  is  such  that  the  Couette  flow  is 

n 

laminar,  then  the  initial  turbulence  will  be  damped  out  and  the 
entire  interface  flow  will  effectively  be  laminar. 

Listed  in  Table  3 are  the  boundary  layer  characteristics 
expected  at  the  entrance  to  the  gap.  Before  discussing  these 
results,  the  nomenclatare  of  the  table  quantities  must  be  de- 
fined. The  column  listed  as  Ax  represents  the  distance  along  the 
rail  from  the  strong  shock  to  the  gap  entrance.  The  tabulated 
values  of  Ax  are  calculated  based  on  a zero  gap  height,  h,  and 
for  a projectile  diameter  of  41  mm.  If  h becomes  a significant 
fraction  of  the  cone  radius,  then  Ax  will  become  smaller  than 
the  tabulated  values.  As  indicated  in  Table  3,  the  values  of 
Ax  for  a 45®  half-angle  cone  are  estimated  values  based  on  the 
laser  shadowgraph  presented  in  Figure  30  of  Reference  1.  The 
precise  location  of  the  juncture  of  the  conical  shock  and  the 
strong  shock  cannot  be  easily  calculated.  Therefore,  recourse 
was  made  to  estimation  from  the  laser  shadowgraph.  The  kine- 
matic viscosities  v at  the  rail  tabulated  in  Table  3 were 

w 

calculated  at  a temperature  of  300 “K  using  the  pressure  at  either 

point  3 or  4 (defined  in  Figure  6)  as  the  case  may  be.  The 

-4 

dynamic  viscosity  of  air  at  300®K  was  taken  to  be  1.8  x 10 
gm/cm-sec  and  the  perfect  gas  law  was  used  to  calculate  the 
density  of  the  air  at  the  rail.  The  column  listed  as  y^  gives 
the  laminar  boundary  layer  thickness  based  upon  the  point  in 
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TABLE  3.  BOUNDARY  LAYER  CHARACTERISTICS 
GUIDE-RAILS . 


a 

cone  angle 

"l 

•■i 

atm 

Ax 

cm 

VwsXlO 

cm^/sec 

,/v** 

V ,.xlO 

W^l 

cm^/sec 

”l 

"4 

y^xlO^ 

cm 

■■ 

1.0 

.457 

3.0 

41.9 

1.44 

2.36 

.1 

.457 

30.0 

419 

1.45 

7.50 

30® 

12 

1.0 

.401 

1.87 

?9.6 

1.47 

2.27 

30® 

12 

.1 

.401 

18.4 

295 

1.50 

7.09 

30® 

14 

1.0 

.348 

1.25 

22.1 

1.49 

2.14 

30® 

14 

.1 

.348 

12.2 

220 

1.51 

6.57 

30® 

15 

1.0 

.325 

1.05 

19.3 

1.49 

2.07 

30® 

15 

.1 

.325 

10.1 

■ 192 

1.51 

6.37 

45® 

10 

1.0 

♦ 

12.3 

21.8 

2.78 

4.00 

45® 

10 

.1 

122 

217 

2.84 

12.91 

45® 

12 

1.0 

.7 

8.5 

15.2 

2.88 

4.09 

45® 

12 

.1 

(estimated 

from 

84.1 

151 

2.96 

13.10 

45® 

14 

1.0 

Fig.  30 

6.22 

11.1 

2.93 

4.51 

45® 

14 

.1 

of  Ref.l) 

1 

61.9 

ill 

3.01 

13.86 

45® 

15 

1.0 

1 

5.4 

9.77 

2.94 

4.39 

45® 

IS 

.1 

53.7 

96.8 

3.03 

14.11 

FOR  CONE  PROJECTILES  TRAVERSING 


gm/cin- 

&ec 


P4*10 

3 

gm/cm 


**2. 

cm  /sec 


xlO 
Ax 


-5 


Re^  XlO 


MW. 


7.43 

4.602 

.162 

1.31 

2.1 

28.91 

7.68 

.4594 

1.67 

.135 

.21 

28.80 

9.11 

4.905 

.106 

1.35 

3.6 

20.43 

8.39 

.5007 

1.68 

.165 

.36 

27.85 

9.67 

5.234 

.185 

1.45 

5.58 

27.37 

9.12 

.5435 

1.68 

.176 

.56 

26.41 

10.4 

5.374 

.194 

1.40 

6.85 

26.75 

9.70 

.5596 

1.73 

.171 

.69 

25.70 

8.92 

6.130 

.146 

18.97 

4.04 

28.17 

8.63 

.6402 

1.35 

2.15 

.404 

27.56 

10.2 

6.661 

.154 

23.2 

6.96 

26.71 

9.77 

.7057 

1.38 

2.75 

.670 

25.75 

]1.76 

6.920 

.170 

25.5 

11.1 

25.19 

11.29 

.7281 

1.55 

2.95 

1.11 

24.31 

12.48 

6.980 

.179 

26.10 

13.5 
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the  boundary  layer  where  the  velocity  is  within  one  percent 
of  the  free- stream  value.  Values  for  y^  were  calculated  from 
results  presented  in  References  17,  18,  and  19.  The  boundary 
layer  thickness  was  estimated  using  the  following  relation; 


^6 


(31) 


G was  taken  to  be  2.0  based  primarily  on  the  results  presented  in 

Reference  17  for  values  of  Uj^/U^  ranging  from  2 to  4.  The  values 

of  T^^/T^  for  30“  half-angle  cones  do  not  correspond  to  values 

of  T,,./T.  obtained  for  a single  normal  shock  so  that  G can  be 
^4  4 

off  by  perhaps  as  much  as  20  percent.  The  Reynolds  number  R^^^ 
listed  in  Table  3 is  defined  in  the  following  way; 


(Ui  - Ax 

*^Ax  U— 


(32) 


This  is  the  appropriate  Reynolds  number  to  use  for  determining 

the  transition  Reynolds  number  for  boundary  layers  formed  by  a 

moving  shpck  wave  (References  19,  20,  21).  The  dynamic  viscosity 

V.  was  determined  from  Reference  22. 

4 

The  last  column  of  Table  3 gives  the  Reynolds  number  of 
the  interface  flow,  Rej^,  based  on  the  following  definition; 


(33) 


where  is  equal  in  magnitude  to  U^.  A value  of  h of  25.4  ym 
was  used  to  define 

The  results  presented  in  Table  3 provide  the  following 
insights;  (1)  the  thickness  of  the  boundary  layer  for  bar 

is  on  the  order  of  25  ym,  (2)  if  the  gap  height  is  significantly 
greater  than  25  ym,  a uniform  mean  flow  region  will  exist  and 
the  problem  becomes  similar  to  that  treated  by  Korkegi  and  Briggs, 
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(3)  the  boundary  layer  Reynolds  number,  Re^^,  for  30®  half -angle 
cones  is  well  below  the  transition  Reynolds  number  range  of 
1 to  10  X 10®  for  this  type  of  boundary  layer”’  “ so  that  the 
boundary  layer  on  the  rail  at  the  entrance  to  the  interface  flow 
region  is  always  laminar,  and  (4)  the  interface  flow  Reynolds 
number  based  on  gap  height  is  at  least  an  order  of  magnitude 
greater  than  the  commonly  accepted  value  of  transition  Reynolds 
number  of  1 to  4 x 10^.  Thus,  only  for  gap  heights  smaller  than 
2 . 5 urn  could  completely  laminar  Couette  flow  be  expected  in  the 
interface. 

Thus,  the  following  conclusions  can  be  drawn  based  on 
the  proceeding  analysis: 

(1)  Laminar  Couette  flow  should  not  be  expected  to 
occur  unless  gap  heights  of  less  than  2.5  urn 
are  to  be  considered. 

(2)  For  gap  heights  in  the  range  of  2.5  um  to  about 
25  ym  turbulent  Couette  flow  can  be  expected  over 
the  entire  interface  flow  regime. 

(3)  For  gap  heights  significantly  greater  than  about 
25  ym,  the  flow  in  the  initial  portion  of  the 
interface  region  will  contain  a basically  uniform 
core  which  implies  that  gas  dynamic  considerations 
need  to  be  included  in  the  interface  flow  analysis. 

The  interface  flow  situation  then  becomes  one  of  two  separate 
boundary  layers  developing  one  on  the  rail,  and  one  on  the 
projectile  with  a bulk  flow  between  the  two  boundary  layers. 

The  behavior  of  the  bulk  flow  becomes  quite  sensitive  to  the 
shapes  of  the  developing  boundary  layers.  For  large  gap  heights, 
the  bulk  flow  will  also  be  influenced  by  the  circumferential 
pressure  gradient  between  the  center  of  the  rail  and  the  edge 
of  the  rail.  Because  a circumferential  pressure  gradient  was 
not  included  in  Korkegi  and  Briggs  analysis  of  a hypersonic 
slipper  bearing,  the  direct  application  of  their  approach  to  the 
present  problem  is  not  possible. 
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3.5.3  Simplified  Fluid  Dynamic  Models 

Based  on  the  results  of  the  previous  section,  it  is 
clear  that  three  basic  flow  models  are  needed  to  describe  the 
interface  phenomena  in  order  to  provide  complete  coverage  of 
the  entire  Reynolds  number  range  from  0 to  ».  For  the  Reynolds 
number  range  of  0<Rej^<lxl0  , the  laminar  Couette  flow  model  is 
used.  For  the  Reynolds  number  range  of  10^  to  10^  the  turbulent 
Couette  flow  mode  of  Korkegi  and  Briggs^*’*  is  used.  For 
Reynolds  numbers  greater  than  1x10^  the  projectile  boundary 
layer  is  presumed  to  be  a fully  turbulent  boundary  layer  acting 
as  if  no  rail  were  present.  The  methods  used  to  implement  the 
above  three  models,  including  the  effects  of:  non-parallel 
boundaries;  ablation  or  wear;  pressure  buildup;  and  the  forces 
and  moments  on  the  projectile,  are  discussed  in  Section  IV. 

3.6  PROJECTILE  MATERIAL  DEFORMATION 

During  the  initial  time  period  of  this  investigation  it 
was  assumed  that  elastic  deformation  of  the  projectile  material 
would  not  be  important.  The  gap  heights  were  expected  to  be  on 
the  order  of  25.4  ym  (0.001  in.).  Thus,  the  maximum  possible 
angular  rotation  of  the  projectile  and  consequently  the  normal 
forces  and  aerodynamic  moments  associated  with  rotation  of 
the  projectile  were  expected  to  be  insignificant  compared  to  the 
normal  forces  and  moments  due  to  interfacial  forces.  However, 
as  test  results  became  available,  it  was  apparent  that  initial 
gap  heights  were  much  greater  than  25  ym  (see  footnote  on  page  3) . 
Analyses  of  the  forces  required  to  keep  the  projectile  within  the 
confines  of  the  rails  with  the  larger  than  anticipated  rotation 
of  the  projectile  showed  that  pressures  in  the  interfacial  region 
were  such  that  elastic  deformation  of  the  projectile  would  be  an 
important  phenomenon.  Thus,  the  fluid  dynamic  processes  in  the 
interfacial  region  are  not  just  simply  hydrodynamic,  but  rather 
elastohydrodynamic  in  nature. 
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SECTION  IV 

COMPUTER  SIMULATION  OF 
PROJECTILE-RAIL  INTERACTION  PHENOMENA 

In  this  section  we  describe  the  methods,  procedures,  and 
equations  that  are  used  in  the  Projectile  Dynamics  Computer 
Program  to  compute  the  forces,  moments,  and  ablation  of  the  pro- 
jectile due  to  the  presence  of  a rail.  The  computations  are 
carried  out  individually  for  the  top  rail,  the  bottom  rail,  and 
the  two  side  rails  and  then  the  net  quantities  are  determined. 
This  section  also  describes  the  methods  used  to  determine  those 
aerodynamic  forces  and  moment  coefficients  which  are  calculated 
in  the  program.  This  section  consists  of  three  subsections. 
Section  4.1  describes  how  the  zero  angle-of -attack  total  drag 
coefficient,  Cq^  is  calculated.  This  quantity  includes  the  con- 
tribution of  the  shock  interaction  drag  to  the  total  drag,  but 
does  not  include  drag  due  to  interface  forces.  Section  4.2 
describes  subroutine  INTER,  in  which  the  interface  effects  on 
the  projectile  are  determined.  These  effects  are  the  viscous 
shear  forces,  the  normal  forces  (pressure  forces),  and  the 
ablation  rates.  The  moments  on  the  projectile  due  to  both  inter- 
face forces  and  shock  interaction  forces  are  also  computed  by 
INTER. 

4.1  COMPUTATION  OF  Cn 

0 

In  this  section  we  describe  the  methods  used  to  determine 
the  drag  coefficient  which  includes  drag  due  to  aerodynamic 
and  shock  interaction  effects.  In  essence,  this  section  is  a 

description  of  subroutine  DRAGCO.  The  factor  Cn  is  defined  as; 

0 


+ + C, 


+ C. 


'SI 


5e 


(34) 
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The  first  three  components  of  Cn  are  standard  aerodynamic  drag 

0 

coefficients;  namely,  pressure  drag  on  the  forebody  of  the  pro- 
jectile, viscous  drag  on  the  entire  projectile,  and  base  drag. 

The  reference  area  used  to  define  Cn  is  the  base  area  of  the 

0 

projectile. 


4.1.1  Pressure  Drag,  Cn 

P 

Two  different  equations  for  Cj^^  are  used  in  subroutine 
DRAGCO.  Both  are  valid  for  blunted  or  sharp  cones.  The  first 
equation  for  Cn  is  taken  from  Reference  25  and  is  valid  for  Mach 
numbers  in  the  range  of  5 to  25,  cone  half-angles  up  to  15®  and 
blunted  ratio,  C,  defined  as  the  nose-tip  radius,  r^,  divided  by 
the  base  radius  of  the  cone,  rj^,  from  0 to  0.3; 

c c c 

\ = Cl  Sin(5)  ^ + 3-57(5)  ^ ' (35) 


where ; 

= 1.944  + 1.872/M  - 17.0/M^  + 38.19/m3, 

C2  = 1.931  + 0.8635/M  - 8.063/M^  + 12.21/m3, 

C3  = 11.433  + 34.96/M  - 921. 5/M^  + 2607. 3/M^, 

C,  = 0.5359  + 0.09964/M  + 10.769/m3  - 104.21/m3 

^ 4 

+ 209.43/M^,  and 

C,  = 3.296  + 2.997/M  - 14.378/K^  + 154.67/m3. 


If  the  above  equation  is  not  applicable  then  the  following,  less 
accurate  equation  for  Cjj  , based  on  Newtonian  theory  of  hyper- 
sonic drag  is  used 

Cjj  = — ^(1  - ^ cos^o)  (2  sin^a)  + cos^a^  , (36) 

where  C for  Mach  numbers  greater  than  6 has  the  value  of 
Pmax 

1.83  for  air. 
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4.1.2  Viscous  Drag,  Cn 

The  viscous  drag  is  the  total  of  the  contributions  of 
skin  friction  drag  on  the  projectile  forebody  and  the  cylindrical 
portion  of  the  projectile.  Both  of  these  contributions  are  com- 
puted from  the  same  basic  skin  friction  coefficient  formulation 
based  on  a turbulent  flow  skin  friction  formula  for  incompressi- 
ible  flow  and  corrected  to  compressible  flow  conditions  by  the 
reference  enthalpy  method. “ 

The  drag  coefficient  for  the  cylindrical  portion  of  the 
projectile  then  becomes 


cot  0 , (37) 

where  Rejj,  is  the  Reynolds  number  based  on  free  stream  conditions 
and  the  slant  length  of  the  cone.  H*  is  the  Eckert  reference 
enthalpy  given  by;** 

^ = 0.5  + 0.5  (T^g)  (T^/Tg)  + 0.374  (38) 

The  subscript  e in  the  above  equations  denotes  conditions  evalu- 
ated just  outside  the  boundary  layer  on  the  cone.  This  would 
correspond  to  conditions  evaluated  at  state  2 in  Table  2.  Con- 
ditions at  state  2 are  evaluated  in  subroutine  GASDYN.  The 
subscript  « in  the  above  equations  refers  to  free  streeun  condi- 
tions and  this  denotes  conditions  at  state  1 in  Table  2.  B in 

c 

the  above  equations  is  a blowing  correction  to  account  for  mass 
addition  to  the  boundary  layer  due  to  ablation.  B^  is  evaluated 
from  a curve  fit  to  data  given  in  Figure  18  of  Reference  13.  The 
skin  friction  drag  coefficient  for  the  cylindrical  portions  of 
the  projectile  is  given  by; 
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(39) 


where  B denotes  the  length  of  the  projectile  afterbody,  Rob  de- 
notes the  Reynolds  number  based  on  the  length  B,  and  subscript 
5 denotes  conditions  existing  just  outside  the  boundary  layer  on 
the  projectile  afterbody.  Conditions  at  state  5 are  evaluated 
in  the  Prandtl-Meyer  expansion  subroutine  PMEXP.  The  last  term 
in  Equation  (39)  accounts  for  the  use  of  the  base  area  of  the  pro- 
jectile as  the  reference  area  in  defining  the  drag  coefficient. 


4.1.3  Base  Pressure  Drag,  Cp 

B 

The  following  expression  for  base  pressure  drag  was  taken 
from  Reference  13 ; 


B nj  bee 

where  Rj^/Pg  is  evaluated  from  a curve  fit  to  data  given  in 
Figure  19  of  Reference  13.  Pj^  is  the  base  pressure  and  P^  is 
defined  as  in  Section  4.1.2. 


4.1.4 


Shock  Interaction  Drag. Cn 
^SI 

The  basic  expression  used  to  define  C^gj  is  given  by; 


CD 


SI 


4 W^Axg^sin  P(Psi-Pg) 
l/2p„  uj  A. 


(41) 


AXgj  is  the  distance  from  the  base  of  the  conical  forebody  to  the 
point  on  the  cone  where  the  rail  shock  intersects  the  cone  sur- 
face, measured  along  the  cone  surface.  P__  represents  the 
average  pressure  existing  on  the  cone  surface  over  the  rail 
shock  interaction  region  of  length  Ax_t.  sin  a in 

Equation  (41)  represents  the  projected  area  normal  to  the  x axis 
over  which  the  average  pressure  increase  due  to  the  rail  shock 
interaction  (P-t-P.^^  acts.  Note  that  P„t“P.,.  used  in  Equation 
(41)  rather  than  Pg^  since  the  drag  force  in  the  x-direction  due 
the  pressure  P^  has  already  been  accounted  for  in  the  expression 
for  Cj3p.  The  factor  4 in  Equation  (4  0)  accounts  for  the 
existence  of  four  rails.  A correction  factor  having  a value 
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greater  than  unity  should  probably  be  applied  to  Equation  (40) 
whenever  the  rail-projectile  shock  interaction  is  a weak  shock 
interaction  to  account  for  the  fact  that  the  pressure  immediately 
downstream  of  the  point  of  intersection  of  the  shock  with  the  pro- 
jectile is  higher  than  the  pressure  in  the  region  adjacent  to  the 
shock  interaction  region.  For  a weak  shock  interaction  the 
velocity  of  the  gas  in  the  high  pressure  zone  is  still  supersonic. 
Thus,  a complicated  three-dimensional  shock  system  will  be  pre- 
sent which  will  cover  an  area  somewhat  greater  than  that  given  by 
W^AXgj.  A first  approximation  to  this  correction  factor  could  be 
obtained  from  a loose  analog  to  the  two-dimensional  problem  of  two 
parallel  supersonic  streams  at  different  pressures  suddenly  con- 
fronting each  other.  An  example  of  this  type  of  problem  is  given 
on  page  480  of  Reference  26  from  which  it  is  clear  that  the  high 
pressure  region  will  expand  outward  toward  the  low  pressure  region 
and  this  expansion  occurs  such  that  the  average  pressure  times 
the  affected  area  is  larger  then  the  pressure  of  the  high  pressure 
stream  times  just  the  area  covered  by  the  high  pressure  stream. 

The  determination  of  this  correction  factor  was  not  attempted  in 
the  course  of  this  effort. 

The  geometric  conditions  for  an  attached  shock  on  a sharp 
cone  are  shown  in  Figure  12.  The  reflection  angle  of  the  shock 
reflecting  off  the  rail,  0^  in  Figure  12,  is  determined  in  sub- 
routine GASDYN  for  a sharp  cone,  -^^gj  is  calculated  from  the 
following  formula  obtained  from  a geometric  analysis  of 
Figure  12; 

r.  (r.  + h)  sin(0  + 6 ) 

^SI  sin  a sin  0 sin (a  + 0 ) ’ ^ ' 

s r 

The  shock  interaction  pressure  Pg^  is  a quantity  which 
has  not  been  defined  rigorously  in  the  present  effort:  In  sec- 
tion 3.3.1  we  discuss  the  basic  difficulty  in  trying  to  define 
the  exact  nature  of  the  rail  shock  interaction.  Certainly  the 
pressure  existing  over  the  length  ^Xg^  is  not  constant  since  the 
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Figure  12.  Sketch  Showing  Definition  of  Shock  Interaction 
Distance 

flow  will  expand  outward  from  the  high  pressure  region.  Then 
there  is  the  question  concerning  the  number  of  weak  reflected 
shocks  which  can  exist  between  the  rail  and  the  projectile  for 
those  situations  where  the  reflected  shock  from  the  rail  is  a 
weak  shock.  The  present  best  estimate  of  the  authors  is  that, 
for  slender  cones,  the  rail  shock  impinging  on  the  cone  will 
experience  a more  or  less  regular  reflection  at  the  cone  sur- 
face, but  the  shock  reflecting  off  the  cone  will  probably  not 
penetrate  back  to  the  rail  since  it  will  be  diminished  in 
strength  as  it  moves  into  gas  which  itself  is  expanding  out- 
ward in  a direction  parallel  to  the  rail  surface  but  perpendi- 
cular to  the  direction  of  motion  of  the  projectile.  On  this 
basis,  perhaps  the  best  approach  to  follow  in  defining  for 
slender  cones  would  be  to  use  an  average  pressure  obtained  for 
the  case  of  no  reflection  from  the  cone  surface  and  the  pressure 
obtained  for  the  case  of  a regular  reflection  at  the  cone  surface. 
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The  program  is  set  up  such  that  the  user  can  specify  the 
total  number  of  reflected  waves  which  he  wants  to  consider  (up 
to  a total  of  4).  In  determining  the  value  of  Pgjr  the  prograim 
will  use  the  pressure  existing  behind  the  first  reflected  shock 
(conditions  behind  the  first  reflection  from  the  rail,  e.g., 
in  Figure  6}  if  the  user  calls  for  one  reflected  shock,  and  will 
use  the  pressure  existing  behind  the  second  reflected  shock  if 
two  reflected  shocks  are  specified.  If  more  than  two  reflected 
shocks  are  called  for,  the  program  will  compute  the  thermodynamic 
and  flow  conditions  behind  each  of  the  reflected  shocks,  but  will 
use  the  pressure  behind  the  second  shock  as  Pgj* 

In  concluding  this  paragraph,  it  should  be  mentioned 
that  Appendix  B contains  a description  of  the  input  data 
required  to  operate  the  progrcun  and  a description  of  the  indi- 
vidual subroutines  which  comprise  the  total  program. 

4.2  FORCES,  MOMENTS,  AND  ABLATION  DUE  TO  THE  INTERFACE 

In  this  section  we  describe  subroutine  INTER.  The 
sequence  of  computations  in  this  subroutine  is  presented  below. 

(1)  Gap  heights  at  the  front  and  rear  of  the  inter- 
face on  each  rail  (top,  bottom,  and  side)  are 
read  from  COMMON.  An  average  gap  height  and  a 
minimum  gap  height  for  each  of  these  three 
interfaces  are  then  computed  and  stored. 

(2)  Reynolds  numbers  based  on  these  average  and 
minimum  gap  heights  are  then  computed  and 
stored. 

(3)  Viscous  shear  forces  (tangential  forces  on  the 
projectile)  are  then  computed  followed  by  com- 
putation of  ablation  rates.  One  of  the  three 
fluid  dynamic  models  discussed  in  Section  3.5 
is  used  to  determine  an  average  skin  friction 
coefficient  for  the  interface  depending  on  the 
value  of  the  Reynolds  number  based  on  gap 
height  or  (if  the  projectile  axis  is  not  par- 
allel to  the  rail)  the  Reynolds  number  based  on 
the  minimum  gap  height.  Through  the  use  of 
Reynolds  analogy  for  heat  and  momentum  trans- 
fer, an  appropriate  heat  transfer  coefficient  is 
then  computed  followed  by  computation  of  the 
ablation  rate. 
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(4)  Normal  forces  due  to  the  pressure  distribution 
in  the  interfaces  are  computed  along  with  the 
moments  caused  by  these  normal  forces. 

(5)  Normal  forces  due  to  the  rail  shock  interaction 
are  computed  along  with  the  moments  caused  by 
both  the  normal  and  tangential  components  of  the 
rail  shock  interaction  forces. 

(6)  The  net  forces  in  the  tangential  and  normal  direc- 
tions are  then  computed  along  with  the  net  moment 
about  the  pitching  axis. 

The  reader  will  observe  that  some  effects  of  the  rail  shock  inter 
action  are  included  in  this  subroutine.  The  reason  for  this  is 
that  the  net  normal  force  and  the  net  moment  about  the  pitching 
axis  both  are  sensitive  to  the  position  of  the  projectile  re- 
lative to  the  rails.  From  Equation  (41)  it  is  observed  that 
Ax__,  and  thus  the  shock  interaction  force,  is  a function  of  gap 
height.  For  slender  cones  in  particular,  ^Xgj  is  sensitive 
enough  to  small  changes  in  the  gap  height,  h,  that  the  effect  of 
variations  in  h must  be  included  in  determinimg  the  net  normal 
force  and  the  net  moment.  The  net  tangential  force,  however,  is 
not  sensitive  to  h and  this  effect  is  already  included  in  the 
total  zero-angle-of -attack  drag  coefficient. 

Track  G projectiles  wear  data  became  available  near 
the  end  of  this  effort.  Relatively  large  rotations  of  the  pro- 
jectile were  found  to  occur  (5“  to  7.5“  at  6000  m/s  and  1 atm.) 
The  rotation  of  the  projectile  as  well  as  the  gap  height  affect 
AXgj.  Figure  13  shows  how  iXg^  is  defined  when  the  projectile 
is  rotated  through  an  angle  6.  Pg^  as  well  AXg^  depends  on  h 
and  0 . The  variation  of  both  A^g ^ and  Pg^  with  h and  6 were 
included  in  the  procedures  for  determining  the  overall  net 
forces  and  moments  on  the  projectile.  The  details  of  these 
procedures  are  found  in  Section  4.3. 

In  the  following  paragraphs,  the  equations  used  to  re- 
flect the  physics  of  the  three  fluid  dynamic  models  are 
discussed. 
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Figure  13.  Sketch  Showing  Definition  of  Shock  Interaction 

^^Sl  when  Projectile  is  Rotated  through 
an  Angle  0 . ^ 
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3 

4.2.1  Laminar  Couette  Flow  Model  (Re^<10  ) 

For  determining  the  viscous  shear  forces  and  ablation 
rate,  a local  two-dimensional  solution  to  Equations  (13)  and  (16) 
is  used.  Based  on  a thorough  study  of  gas  bearing  literature,  we 
concluded  in  Section  3.5  that  the  pressure  in  the  interface  could 
be  expected  to  remain  reasonably  constant  at  the  entrance  value 
over  the  entire  length  of  the  interface.  Some  pressure  variation 
can,  of  course,  be  expected  if  the  surface  of  the  interface  is 
not  parallel  to  the  rail  surface  and  if  the  rate  of  mass  addition 
to  the  interface  flow  due  to  ablation  reaches  of  the  same  order 
as  the  flow  of  air  pumped  into  the  interface  due  to  motion  of  the 
projectile.  Both  of  these  effects  are  included  in  determining 
the  normal  forces  on  the  projectile.  The  effect  of  variation 
of  the  gap  height  along  the  interface  is  taken  into  account  in 
determining  the  viscous  shear,  but  not  to  determine  the  ablation 
rate. 

The  skin  friction  coefficient  for  laminar  Couette  flow, 
c^,  based  upon  the  assiimption  that  there  is  roughly  equal  heat 
transfer  to  the  rail  and  to  the  projectile  is  given  by;“ 

Cj  = 2(1  + 0.0233  M^^)/Rej^  , (43) 

and  the  heat  transfer  rate,  q^,  is  given  by; 

^ = 0.25  CjP„  vj  , (44) 

where  c^  is  determined  from  Equation  (43)  using  the  average  gap 
height.  The  net  tangential  force  due  to  the  interface  is  given 
by; 


^tan  “ \ ? 


w 


where ; 
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•^00  oo 


(45) 


(46) 
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if  the  projectile  surface  is  parallel  to  the  rail.  If  the  pro- 
jectile surface  is  not  parallel  to  the  rail,  then  the  tangential 
force  is  determined  by  integrating  Equation  (42)  over  the 
length  of  the  interface  in  order  to  obtain  an  average  value  of 
Cj  which  is  then  substituted  into  Equation  (46) , in  order  to 

determine  f*..„. 

t:an 

The  normal  force  due  to  pressure  is  taken  to  be  the  value 
of  the  pressure  multiplied  by  the  area  given  by  W^B  plus  an 
additional  force  due  to  pressure  buildup  caused  by  ablation  and  a 
contribution  due  to  gas  bearing  effects  if  the  surfaces  are  not 
parallel.  The  additional  pressure  due  to  ablation  is  determined 
by  applying  the  simple  model  developed  in  Reference  1 for  the 
case  of  strong  ablation.  In  Reference  1,  the  average  pressure  in 
the  interface  due  to  ablation  is  determined  by  comparing  the  mass 
flow  at  the  exit  of  the  interface  to  the  mass  ablation  rate.  It 
was  assumed  in  Reference  1 that  the  pressure  would  buildup  in  the 
interface  such  that  the  density  of  the  ablated  gas  vapor  would 
become  large  enough  to  allow  the  full  ablation  mass  flow  to  pass 
through  the  exit  area  of  the  interface.  This  is  a grossly 
simplified  model,  but  was  considered  to  be  adequate  for  the  pre- 
sent purpose.  The  model  of  Reference  1 was  implemented  in  the 
program  by  fitting  a curve  to  Figure  3 of  Reference  1 to 
describe  a non-dimensional  mass  loss  quantity  (^  , as  a function 
of  projectile  velocity.  ♦ depends  upon  the  thermal  conductivity 
and  viscosity  of  the  ablated  vapor,  the  ablation  temperature  of 
the  projectile  material  as  well  as  the  projectile  velocity.  The 
pressure  in  the  interface  due  to  ablation  was  then  deter- 

mined from  Equation  11  of  Reference  1; 


abl 


*a  ^a  % 

h 

00  r 


(47) 


where  R is  the  gas  constant  of  the  ablated  vapor.  The  gas  dyna- 

cl 

mic  properties  for  the  ablated  vapor  were  taken  as  those  given  in 
Reference  1 for  polycarbonate  plastic. 

The  Increase  in  pressure  due  to  non-parallel  interface 
surfaces  is  taken  into  account  in  the  following  fashion.  First, 
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the  bearing  number  is  calculated  for  an  interface.  The  bearing 
number.  A,  is  defined  as; 


A = — - , 

V P4 

where  hj^^  is  the  minimum  gap  height  in  the  interface.  The  bearing 
number  is  an  important  dimensionless  quantity  which  arises  in  the 
hydrodynamic  theory  of  lubrication.  The  following  observation 
has  been  made  from  an  analysis  of  computer  solutions  to  the  com- 
pressible flow  version  of  the  Reynolds  Lubrication  Equation  pre- 
sented in  Reference  8 . For  a wide  variety  of  bearing  shapes  the 
non-dimensional  bearing  load,  defined  as: 


can  be  expressed  as; 


W = 


w 


(49) 


W = 0.006  X A (A  < 100)  , 


or; 


W = 0.6  (A  > 100)  , 


where  W is  the  normal  load  on  the  bearing  and  W is  the  average 
pressure  which  exists  over  the  bearing  surface,  over  and  above 
the  pressure  at  the  inlet  to  the  bearing.  Thus,  the  additional 
normal  force,  to  the  gas  bearing  effect  caused  by  non- 

parallel interface  surfaces  can  be  expressed  as; 


'lube 


= W 


A P 
*£  4 


(50) 


The  area  referred  to  in  Equation  (50)  as  represents 
the  area  over  which  significant  lubrication  effects  can  be  ex- 
pected. The  procedure  used  in  the  subroutine  INTER  to  evaluate 
this  area  is  to  first  consider  the  location  on  the  projectile  where 
the  minimum  gap  height  exists.  If  the  minimum  gap  height^  exists 
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at  the  rear  of  the  model,  the  entire  interface  area  (W^B)  is  used 
to  define  The  center  of  pressure  through  which  the  lubrica- 

tion normal  force  acts  in  this  case  is  taken  to  be  at  a distance 
of  0.7B  from  the  front  of  the  interface.^  If  the  minimum  gap 

height  is  found  to  be  at  the  front  of  the  interface,  then  the 

2 

area  is  defined  as  . For  this  situation,  the  front  part 
of  the  gas  bearing  surface  is  considered  to  be  the  entrance  region 
of  the  interface,  i.e.,  the  bearing  can  be  described  as  a tapered- 
flat  slider  bearing  where  the  taper  angle  is  equal  to  the  half- 
angle of  the  cone.  The  center  of  pressure  is  taken  to  be  at  the 
shoulder  of  the  cone  cylinder. 

4.2.2  Turbulent  Couette  Flow  Model  (10^<Rej <10^) 

In  this  model  fully  merged  turbulent  boundary  layers 
are  presumed  to  exist  at  the  entrance  region  upstrecun  of  the 
interface  region.  The  Reynolds  number  Re^^  is  sufficiently  large 
to  prevent  relaminarization  in  the  interface  itself  and  the  tur- 
bulent Ccuette  flow  model  of  Korkegi  and  Briggs^**  is  used.  The 
model  presented  in  Reference  10. is  developed  for  the  case  of  equal 
heat  transfer  rates  to  the  rail  and  to  the  projectile.  In 
Reference  11,  Korkegi  and  Briggs  developed  a more  general  model 
for  non-equal  heat  transfer  rates,  but  the  effect  of  the  non- 
equal heat  transfer  rates  does  not  significantly  affect  the  skin 
friction  (see  Figure  6 of  Reference  11) . In  addition,  a con- 
siderable savings  in  computational  time  is  realized  by  using  the 
formulation  of  skin  friction  for  equal  heat  transfer  rates.  The 
equation  developed  in  Reference  10  for  skin  friction  is  given  by; 


2ni/2 


(51) 


where ; 
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Equation  (51)  is  shown  as  a function  of  Mach  number  for 
various  values  of  Re^  in  Figure  14.  The  actual  evaluation  of 
Equation  (51)  in  the  computer  program  is  carried  out  in  sub- 
routine TURBCE. 

The  average  pressure  over  the  interface  will  be  a weak 
function  of  Re^^  due  to  the  fact  that  the  mean  velocity  distri- 
bution is  a fully  developed  turbulent  Couette  flow  is  such  that, 
near  the  center  of  the  interface , the  derivative  of  the  mean 
velocity  profile  with  respect  to  z is  zero.  Thus,  the  central 
part  of  the  flow  will  be  influenced  by  the  circumferential  pres- 
sure gradient  and  some  pressure  relief  can  be  expected.  In 
subroutine  INTER,  the  average  pressure  is  taken  to  be; 

P = P3  - (P3-P2)  [log(Rejj)  -3]/2  (52) 

The  center  of  pressure  of  this  average  pressure  is  presumed  to 
act  at  a distance  of  0.4B  from  the  front  of  the  interface  reflect- 
ing the  fact  that  the  pressure  distribution  is  such  that  a higher 
pressure  will  exist  over  the  front  portion  of  the  interface  than 
over  the  rear  portion.  The  heat  transfer  rate  is  determined  by 
substituting  from  Equation  (51)  into  Equation  (44) . 

4.2.3  Turbulent  Boundary  Layer  Model  (Re^^>10^) 

The  basic  assumption  underlying  models  for  turbulent 
boundary  layers  is  that  significant  circumferential  pressure 
relief  will  occur  due  to  the  presence  of  an  inviscid  core  region 
over  the  initial  portion  of  the  interface  zone  and  that  the 

amount  of  pressure  relief  will  increase  as  Re.  becomes  greater 

5 5 “ 

than  10  . The  criterion  of  using  Rej^  = 10  as  the  Reynolds  num- 
ber at  which  this  pressure  relief  will  become  significant  has 
not  been  rigorously  established,  but  rather  represents  our  judge- 
ment made  on  the  basis  of  interface  pressure  distributions 
obtained  by  Korkegi  and  Briggs  for  hypersonic  slipper  bearings 
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Figure  14.  Skin  Friction  Coefficient  Versus  Mach  Number  for 
Turbulent  Couette  Flow  for  Equal  Heat  Transfer 
Rates  at  Rail  and  Projectile  Surfaces. 
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tempered  with  the  fact  that  pressure  relief  will  occur  more 
rapidly  in  the  rail/projectile  problem  than  in  the  hypersonic 
slipper  problem. 

The  skin  friction  for  the  turbulent  boundary  layer  is 
determined  from  the  Karman-Schoenher  implicit  formulation* 
given  by: 


0.242  = log  (Re^^  C^)  , 


(53) 


where  is  the  average  skin  friction  coefficient  and  the 
Reynolds  number  is  based  on  conditions  evaluated  at  the  pressure 
existing  at  the  sides  of  the  interface  and  at  a reference  temper- 
ature corresponding  to  the  Sommer  and  Short  reference  enthalpy 

it 

defined  by  H /H^  given  by; 


H*  ”w  “aw 

= 0.36  + 0.45  ^ + 0.19  ^ 

“e  ”e  “e 


(54) 


A blowing  correction  factor  is  applied  to  the  skin  friction 
Cj.  This  blowing  factor  is  given  by;“* 


»c  - 1 + S.57 


(55) 


is  computed  in  subroutine  TURBFP.  The  tangential  force  is 
then  determined  from  the  following  equation; 


2 


^tan  - 1/=^  “e  <=£  "r  '‘a  ®c 


(56) 


and  heat  transfer  rate  is  determined  from; 


‘Jw  = Pe  ^e  2^  ^“aw-“wJ  "r  '^a  “c 


(57) 


The  average  pressure  is  given  by; 


P = P2  - (P2''^5^  [log(Rej^)  - 5]/5  . 


(58) 
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The  center  of  pressure  is  assumed  to  act  at  the  midpoint  of  the 
afterbody.  This  completes  the  discussion  of  the  present  version 
of  subroutine  INTER.  It  should  be  mentioned  that  a simplified 
interface  model,  discussed  next,  was  developed  early  in  the  pre- 
sent effort  in  order  to  check  out  the  Projectile  Dynamics  Com- 
puter Program. 


4.2.4  Simplified  Interface  Subroutine 


A simple  model  which  estimated  at  least  the  order  of 
magnitude  of  the  forces  and  moments  on  the  projectile  due  to  the 
interface  was  developed  at  the  beginning  of  the  program.  From 
Equations  (43)  and  (46)  it  is  observed  that  for  high  projectile 
Mach  numbers  and  small  gap  heights,  the  shear  stress  can  be 
approximated  by; 


T 


W 


C U' 


(59) 


where  c is  a constant  of  proportionality.  For  small  inclinations 
of  the  projectile  relative  to  the  rail,  the  pressure  due  to  gas 
bearing  effects  and  lubrication  effects  is  inversely  proportional 
to  the  square  of  the  gap  height,  i.e.; 


P = e/h^  , 


(60) 


where  e is  another  constant  of  proportionality.  These  two  equa- 
tions were  used  to  construct  a simple  interface  model.  The 
details  of  this  model  are  presented  in  Appendix  A.  This  model 
has  some  potential  usefullness  insofar  as  it  is  very  easy  to  use 
requiring  only  that  the  two  constants  c and  e be  specified. 

This  model  could  conceivably  be  used  to  represent  the 

forces  required  to  stabilize  an  unstable  model  since  no  actual 

2 

physical  collision  with  the  rail  can  occur  because  of  the  1/h 
behavior  of  the  pressure.  Thus,  the  maximum  normal  force  required 
to  prevent  physical  contact  between  the  projectile  and  the  rail 
can  be  estimated  with  the  simple  model. 
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4.3  INFLUENCE  OF  GAP  HEIGHT  AND  MODEL  ROTATION  ON  SHOCK 
INTERACTION  LENGTH  AND  PRESSURE 

The  shock  interaction  length  is  obtained  by  noting 

the  following  geometric  relations  which  apply  to  Figure  13. 


e = h - ALNSI  sin  6 (61) 

a = e/tan  (0  - 0)  (62) 

s 

c = h/tan  (0^  - 0)  (63) 

d = ALNSI  cos  0 - a - c (64) 


Then  the  application  of  the  law  of  sines  to  the  triangle  defined 
by  the  angles  Qg  “ and  180®  - (0^.  + 0^,  - 20)  with 

two  of  the  sides  being  d and  yields; 

sin  (0^  - 0) 

^^SI  ^ ^ sin  (0J.  + - 20)  * 

Note  that  this  result  is  based  on  the  assumption  that  the  bow 
shock  rotates  with  the  projectile. 

The  length  ALNSI  is  either  computed  or  read  in  as  input 
data  (see  Appendix  B) . 

The  pressures  at  states  2,  3,  or  4 in  Figure  13  depend 
on  0.  As  can  be  seen  from  Figure  13,  if  the  nose  tip  rotates 
toward  the  bottom  rail,  the  lower  conical  surface  of  the  pro- 
jectile is  inclined  at  a smaller  angle,  (0^  - 0),  to  the  on- 
coming flow  and  the  strength  of  the  lower  portion  of  the  bow 
shock  will  be  less  than  the  shock  strength  when  0=0. 
Correspondingly,  the  upper  conical  surface  will  be  inclined  at 

angle  (0  + 0)  to  the  oncoming  flow  and  the  strength  of  the 

c 

upper,  portion  of  the  bow  shock  will  be  greater  than  the  0=0 
value. 

For  small  turning  angles,  0,  the  pressure  change  can  be 
expressed  as; 

AP  = P”  . ..  (±0)  , (66) 

Jv?  - 1 
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where  the  plus  sign  is  used  when  the  flow  is  compressed.  Using 
this  result  the  pressure  at  state  2 with  the  projectile  rotated 
can  be  approximated  in  terms  of  the  pressure  at  state  2 when 
6 = 0 in  the  following  manner; 


^2(0?^O)  ~ ^2(0=0) 


+ (±0)  , 


(67) 


where  the  plus  sign  is  used  if  the  surface  becomes  more  steeply 
inclined  to  the  oncoming  flow  through  projectile  rotation  and 
the  negative  sign  is  used  if  the  surface  becomes  less  steeply 
inclined.  Similar  reasoning  leads  to  the  following  expressions 
for  the  pressures  at  states  3 and  4; 


^3(0?«O) 

^4(0;^O) 


^2(0;^O) 

^3(6?^0) 


+ 


(±0)  , 


(±20) . 


(68) 

(69) 
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SECTION  V 

COMPUTER  ANALYSIS  OF  PROJECTILE  DYNAMICS 

The  interface  models  described  in  Section  IV  account 
for  the  effects  of  the  longitudinal  motion  of  the  projectile. 

The  projectile  can  also  translate  perpendicular  to  its  axis 
and  can  rotate  about  axes  normal  to  its  principal  longitudinal 
axis  (pitch  and  yaw) . These  transverse  motions  have  a direct 
effect  on  the  interface  regions  and  therefore  contribute  to  the 
dynamic  behavior  of  the  interfaces  with  the  rails.  A computer 
program  describing  projectile  dynamics  was  developed  to  under- 
stand these  effects  which  includes  descriptions  of  longitudinal, 
translational,  and  rotational  dynamics  of  the  projectile.  The 
interface  behavior  is  treated  as  coupling  between  projectile  dy- 
namics and  interface  behavior. 

This  model  of  projectile  dynamics  thus  provides  a complete 
description  of  the  operation  of  a guide-rail  system  with  hyper- 
velocity  projectiles.  All  l«nown  components  of  the  drag  induced 
on  the  projectile  are  computed  as  are  the  wear  and  wear  rates 
generated  on  the  cylindrical  surface  of  the  projectile  by  inter- 
action with  the  guide-rails.  Finally,  the  model  predicts  the 
stability  of  the  projectile  in  flight  along  the  rails  and  establishes 
whether  or  not  the  projectile  contacts  the  rails  in  flight.  The 
remainder  of  this  section  describes  the  computational  architecture 
of  the  projectile  dynamics  program  and  presents  results  from  a 
limited  initial  comparison  between  computed  results  and  experi- 
mental measurements  from  operation  of  the  Tracks  K and  G at  AEDC . 

Some  results  from  a scaling  study  are  presented  at  the  end  of 
this  section. 
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5.1  OVERVIEW  OF  PROJECTILE  DYNAMICS  PROGRAM 

A generalized  flow  chart  of  the  projectile  dynamics 
program  is  presented  in  Figure  15.  Operation  of  the  program  is 
initiated  by  inputting  the  following  parameters  to  an  initialization 
subroutine:  dimensions  of  the  rail  system  and  projectile;  material 
characteristics  of  the  projectile;  atmospheric  conditions  in  the 
ballistic  range;  aerodynamic  coefficients  of  the  projectile; 
irregularities  of  the  rail  system;  initial  position,  orientation, 
and  velocity  of  the  projectile;  and  the  variety  of  commands  used 
to  select  various  computational  options  within  the  remainder  of 
the  program  structure.  The  initialization  subroutine  then 
calculates  a variety  of  parameters  used  for  subsequent  compu- 
tations such  as  initial  gap  height  between  the  projectile  and  each 
of  the  rails,  and  initial  orientation  of  the  model  and  its  vel- 
ocity vector  with  respect  to  the  upstream  end  of  the  rail  system. 

The  main  section  of  the  program  is  an  executive  routine 
which  coordinates  the  operation  of  an  array  of  subroutines  that 
accomplish  the  actual  computations.  Aerodynamics  subroutines 
(called  "DRAG"  and  "DRAGCO")  are  available  which  translate  the  aero- 
dynamic coefficients  of  the  projectile  together  with  the  model 
velocity  and  range  atmospheric  parameters  into  aerodynamic  forces 
and  moments  applied  to  the  model.  Provisions  are  available  within 
subroutine  DRAGCO  to  treat  aerodynamic  coefficients  as  functions 
of  Mach  number  and  other  independent  variables.  The  aerodynamic 
subroutine  DRAGCO  also  accounts  specifically  for  shock  inter- 
ference effects  due  to  the  presence  of  the  rail  system.  These 
routines  are  described  generally  in  Section  IV  and  specifically 
in  Appendix  B. 

The  projectile  dynamics  computer  program  has  available 
two  separate  subroutines  for  analyzing  the  interface  regions  be- 
tween the  model  and  the  guide-rails.  Both  are  called  "INTER"  and 
are  described  generally  in  Section  IV  and  are  specified  in 
Appendix  B.  The  first  is  an  extremely  simplified  model  con- 
taining virtually  no  physics  that  was  constructed  to  provide 
credible  output  data  from  the  interface  regions  with  minimal 
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Figure  15.  Overview  of  Dynamic  Rail/Projectile  Modeling  Program. 
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computational  effort.  This  subroutine  was  used  during  develop- 
ment of  the  remainder  of  the  overall  computer  program  and  is  still 
useful  for  debugging  tasks.  The  second  interface  subroutine 
contains  all  of  the  physics  of  the  projectile  rail/interface 
region  that  the  authors  currently  understand,  and  is  the  model 
that  should  be  used  for  making  actual  predictions  of  rail/projectile 
performance.  The  complete  interface  subroutine  computes  all 
known  forces  arising  within  the  interface  region  between  the 
projectile  and  the  rails  and  also  computes  wear  of  the  projectile. 

The  subroutine  that  employs  rigid  body  equations  of  motion 
to  describe  the  position,  velocity,  orientation,  and  the  pitch 
rate  of  the  model  at  any  instant  after  the  model  has  engaged 
the  rail  system  is  called  "RK"  in  the  main  computer  program. 

This  subroutine  accepts  inputs  from  the  other  subroutines 
describing  forces  and  torques  applied  to  the  projectile  and 
outputs  to  appropriate  subroutines  projectile  positional  in- 
formation required  for  their  operation.  Subroutine  RK  is  des- 
cribed in  a subsequent  paragraph  of  this  section  and  is  specified 
in  detail  in  Appendix  B. 

The  rail  system  subroutines  called  "ORIENT"  and  "CONFIG" 
contain  equations  of  the  rail  surfaces  relative  to  an  earth- 
fixed  coordinate  system.  These  routines  are  described  in  Appen- 
dix B.  CONFIG  describes  irregularities  of  the  rail  system, 
such  as  steps,  bends,  convergences/divergences,  etc.  The  in- 
stantaneous separation  distances  between  the  model  and  each  of 
the  rails  are  computed  in  ORIENT  based  on:  the  position  of  the 
rails  as  defined  by  subroutine  CONFIG;  the  position  of  the  pro- 
jectile as  defined  by  the  equations-of-motion  subroutine;  and 
the  amount  of  wear  the  projectile  has  experienced  as  defined  by 
the  interface  subroutine.  The  projectile  position,  orientation, 
and  wear  are  all  outputted  at  time  or  space  intervals  specified 
as  input.  Computation  cycles  of  the  main  computer  routine  are 
repeated  until  the  projectile  either  reaches  the  end  of  the  rail 
system  or  exceeds  a predefined  performance  envelope  which  leads 
to  the  premature  termination  of  its  flight. 
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5.2  EQUATIONS  OF  MOTION  SUBROUTINE- 

The  three-dimensional  equations  of  motion  of  the 
projectile  are; 


M ic  = - cos  Y 
p A 


Dj  cos  0 + sin  Y + Lj  sin  6, 


(70) 


Mz  = -Mg-D-  sin  y - Dt  sin  6 + cos  y + L,.  cos  0,  (71) 

D A I A I 


I 9 = T„  + T_, 
y A I 


(72) 


The  x-z  coordinate  system  is  an  earth-fixed  system  with  the 
x-axis  perpendicular  to  the  earth's  gravitational  vector  at  the 
front  end  of  the  rail  system.  The  z-axis  is  parallel  to  the 
gravitational  vector  and  is  positive  upward.  Figure  16  illustrates 
the  directions  of  the  forces  acting  on  the  projectile.  The 
aerodynamic  drag  is  in  the  opposite  direction  to  the  velocity 


Figure  16.  Forces  Acting  on  the  Projectile. 
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vector  and  the  aerodynamic  lift  is  perpendicular  to  the  velocity 
vector.  The  frictional  force  arising  from  the  projectile-rail 
interfaces  is  parallel  to  the  projectile  body  and  the  lift 
arising  from  the  interfaces  is  perpendicular  to  the  projectile 
body.  The  angle  y is  defined  as; 

^-14 

Y = tan  v ■ 

X 

(73) 

acting 

The  magnitudes  of  the  aerodynamic  forces 
on  the  projectile  are  given  by; 

and  torque 

= 1/2  p SV^  Cp, 

(74) 

= 1/2  p SV^  Cj^, 

(75) 

Ta  - 1/2  P SCV^  C^. 

(76) 

The  magnitude  of  the  velocity  vector  is  given  by  the  equation; 

, o 1/2 

V = + 4^)  . (77) 


The  aerodynamic  coefficients  are  given  by  the  equations; 


° *^0  a SI 


C=c  +C  a + ^C  +^C+C^. 

m m^j  m^  2V  m^  2V  m^  mgj 


The  angle  of  attack  and  its  time  derivative  are  given  by; 


(78) 

(79) 

(80) 


a = 9 - Yf 


(81) 
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and; 

d = q + g/V  cos  y - pVSCj^/(2Mp)  . (82) 

The  above  equations  define  the  motion  of  the 

projectile.  The  values  of  the  various  coefficients  in  the  equations 

for  C„,  C, , and  C are  computed  in  subroutine  DRAG.  The  values 
D L in 

for  Dj,  Lj,  and  T^  are  computed  in  subroutine  INTER.  The  forces 
and  moments  in  the  x,z,  and  6 directions  are  computed  in  sub- 
routine F from  information  determined  by  subroutines  DRAG  and 
INTER.  The  equations  of  motion  are  integrated  with  a fourth- 
order  Runge-Kutta  scheme  in  subroutine  RK. 

5.3  ELASTIC  DEFORMATION  SUBROUTINE 

As  was  mentioned  in  Section  3.6,  only  rigid  body  behavior 
of  the  projectile  was  initially  considered.  This  approach  was 
found  to  be  inadequate  when  the  rigid  body  computer  model  was 
used  to  simulate  Track  K experiments.  It  was  found  that  the 
slightest  perturbation  to  the  projectile  would  always  result 
in  the  model  colliding  with  a rail.  The  pressures  developed  in 
the  interface  due  to  gas  bearing  and  ablation  effects  were  too 
small  to  keep  the  projectile  cylindrical  surface  from  con- 
tacting the  rail.  In  order  to  estimate  the  forces  and  pressures 
required  to  prevent  projectile/rail  contact  a linear  elastic 
analysis  of  projectile  deformation  was  conducted.  This  analysis 
showed  that  the  forces  required  to  deform  the  projectile  and 
prevent  contact  were  very  large.  If  this  force  arises  from  gas 
bearing  effects,  the  interfacial  pressures  necessary  to  account 
for  such  high  forces  imply  very  small  gap  heights  and  corres- 
pondingly small  Reynolds  numbers.  The  gap  heights  required  are 
of  the  order  of  the  surface  asperities.  When  the  importance  of 
elastic  deformation  was  realized,  an  elastic  collision  subroutine, 
ELASTC  , was  added  to  the  program.  This  subroutine  is  based 
on  a simple  linear  elastic  relation  of  the  form; 
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(83) 


where  is  the  local  force  per  unit  area  due  to  a local 

elastic  deformation  6 and  k is  an  elastic  coefficient. 

e*  e 

Subroutine  ELASTC  is  called  whenever  subroutine  ORIENT  pro- 
duces a negative  separation  distance  between  the  projectile  and 
a rail,  based  upon  a rigid  body  calculation.  A linear  relation 
between  deformation  required  to  reduce  the  negative  separation 
distance  to  zero  and  distance  along  the  projectile  surface  in  the 
x-direction  is  presumed.  The  normal  force  required  to  deform 
the  projectile  surface  so  that  the  projectile  surface  does  not 
pass  outside  the  boundary  of  the  rail  is  then  determined  by  in- 
tegrating Equation  (83)  over  the  region  where  negative  separation 
distances  would  occur  if  the  projectile  were  perfectly  rigid. 

The  equations  for  determining  the  normal  forces  and  moments  due 
to  elastic  collisions  are  developed  in  Appendix  C. 

Whenever  a collision  occurs  with  a rail,  an  additional 
mechanism  for  heat  transfer  exists  due  to  the  sliding  friction. 
The  sliding  friction  force,  is  given  by; 


sf 


(84) 


where  F is  the  total  normal  force  due  to  elastic  deformation, 
n 

The  energy  dissipation  rate  due  to  the  sliding  friction  is  given 
by; 


^sf 


= F „U 
sf  “ 


(85) 


which  is  taken  to  be  the  rate  of  energy  transfer  to  the 
projectile  surface.  The  ablation  rate  due  to  sliding  friction 
is  then; 


m 


sf 


= u 


U /L. 
n “ 


(86) 
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5.4  METHOD  USED  TO  CALCULATE  WEAR  ON  A PROJECTILE 

Mass  is  removed  from  the  projectile  through  two 
mechanisms,  ablation  due  to  aerodynamic  friction  and  ablation 
due  to  sliding  friction.  Neither  of  these  two  mechanisms  lead 
to  a uniform  surface  recession  rate  over  the  contact  area  be- 
tween a rail  and  the  projectile.  A simplified  method  was  used 
to  approximate  the  change  in  shape  of  the  wear  surfaces  of  the 
projectile.  The  top  and  bottom  wear  surfaces  of  the  pro- 
jectile are  each  approximated  by  two  straight  line  segments  as 
shown  in  Figure  17.  These  straight  line  segments  are  defined 
by  tracking  the  surface  recession  at  three  points  on  each  of 
the  top  and  bottom  wear  segments:  the  front;  the  midpoint;  and 
the  rear  of  the  cylindrical  portion  of  the  projectile.  The  wear 
surface  at  the  side  rail  contact  with  the  projectile  is  approxi- 
mated by  a single  straight  line  segment  defined  by  tracking 
the  surface  recession  at  the  front  and  rear  of  the  cylindrical 
portion  of  the  projectile.  The  two  side  rail  contact  surfaces  on 


Figure  17.  Sketch  Showing  How  the  Instantaneous  Shape  of  the 
Projectile  is  Defined  to  Account  for  Wear. 
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the  projectile  are  assiimed  to  wear  at  the  same  rate  whereas  the 
top  and  bottom  wear  surfaces  will  not,  in  general,  wear  symmetrically. 
Thus,  the  surface  recession  rate  is  computed  at  eight  different 
points  on  the  projectile. 

Note  in  Figure  17  that  the  front  segments  of  the  top  and 
bottom  wear  surfaces  extend  onto  the  conical  forebody.  Track 
G experience  with  shots  at  high  projectile  velocities  (5000  m/s 
and  higher)  at  range  pressures  of  0.5  atm  or  greater  showed  large 
amounts  of  wear.  When  attempts  were  made  to  simulate  these 
shots  with  an  early  version  of  the  computer  program  which  did 
not  consider  the  extension  of  the  wear  surface  onto  the  conical 
forebody,  the  predicted  projectile  wear  was  so  high  that  the 
projectile  could  have  tumbled  down  the  range.  The  increase  in 
the  length  of  the  wear  surface  as  the  shoulder  of  the  cylindrical 
portion  of  projectile  wears  is  called  "increase  in  wheelbase  due 
to  projectile  wear."  In  Figure  18  this  increase  in  wheelbase 
due  to  wear  is  denoted  by  AB.  The  radial  location  of  the  front 


Figure  18.  Sketch  Showing  Increase  in  Wheelbase  Due  to  Wear,  AB. 
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edge  of  this  extended  wear  surface  on  the  top  of  the  projectile 
is  denoted  by  r^^  in  Figure  18.  The  axial  position  of  r^^ 
and  its  counterpart  on  the  bottom  wear  surface,  rj^^  are  contin- 
uously tracked  in  the  computer  program.  In  particular,  this  point 
is  used  as  one  of  the  four  points  on  each  of  the  top  and  bottom 
wear  surfaces  in  subroutine  ORIENT  where  the  separation  dis- 
tance between  the  projectile  and  the  rail  is  computed.  Likewise 
the  normal  force  due  to  contact  with  the  rail  is  computed  at  this 
point  of  the  separation  distance  at  this  point  becomes  negative. 
The  increase  in  wheelbase  due  to  wear  serves  to  increase  the 
moment  arm  of  the  elastic  deformation  normal  force  when  contact 
with  the  rail  occurs  at  the  front  of  the  projectile.  Thus, 
the  magnitude  of  the  normal  force  and  wear  experienced  by  the 
projectile  during  such  a collision  with  the  rail  is  less  than  the 
normal  force  and  wear  computed  when  the  increase  in  wheelbase  due 
to  wear  is  not  included  in  the  analysis. 

5.5  ■ COMPARISON  OF  COMPUTED  RESULTS  WITH  LIMITED  EXPERIMENTAL 

DATA 

Limited  data  for  evaluating  the  predicted  results  have 
become  available  recently  as  a result  of  both  dynamic  measurements 
made  during  range  firings  at  Tracks  K and  G and  observations  of 
models  recovered  after  negotiating  the  hypervelocity  tracks. 

5.5.1  Comparisons  with  Track  K Results 

Track  K at  AEDC  was  outfitted  with  a set  of  hypervelocity 
guide  rails  approximately  30  m long.  This  system  became  oper- 
ational and  produced  significant  amounts  of  experimental  data 
during  the  time  that  the  study  described  in  this  report  was 
being  conducted.  Several  firings  of  projectiles  with  30“  and 
45“  half-angle  cone  sections  were  carried  out  where  good  quality 
data  was  achieved.  More  recently,  a 10“  half-angle  cone  pro- 
jectile was  launched  successfully  in  Track  K.  No  wear  of  the 

cylindrical  surfaces  of  launched  models  caused  by  interactions 
with  the  guide  rails  was  observed  experimentally  in  Track  K. 
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This  result  is  in  agreement  with  our  predictions  that  wear  during 
30  m of  travel  would  be  slight.  Unfortunately,  this  verification 
shed  relatively  little  light  upon  the  question  of  the  validity  of 
our  analysis. 

An  indirect,  but  potentially  important  evaluation  was 
carried  out  involving  the  total  drag  experienced  by  hypervelocity 
projectiles  moving  along  guide-rails.  The  Track  K data  indicate 
that  total  projectile  drag  was  30  percent  to  50  percent  greater 
than  would  have  occurred  had  the  projectiles  been  launched  into 
an  open  range.  The  computer  model  described  in  this  section  is 
capable  of  predicting  total  drag  forces  on  projectiles  flying 
along  hypervelocity  guide  rails,  so  a basis  of  comparison  was 
possible.  Shot  No.  2968  at  the  AEDC  Track  K involved  launching 
a 30®  half-angle  cone  at  a velocity  just  over  3.87  km/s.  The 
remainder  of  the  information  known  to  the  authors  about  Shot  No. 
2968  is  presented  in  Table  4.  The  total  drag  coefficient  as 
determined  from  velocity  reduction  data  was  0^=0.788.  This 
shot  was  used  as  the  basis  for  a computational  run  of  the 
.Projectile  Dynamics  Program  from  which  a total  drag  coefficient 
of  0.789  was  computed.  The  close  agreement  between  the  measured 
and  predicted  value  (just  over  0.1  percent  deviation)  is  almost 
certainly  fortuitous,  although  the  comparison  gives  good  reason 
to  hope  that  at  least  the  aerodynamic  phase  of  the  program 
describes  practical  situations  accurately. 

Also  shown  in  Table  4 are  results  from  a single  firing 
of  a 10®  half-angle  cone  at  a velocity  of  2.35  km/s.  A total 
drag  coefficient  of  0.199  was  measured,  while  the  results  of  a 
computation  using  identical  input  parameters  predict  a drag 
coefficient  of  0.170.  This  value  for  total  drag  was  the  result 
of  allowing  the  reflected  weak  shock  from  the  rail  to  re-reflect 
from  the  cone  once.  The  predicted  value  for  the  total  drag  co- 
efficient was  0.140  when  the  reflected  shock  from  the  rail  was 
not  allowed  to  re-reflect  from  the  conical  surface.  When  the 
rail  reflected  shock  was  presumed  to  be  a strong  oblique  shock, 
a total  drag  coefficient  of  0.633  was  obtained  computationally. 
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TABLE  4.  COMPARISON  OF  COMPUTED  TOTAL  DRAG  COEFFICIENTS 
WITH  EXPERIMENTAL  DATA  FROM  TRACK  K AT  T^DC. 


30“  CONE  PROJECTILE 


Shot  No.  2968 
Projectile  Type 

Projectile  Diameter  Before 
Launch 

Projectile  Diameter  After  Launch 

Projectile  Weight  Before  Launch 

Projectile  Velocity 

Range  Pressure 

Range  Temperature 

Measured  Total  Drag  Coefficient 

Predicted  Total  Drag  Coefficient 


30“  half-angle  cone,  Model 
No.  4927 

1.617"  (41.07  mm) 

1.600  to  1.604"  (40.64  to  40.74  nun) 
121.61  g 

12,317  ft/s  (3.755  km/sec) 

103.7  torr  (1.36  x 10^  pa) 

66“F  (18.9“C) 

0.788 

0.789 


10“  CONE  PROJECTILE 


Shot  No.  3014 


Projectile  Type 

Projectile  Bluntness  Ratio 

Projectile  Diameter  Before  Launch 

Projectile  Diameter  After  Launch 

Projectile  Weight  Before  Launch 

Projectile  Velocity 

Range  Pressure 

Range  Temperature 

Measured  Total  Drag  Coefficient 

Predicted  Total  Drag  Coefficients: 


10“  half-angle  cone 
0.167 

1.625"  (nom)  (41.3  mm) 


152.6  gm 

7700  ft/s  (2.348  km/sec) 
711.8  torr  (9.36  x 10^  pa) 
70“F  (nom)  (21.1“C) 

0.199  + .003 


weak  reflected  shock  reaches  cone  surface:  0.170 
weak  reflected  shock  doesn't  reach  cone  surface:  0.140 
strong  reflected  shock:  0.633 
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We  feel  that  the  later  computational  result  is  especially  sig- 
nificant since  the  large  deviation  from  the  measured  total  drag 
coefficient  demonstrates  that  shock  waves  reflected  from  the 
rails  by  slender  cone  projectiles  are  not  strong.  The  viscous 
drag  computationally  predicted  to  arise  from  the  rail  interface 
region  for  the  10“  half-angle  cone  was  only  0.6  percent  of  the 
total  drag,  and  therefore,  does  not  figure  to  any  significant 
extent  in  the  difference  between  predicted  and  measured  results. 

A second  factor  of  considerable  importance  became  obvious 
during  the  computation  of  the  10“  cone  trajectories.  The  pro- 
jectile is  highly  unstable  aerodynamically  and  our  analysis  which 
included  deformation  considerations  indicates  that  the  projectile 
will  contact  the  rails.  Careful  consideration  of  the  predicted 
results  lead  us  to  believe  that  this  prediction  is  correct,  and 
therefore,  that  the  projectile  fired  in  Track  K indeed  did  contact 
the  rails  repeatedly  during  its  passage  along  them.  No  unusual 
gouging  or  wear  were  noted  after  the  firing  either  on  the  recovered 
projectile  or  the  rails.  These  considerations  lead  us  to  believe 
that  our  initial  ass\imption  that  gouging  always  occurs  when  ^bhe 
projectile  touches  the  rails  during  passage  is  an  oversimplification. 
Clearly,  more  work  is  needed  to  establish  a proper  gouging  criterion. 
In  summary,  we  have  strong  indications  that  projectiles  can  contact 
rails  without  undue  wear  or  gouging  at  least  under  some  circum- 
stances. Computed  values  of  total  drag  coefficients  correlate 
well  with  measurements  for  flights  of  projectiles  with  30“  half- 
angle cones . Deviations  between  predicted  and  measured  total  drag 
coefficients  are  somewhat  larger  for  projectiles  with  slender  cone 
forebodies.  These  deviations  are  probably  caused  by  either  under- 
estimation of  the  conical  area  of  the  projectile  affected  by  shock 
waves  reflected  from  the  rails,  or  by  drag  increases  caused  by 
projectile/rail  contact.  Finally,  comparisons  between  computational 
and  experimental  results  for  the  10“  cone  flights  strongly  indicate 
that  shock  waves  reflected  from  the  rails  are  weak  and  indicate 
somewhat  less  certainly  that  reflected  shock  waves  from  the  rails 
actually  re-ref lect  from  the  conical  surface. 
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5.5.2  Comparison  to  Track  G Results 

Track  G projectiles  have  typically  been  composite 
structures.  The  cylindrical  portion,  including  the  primary 
wear  surfaces,  have  been  made  from  polycrabonate . An  internal 
aluminum  strut  supports  the  nose  tip  and  provides  a conical 
surface  to  which  a carbon-epoxy  heat  shield  is  bonded.  Some 
projectiles  have  forebodies  which  are  essentially  hemisphere- 
cylinder-conical  configurations  while  others  have  forebodies 
which  are  essentially  blunted  cones.  Table  5 lists  wear  mea- 
surements and  instantaneous  angular  rotation  measured  from  flash 
radiographs  or  laser  photographs  at  the  last  observation  point 
where  photographs  were  available.  The  half-angle  of  the  conical 
portion  of  the  forebody  on  all  of  the  projectiles  listed  in 
Table  5 was  25  degrees. 


TABLE  5.  WEAR  AND  ANGULAR  ROTATION  OBSERVATIONS  FOR  SELECTED 
TRACK  G SHOTS 


Ifelocity 

Shot  No.  (m/s) 

Range 

Pressure 

(torr) 

Max.  Groove 
Depth* 

(mm) 

Min.  Groove 
Depth* 

(mm) 

Angle  of 

Rotation***  Type  of 
(degrees)  Model 

4836 

2870 

233 

0.89 

0.41 

-1.5 

Erosion 

4882 

4880 

570 

3.30 

1.27 

8 

Transition 

4885 

4880 

100 

2.06 

1.37 

2.2 

TCNT** 

4886 

4880 

150 

2.29 

1.02 

Erosion 

4880 

4880 

100 

1.78 

0.89 

>1 

Erosion 

4884 

4880 

570 

3.18 

1.27 

5 

Transition 

*Based 

on  post-launch  diameter 

of  62.5  mm 

(0.51  mm  wear 

assumed  in 

launcher) . 

**TCNT — transpiration  cooled  nose  tip. 


***Angle  of  rotation  at  last  x-radiograph  or  laser  photography  station — 
not  necessarily  the  maximum  angle  of  rotation  of  the  projectile. 
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A number  of  computer  simulation  runs  were  made  for  a 
Track  G projectile.  The  purpose  of  these  runs  was  primarily 
to  determine  values  of  the  sliding  friction  coefficient  il  and 
the  heat  of  ablation  L which  would  yield  good  correlations 
between  computed  and  measured  wear,  ii  was  varied  over  a range 
from  0.0001  to  0.02  and  L was  varied  from  3.8  x 10®  joules/kg 
to  6.4  X 10®  joules/kg.  Typical  results  of  these  computations 
are  shown  in  Table  6.  The  value  of  the  elastic  coefficient  k^  = 

2.6  X 10^^  N/m^  ('\'0.78  x 10®  Ib^/in^)  was  obtained  from  Dr.  J. 

Hicks  of  ARO/VKF.  He  used  a finite  element  structural  analysis 
technique  to  determine  the  stiffness  of  a typical  Track  G com- 
posite model.  His  analysis  showed  that  the  elastic  coefficient 
for  Track  G models  could  be  represented  by; 

k = 2.05  X 10_^  (N/m^),  (87) 

e 

where  is  the  width  of  the  contact  surface  (rail  width  in  meters) . 
The  rail  width  of  Track  G rails  is  8 mm  (0.31  in.).  The  geometric 
parameters  of  the  projectile  used  in  the  computer  simulation  runs 
were  as  follows:  a = 25®,  = 0.3,  d^  = 62.5  mm  (2.46  in.), 

B = 63.5  mm  (2.50  in.).  An  initial  gap  height  of  0.5  mm  (0.020  in.) 
was  used  reflecting  Track  K and  Track  G experience  where  about 
0.5  mm  of  wear  was  found  to  have  occurred  in  the  launcher . Based 
on  the  results  shown  in  Table  6 it  was  concluded  that  u = 0.0075 
and  L = 4.2  X 10®  joules/kg  produced  acceptable  agreement  between 
computed  and  measured  wear  for  Track  G projectiles. 

An  important  conclusion  that  can  be  drawn  from  the  results 
presented  in  Table  6 is  that  the  maximum  wear  experienced  by  the 
projectile  decreases  as  the  wavelength  of  sinusoidal  rails  de- 
creases, all  other  parameters  being  held  constant.  Note  that  the 
maximum  wear  for  a wavelength  of  24.4  m (80  ft)  is  nearly  the 
same  as  for  the  corresponding  straight  rail  case. 
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TABLE  6. 


COMPUTED  WEAR  AND  MAXIMUM  ANGLE  OF  ROTATION  FOR  TRACK 
G MODEL  FOR  VARIOUS  VALUES  OF  HEAT  OF  ABLATION  AND  ^ 
COEFFICIENT  OF  SLIDING  FRICTION.  k = 2.6  x lO-*--*-  N/m 


\felocity 

(m/s) 

Range 

Pressure  L 

(torr)  MJ/kg 

Max.  Groove 
Depth 

VI  (mm) 

Min.  Groove 
Depth 
(mm) 

Max.  Angle 
of  Rotation 
(degrees) 

Straight  Rails 

2870 

233 

4.2 

0.0005 

0.26 

0.23 

-1.3 

4880 

100 

4.2 

0.0005 

0.48 

0.43 

2.0 

4880 

150 

4.2 

0.0005 

0.90 

0.58 

2.5 

4880 

150 

6.4 

• 0.0005 

0.58 

0.45 

-2.0 

4880 

570 

4.2 

0.0005 

4.52 

1.50 

-6.6 

4880 

570 

3.8 

0.0005 

5.00 

1.61 

-6.8 

4880 

760 

4.2 

0.0005 

5.44 

1.76 

-6.8 

Sinusoidal  Rails, 

Wavelength 

♦H 

II 

(20  ft).  Amplitude  = 0.19  mm  (0.0075  in.! 

4880 

570 

4.2 

0.0010 

2.02 

1.50 

1.8 

Sinusoidal  Rails, 

Wavelength  = 6.1  m 

(20  ft).  Amplitude  = 0.38  mm  (0.015  in.) 

4880 

150 

4.2 

0.0005 

1.07 

0.58 

2.6 

4880 

150 

4.2 

0.0010 

1.77 

0.58 

-5.8 

4880 

570 

4.2 

0.0005 

1.95 

1.50 

3.9 

4880 

570 

4.2 

0.0010 

2.39 

1.50 

-5.4 

Sinusoidal  Rails, 

Wavelength  = 12.2  m 

(40  ft)  , 

Amplitude  =0.38 

mm  (0.015  in. 

4880 

570 

4.2 

0.0005 

2.11 

1.50 

5.8 

4880 

570 

4.2 

0.00075 

4.80 

1.50 

-6.7 

4880 

570 

4.2 

0.0010 

6.34 

1.52 

-6.5 

Sinusoidal  Rails, 

Wavelength  = 24.4  iti 

1 (80  ft) , 

Amplitude  = 0.38 

mm  (0.015  in. 

4880 

570 

4.2 

0.0005 

4.44 

1.51 

-6.6 
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5.6  RESULTS  OF  AN  INITIAL  SCALING  STUDY 

Much  larger  guide-rail  systems  than  these  installed  in 
Tracks  K and  G are  currently  being  contemplated  at  AEDC.  A 
primary  concern  for  estimating  feasibility  of  such  track 
facilities  is  to  determine  how  projectile  wear  will  scale  with; 
projectile  diameter  and  velocity;  track  length;  rail  width;  and 
range  pressure. 

The  projectile  dynamics  computer  program  has  been  used 
to  conduct  an  initial  scaling  study  to  support  such  an  investi- 
gation. Some  results  obtained  during  the  course  of  this  study 
are  presented  in  this  section. 

The  values  of  u and  L used  for  the  scaling  study  were 
those  which  lead  to  good  agreement  between  computed  and  measured 
wear  for  Track  G projectiles,  namely  w = 0.00075  and  L = 4.2  x 
10®  J/kg.  The  value  of  k used  in  the  scaling  study  was  determined 
from  Equation  87.  The  basic  configuration  of  the  model  used  for 
the  scaling  study  was  provided  by  AEDC  staff  and  is  shown  in  Figure 
19.  Nominal  projectile  diameters  of  12.7  cm,  20.32  cm,  and  30.48  cm 
(5,  8,  and  12  inches  respectively)  were  considered  with  cone 
half-angles  of  8,  10, and  15  degrees.  Dimensional  and  other 
properties  of  these  projectiles  are  presented  in  Table  7.  The 
initial  gap  height  used  for  the  scaling  study  was  linearly  scaled 
with  projectile  diameter  from  the  Track  G value  of  0.5  mm.  The 
initial  lateral  velocity  of  the  projectile  was  fixed  at  0.6  m/s 
(2  ft/sec) . 

Some  typical  wear  histories  for  12.7  cm  diameter,  10° 
half-angle  cone  projectiles  are  shown  in  Figures  20,  21,  and  22. 

In  these  figures,  the  maximum  wear  on  each  of  the  wear  surfaces 
is  given  as  a function  of  distance  traveled  by  the  projectile. 

DRTF  represents  the  change  in  radius  (wear)  at  the  front  of  the 
top  wear  surface,  DRBF  represents  the  change  in  radius  at  the 
side  front  wear  surface.  No  collisions  occur  at  the  side  rails 
so  that  wear  at  the  side  surface  is  due  solely  to  aerodynamic 
heating  whereas  wear  at  the  top  and  bottom  surfaces  is  due  to 
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Figure  19.  Configuration  of  Composite  Projectile  Used  for  Scaling 
Study. 


DISTANCE  (feet) 


Figure  20.  Wear  History  of  a Nominal  13  cm  (5  inch)  Diameter 

Projectile.  Straight  Rails  with  Rail  Width  Linearly 
Scaled  from  Track  G. 
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TABLE  7.  PROPERTIES  OF  PROJECTILES  USED  IN  SCALING  STUDY 


Cone  Half 


Diameter 
(in)  (m) 

Angle 

(deg) 

Mass 

(kq) 

■"f 

(m) 

a 

(m) 

ro 

5.0 

0.127 

8 

4.41 

-0.0138 

0.1662 

0.0666 

5.0 

0.127 

10 

4.05 

+0.0054 

0.1470 

0.0446 

5.0 

0.127 

15 

3.55 

+0.0294 

0.1230 

0.0248 

8.0 

.203 

8 

18.06 

-0.0222 

0.2660 

0.6983 

8.0 

.203 

10 

16.60 

+0.0086 

0.2352 

0.4679 

8.0 

.203 

15 

14.56 

+0.0471 

0.1967 

0.2598 

12.0 

.305 

8 

60.95 

-0.0333 

0.3990 

5.301 

12.0 

.305 

10 

56.02 

+0.0130 

0.3528 

3.553 

12.0 

.305 

15 

49.13 

+0.0707 

0.2951 

1.973 

*Negative  means  center  of  gravity  is  located  on  forebody. 


both  sliding  friction  and  aerodynamic  heating.  Figure  20  shows 
that  the  projectile  has  nosed  down  and  is  sliding  on  the  front 
of  the  bottom  wear  surface  after  61  m (200  ft)  of  travel  and  re- 
mains in  that  position  during  the  remainder  of  its  travel.  This 
is  a phenomenon  which  was  observed  in  Track  G projectiles  at 
range  pressures  between  0.75  and  1.0  bars.  Figure  21  shows  that 
the  influence  of  periodic  rail  unevenness  is  negligible  when  the 
wavelength  is  30.5  m (100  ft).  This  is  consistent  with  the  re- 
sults given  in  Table  6 for  a Track  G projectile  where  it  was  ob- 
served that  the  wear  for  a wavelength  of  24.4  m (80  ft)  was  nearly 
identical  to  that  for  a corresponding  straight  rail  case.  However, 
smaller  wavelengths  can  cause  the  wear  to  be  more  evenly  distributed 
between  the  top  and  bottom  surfaces  while  at  the  same  time  re- 
ducing the  maximum  wear  experienced  by  the  projectile  as  is  evident 
from  the  results  given  in  Table  6.  Apparently  the  rail  unevenness 
can  prevent  the  projectile  from  sliding  on  just  one  rail.  This 
result  suggests  that  periodic  misalignment  should  purposely  be 
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designed  Into  long  tracks  to  keep  the  projectile  from  sliding 
along  a single  rail.  The  only  difference  introduced  between 
track-projectile  configurations  of  Figure  20  and  Figure  22  is 
that  the  rail  width  was  doubled  for  Figure  22.  The  rail  width 
used  in  Figure  20,  1.59  cm,  (0.625  inches)  represents  a linear 
scaling  from  the  current  Track  G track  configuration,  A com- 
parison of  Figures  20  and  22  shows  that  doubling  the  rail  width 
does  not  change  the  wear  due  to  aerodynamic  heating,  but  does 
reduce  the  wear  due  to  sliding  friction.  Of  course,  any  increases 
in  rail  width  will  increase  overall  projectile  drag  by  increasing 
shock  reflection  effects. 

The  results  of  the  initial  scaling  study  presented  here 
when  considered  together  with  the  comparison  of  the  predictions 
with  Tracks  K and  G measurements  indicate  that  the  Projectile 
Dynamics  Computer  Program  is  now  a useful  tool  for  establishing 
the  feasibilities  of  and  deisgn  for  very  large  track  systems. 
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SECTION  VI 

SUMMARY,  CONCLUSIONS,  AND  RECOMMENDATIONS 
6 . 1 SUMMARY 

The  program  undertaken  for  Arnold  Engineering  Development 
Center  and  described  in  this  report  involved  the  study  of  inter- 
faces between  the  projectile  and  rails  in  hypervelocity  guide-rail 
systems.  The  objective  of  the  study  was  to  provide  an  analysis  of 
the  interface  region  between  the  projectile  and  the  guide-rail  which 
would  be  capable  of  determining  the  friction  force  exerted  on  the  pro- 
jectile, the  wear  experienced  by  the  projectile  (and  possibly 
by  the  rail) , and  the  conditions  required  to  produce  gouging  of 
the  rails  during  projectile  passage.  The  analysis  must  be 
applicable  both  to  rail  systems  currently  in  operation  and  to 
future  systems  which  are  planned  to  be  far  larger  than  any 
currently  existing.  The  program  should  provide  capability 

for  treating  generalized  projectile  shapes,  although  special  em- 
phasis must  be  placed  on  the  treatment  of  slender  conical  fore- 
bodied  projectiles.  The  project  described  in  this  report  was 
successful  in  meeting  all  but  one  of  these  goals  .(the  prediction  of 
gouging) . In  addition,  the  project  has  produced  an  analysis  (and  an 
efficient  computer  program)  capable  of  describing  many  other  aspects 
of  hypervelocity  guidance  rail  operation. 

The  initial  analysis  of  the  projectile  guide  rail  interface 
was  based  on  two  assiimptions:  (1)  when  the  rail  system  is  op- 
erating properly,  the  interface  region  is  a self-pumped  gas 
bearing  which  may  or  may  not  be  affected  by  ablation  products 
from  the  projectile  and;  (2)  physical  contact  between  the  pro- 
jectile and  the  rails  is  a necessary  and  sufficient  condition  to 
produce  gouging  of  the  rails.  The  first  assumption  together  with 
fundamental  gasdynamic  analyses  have  been  used  to  identify  three 
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possible  types  of  gas  flow  within  the  interface  region.  These 
types  of  flow  were  identified  by  a Reynolds  niimber  calculated  on 
the  basis  of  the  average  separation  between  the  projectile  cuid 
the  rail.  Operation  in  two  of  these  regimes  can  be  expected 
for  practical  hypervelocity  guidance  rail  systems.  Each  of 
the  practically  reached  regimes  was  analyzed  thoroughly  and  com- 
puter subroutines  describing  them  were  generated. 

Early  in  the  project  it  was  realized  that  conditions  in 
the  interface  region  are  closely  coupled  to  the  dynamic  behavior 
of  the  projectile  along  the  guide-rails.  For  this  reason,  a gen- 
eralized analysis  of  hypervelocity  guide-rail  operation  was  de- 
veloped. This  analysis  was  formalized  in  a Projectile  Dynamics 
Computer  Program  which  combines  the  interface  subroutines  with 
other  subroutines  that  describe  the  rigid  body  motion  of  the  pro- 
jectile, the  aerodynamic  forces  experienced  by  the  projectile, 
and  the  detailed  geometry  of  the  hypervelocity  guide  rail. 

Later  in  the  program,  the  importance  of  elastic  deformation  of 
the  projectile  wear  surfaces  was  recognized  and  taken  into  account 
in  the  program.  The  program  can  be  used  to  compute  most  major 
parameters  required  to  describe  the  projectile  passing  along  a 
guide-rail  system.  Exercising  of  this  program  has  provided  con- 
siderable insight  into  the  overall  operation  of  hypervelocity 
guide-rail  systems.  The  relative  importances  of  simple  aerody- 
namic drag,  drag  induced  by  forces  from  the  projectile/guide-rail 
interface  region,  and  forces  produced  by  reflections  of  projectile 
bow  shock  waves  from  the  guide-rails  back  onto  the  projectile 
surfaces  were  evaluated.  The  importance  of  the  influence  of 
normal  forces  developed  by  elastic  deformation  of  the  projectile 
bearing  surface  on  wear  has  been  demonstrated  in  addition  to  the 
wear  produced  by  aerodynamic  heating.  Comparisons  of  some  im- 
portant aspects  of  the  Projectile  Dynamics  Computer  Program  pre- 
dictions were  made  with  experimental  results  obtained  from  hyper- 
velocity guide-rail  operations  at  Tracks  K and  G at  AEDC.  These 
comparisons  showed  good  agreement  between  computed  and  measured 
parameters,  including  total  drag,  wear,  and  rotation  of  projectiles 
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in  flight  along  the  rails.  Finally,  the  Projectile  Dynamics 
Computer  Program  was  used  to  conduct  an  initial  scaling  study  for 
investigating  the  feasibility  of  developing  very  large  hyper- 
velocity guide-rail  systems. 

No  specific  rail  gouging  criterion  was  developed  during 
the  course  of  this  project.  The  assumption  made  at  the  be- 
ginning of  this  study,  namely,  that  physical  contact  between  the 
projectile  and  the  rails  is  a necessary  and  sufficient  condition 
to  produce  gouging  proved  to  be  incorrect.  However,  with  the 
addition  of  elastic  deformation  of  the  projectile-bearing  surface, 
it  should  now  be  possible  to  develop  a normal  force  criterion  for 
gouging  to  occur  such  as  that  briefly  investigated  b^/  Graff  et  al. 

6.2  CONCLUSIONS 

The  project  described  in  this  report  produced  a number 
of  firm  conclusions  which  are  of  importance  to  the  operation  and 
development  of  hypervelocity  guide-rail  facilities.  These  con- 
clusions address  the  following  areas;  forces  and  pressures  gen- 
erated in  the  interface  region  between  a hypervelocity  projectile 
and  a guide-rail,  wearing  and  gouging  characteristics  of  inter- 
face areas;  and  hypervelocity  guide-rail  phenomena  not  directly 
connected  with  the  projectile  rail  interface  regions.  Conclusions 
in  this  last  area  are  important  to  the  overall  guide-rail  tech- 
nology although  they  lie  beyond  the  scope  of  this  project. 

6.2.1  Development  of  Projectile  Dynamics  Computer  Program 

A complete  computer  program  has  been  developed  with 
simulates  almost  all  known  aspects  of  elastically  deformable 
projectiles  traversing  a guide-rail  system.  Aspects  not  specifically 
covered  by  the  program  are  the  elastic  response  of  the  rails  and 
the  ablation  or  erosion  of  the  projectile  forebody  as  it  traverses 
the  rail  system. 

6.2.2  Interface  Flow  Regimes 

Three  separate  types  of  gas  flow  within  the  projectile/rail 
interface  region  have  been  identified.  The  type  of  flow  existing 
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in  any  particular  situation  is  determined  by  the  Reynolds  number 
Re^^  calculated  on  the  basis  of  the  gasdynamic  conditions  within 
the  gap  and  the  gap  separation.  The  three  types  of  flow  are: 

(1)  Reh<10^  leads  to  laminar  Couette  flow  in  the  inter- 
face region; 

3 5 

(2)  10  <Rej^<10  leads  to  turbulent  Couette  flow  in  the  in- 
terface region; 

5 

(3)  Rej^>10  leads  to  turbulent  boundary  layer  flow  in  the 
interface  region. 

Analyses  of  the  second  and  third  flow  regimes  have  been  emphasized 
since  these  are  the  only  ones  experienced  by  systems  of  practical 
importance.  The  first  regime  usually  only  occurs  with  such  small 
gap  separation  that  physical  contact  between  surface  asperities  al- 
most certainly  occurs. 

6.2.3  Fundamentally  Complex  Nature  of  Interface  Flow  Analysis 

Simple  analysis  can  not  be  developed  for  predicting 
projectile/rail  interface  parameters  because  the  frictional  shear 
forces  (and  hence  power  dissipation  within  the  gap)  are  strongly 
dependent  upon  gap  separation,  and  Re^^  which,  in  turn,  is  strongly 
dependent  upon  several  parameters  including  range  gas  density, 
projectile/rail  spacing,  flow  Reynolds  number,  and  the  geometry 
of  the  projectile  forebody. 

6.2.4  Relative  Magnitudes  of  the  Interface  Shear  Forces 

Projectile  drag  forces  arising  in  the  interface  regions 
are  negligible  compared  to  other  drag  forces  acting  on  the 
projectile  for  pro jectile/rail  separations  achievable  in  feasible 
systems . 

6.2.5  Effect  of  Interface  Gas  Pressures 

Internal  gas  pressure  within  the  projectile/rail  interface 
regions  increases  with  decreasing  gap  height,  but  in  general  does 
not  cause  sufficient  restoring  forces  to  prevent  the  projectile 
from  contacting  the  rails. 

6.2.6  Drag  Augmentation  By  Reflection  of  Shock  Waves 

The  predicted  total  drag  forces  for  projectiles  with 
both  blunt  and  slender  cone  forebodies  were  compared  with 
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experimental  results  from  the  Track  K at  AEDC.  Excellent  agreement 
was  found.  Detailed  analysis  of  the  slender  cone  drag  leads  to 
the  firm  conclusion  that  bow  shocks  reflected  from  the  rails  are 
weak  shocks,  and  the  less  firm  conclusion  that  the  reflected 
shock  from  the  rail  will  re-reflect  from  the  conical  surface  of 
the  projectile  and  contribute  significantly  to  the  total  drag  force. 

6.2.7  Importance  of  Rail  Imperfections 

In  minimizing  projectile  wear  one  of  the  more  surprising 
insights  obtained  from  the  initial  scaling  study  was  that  rail 
misalignment  could  lead  to  less  projectile  wear  than  the  wear 
obtained  in  tracks  with  perfectly  straight  rails.  Crooked  tracks 
can  insure  that  the  projectile  will  not  slide  for  too  long 
a time  on  a single  wear  surface.  The  influence  of  crooked  tracks 
is  to  more  evenly  distribute  the  wear  over  all  of  the  bearing 
surfaces.  More  total  mass  is  generally  removed  from  a projectile 
in  a crooked  track  than  in  one  with  straight  rails,  but  the 
maximum  depth  of  the  grooves  and  the  maximum  angular  rotation  of 
the  projectile  is  less. 

6.2.8  Versatility  of  the  Projectile  Dynamics  Computer  Program 

The  Projectile  Dynamics  Computer  Program  is  much  more 
versatile  than  its  application  to  date  would  indicate.  For 
example,  subroutines  treating  gasdynamics,  equilibrium  air  and 
transport  properties  can  be  used  together  as  a general  purpose 
program  for  studying  virtually  any  one  or  two-dimensional  gas- 
dynamic  problem  involving  air  flows  with  temperatures  in  the 
range  of  50“K  to  8,000®K. 

6 . 3 RECOMMENDATI ONS 

The  development  of  a powerful  computer  program  useful  for 
the  design  and  evaluation  of  current  and  future  hypervelocity 
rail  guidance  systems  provides  a variety  of  exciting  possibilities 
whose  exploitation  will  significantly  aid  the  overall  program 
for  developing  improved  ballistic  range  facilities  at  AEDC. 
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Additional  work  is  required  in  three  major  areas  in  order  to 
realize  these  potentials,  however.  First,  a gouging  criterion 
based  upon  determining  the  critical  normal  force  between  a pro- 
jectile and  a rail  which  will  lead  to  gouging  of  the  rail  should 
be  developed.  The  capability  of  determining  the  elastic  de- 
formation of  the  projectile  bearing  surfaces;  a capability  which 
was  added  to  the  program  near  the  end  of  this  project;  makes  it 
possible  to  now  use  this  capability  to  further  investigate  the 
gouging  mechanism.  Even  though  rail  gouging  has  proved  to  not  be 
a serious  problem  in  either  Tracks  K or  G,  this  damage  mechanism  can  be 
expected  to  be  a potentially  serious  problem  under  certain  circum- 
stances where  abnormally  large  normal  forces  are  generated. 

Second,  a study  should  be  performed  to  determine  optimum 
projectile  configurations  for  minimizing  wear.  The  location 
of  the  center-of-gravity,  wheelbase  (or  length  of  the  cylindrical 
portion  of  a projectile)  and  shock  interaction  length  are  some  of 
the  parameters  which  the  projectile  designer  can  control  and  these 
are  parameters  which  can  have  a significant  influence  on  pro- 
jectile wear.  The  Projectile  Dynamics  Computer  Program  can-  be 
used  to  assess  the  relative  wear  characteristics  of  rather 
general  axisymmetric  projectiles  and  with  minor  modifications 
could  also  be  used  to  assess  the  wear  characteristics  of  non- 
axisymmetric  projectiles.  There  exists  the  intriguing  possibility 
that,  with  a suitably  designed  projectile,  the  shock  interaction 
forces  can  be  used  to  stabilize  an  otherwise  unstable  projectile. 

Third,  further  effort  should  be  directed  toward  verifying 
the  true  nature  of  the  rail-shock  interaction  for  small  half-angle 
cones  (8®  to  15") . To  date  no  good  quality  schlieren  or  shadow- 
graph photographs  have  been  made  of  the  rail-shock  interaction 
region  for  cones  with  half-angles  less  than  30®.  To  our  knowledge, 
only  one  10®  half-angle  cone-cylinder  projectile  has  been  success- 
fully launched  in  Track  K and  none  have  been  launched 
in  Track  G.  Also,  apparently  none  of  the  existing  three-dimensional 
method-of-characteristics  computer  programs  have  the  capability 
of  analyzing  an  embedded  shock  system  as  complicated  as  the  rail- 
shock  system.  Good  quality  shadowgraphs  of  the  rail-shock 
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interaction  which  can  be  interpreted  unambiguously  cannot  be  made 
with  the  instrumentation  as  presently  installed  in  Tracks  K or 
G because  the  optical  path  through  the  flow  field  is  not  per- 
pendicular to  a plane  defined  by  two  opposing  rails.  It  is 
recommended  that  attempts  be  made  to  obtain  experimental  ver- 
ification of  the  nature  of  the  rail-shock  interaction  by  some 
means  so  that  the  validity  of  the  simplified  model  used  in  the 
projectile  dynamics  computer  program  can  be  determined  with 
greater  certainty  than  is  possible  from  drag  measurements. 
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APPENDIX  A 

SIMPLE  MODEL  FOR  PROJECTILE-RAIL  INTERFACE 
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SIMPLE  MODEL  FOR  PROJECTILE-RAIL  INTERFACE 


For  the  simple  interface  model  we  have  assumed  that  the  shear 
stress  ( T ) acting  on  the  projectile  is  proportional  to  the  cube  of  the 
velocity  (V)  and  inversely  proportional  to  the  rail -projectile  separation 
distance  (h).  We  have  assumed  the  pressure  acting  on  the  projectile 

is  inversely  proportional  to  the  square  of  the  separation  distance.  Thus 
we  have 


T = 


(A.  1) 


P 

T 


where  c and  e are  proportionality  constants. 


(A.  2) 


The  separation  distance  between  the  rail  and  projectile  at  any  point 
along  the  projectile  is  a function  of  the  transverse  position  of  the  projectile 
in  the  rail  system  and  the  orientation  of  the  projectile.  This  is  illustrated  in 
Figure  16.  The  letter  I denotes  some  position  along  the  projectile  center 
line  with  & = 0 corresponding  to  the  location  of  the  projectile  center -of -gravity. 
For  any  value  of  along  the  projectile  the  separation  distance  from  the  upper 
rail  (h^)  is  given  by 


h = z - z - £ sin  a - 1/2  d cos  0 (A,  3) 

11  P 
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Figure  A-1.  Separation  Distances  (At  Some  Position  £ Along  the  Pro- 
jectile) Between  the  Projectile  and  the  Upper  and  Lower 
Rails . 

where 

Zj  is  the  vertical  coordinate  of  the  upper  rail 

z is  the  vertical  coordinate  of  the  projectile  center  of 

gravity 

d^  is  the  diameter  of  the  projectile  body. 

The  separation  distance  from  the  lower  rail  (h  ) is 

h_  = z - z + X sin  9-  1 d cos  0,  (A.  4) 

2 2 2 P 

where  z^  is  the  vertical  coordinate  of  the  lower  rail. 
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In  order  to  determine  the  motion  of  the  projectile,  the  shear  stress 
and  pressure  must  be  converted  into  forces  and  torques  acting  on  the  projectile. 
The  shear  stress  will  produce  a force  (D^)  which  is  parallel  to  the  projectile 
body.  At  some  point  i along  the  projectile  body  the  shear  stress  acting  on  the 
projectile  is  t , + t _ « where  T is  shear  stress  from  the  upper  rail  and 
T Is  the  shear  stress  from  the  lower  rail.  The  incremental  force  AD 

2 I 

at  point  f is  given  by 


A°I  = 

(T^  + T^)  W^(Af) 

(A.  5) 

1 = - 

cU  ^ 

09 

(A.  6) 

1 

‘^l 

To  = 

c U ^ 

09 

{A.  7) 

2 

^2 

W is  the  width  of  the  rails,  and 
r 

A f is  an  Increment  along  i . 


The  value  of  is  found  by  integrating  the  incremental  forces  over  the 
body  of  the  projectile. 


D 


1 


+ 


W df 
r 


(A.  8) 


where 

i is  the  distance  from  the  center- of- gravity  to  the  forward  end 

^ of  the  projectile  body,  and 

I is  the  distauice  from  the  center- of- gravity  to  the  aft  end 

® of  the  projectile  body. 
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Substitutmg  Equations  (A.  3)  and  (A.  4)  into  Equation  (A.  8)  and  performing  the 
integration  yields: 


cvV 


D. 


sin  9 


/ h h 

I In'  la 


If  2a 


(A.  9) 


where 


If 

ha 

hf 

ha 


z,  - z-f,sine-T-d  cos  0 
1 t ^ P 


- *1 


+ 1 sin  9 - d cos  0 
a 2 p 


z - z_  +f  , sin  0 - -r  d cos  0 
£,  1 ^ P 

z - z_  - f sin  0 - 7 d cos  0 . 
^ 9l  ^ P 


(A.  10) 


The  factor  h ^ is  the  distance  between  the  upper  rail  and  the  forward  end  of  the 

projectile  body;  h is  the  distance  between  the  upper  rail  and  the  aft  end  of  the 
la 

projectile  body;  h_,  is  the  distance  between  the  lower  rail  and  the  forward  end  of  the 
projectile  body;  and  h_  is  the  distance  between  the  lower  rail  and  the  aft 
end  of  the  projectile  body.  It  might  be  noted  that  in  the  limit  as  0 approaches 
zero.  Equation  (A.  9)  becomes 


'0-0 


= cvV  (i  +<  ) (-^  ) IA.11) 

Ic  ^c 


where 

h is  the  distance  between  the  upper  rail  and  the  center  projectile  body, 
Ic 

h^  is  the  distance  between  the  lower  railand  the  projectile  body. 

The  pressure  in  the  interface  zone  exerts  a force  (L^)  perpendicular 
to  the  body  of  the  projectile.  The  pressure  from  the  lower  rail  (P2)  exerts 
a positive  force  on  the  projectile;  the  pressure  from  the  upper  rail  (P^)  exerts 
a negative  force  on  the  projectile.  The  incremental  force  (AL^)  at  some 
point  along  the  projectile  body  is 
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where 


The  totaL  force  Ls 


(P2  - Pj)  W^(Af) 


P = — 
1 . 2 


and 


P = — 
^2 


(A.  12) 
(A.  13) 

(A.  14) 


(A.  15) 


Substituting  Equations  (A.  13)  and  (A.  14)  into  Equation  (A.  15)  and  performing  the 
integration  yields 

I 1 


h = 


eW 
s 


in  9 h 


la  2a 

As  9 approaches  zero  Equation(A,  16)  becomes 


If 


2c 


Ic 


(A.  16) 


(A.  17) 


The  shear  stress  can  exert  a torque  (T  ) on  the  projectile  if  t and 

S X 

T2  are  not  equal  which  occurs  when  h^  is  not  equal  to  h^.  At  some  point 

along  the  projectile  body  the  incremental  torque  (^T  ) Ls 

s 


(A.  18) 


where  exerts  a positive  torque  and  exerts  a negative  torque.  The 
total  torque  due  to  the  shear  stress  is 

^£3  ^ 

= J-'  ' ^ ')  (A- 
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Substituting  Equations  (A.  3)  and  (A.  4)  into  Equation  (A.  19)  and  performing  the 
integration  yields  . 


cV  W d /K  K 

E__P  In  ‘ la-  2a 

2 sin  e 


/ la  2a  \ 


(A.  20) 


It  might  be  noted  that  the  limit  of  as  6 approaches  0 is 

= 4 W d ( f,  + f ) ( - 


®0  - 0 ^ 


r P 


Ic 


).  (A.  21) 


The  pressure  exerts  a torque  (T  ) on  the  projectile  if  P and  P_  are 

p 1 & 

not  equal.  The  incremental  torque  at  some  point  along  the  projectile  body  is 


AT  = (P_  - P.)  f W (Af) 
p 2 1 r 


(A.  22) 


where  P^  exerts  a positive  torque  when  i is  positive  and  P^  exerts  a negative 

torque  when  I is  positive.  The  total  torque  due  to  the  pressure  is 

f • 


= J 


) f W df 


(A.  23) 


- I 


Substituting  Equation  (A.  3)  and  (A.  4)  into  Equation  (A.  23)  and  performing  the 
integration  yields 


sW 


P ■ 2 - 
^ sm  0 


( 


h^.  h 


h,  h 


- — ^)+{z-z, -7d  cosexr^ 

2 2 p h^f 


) 


a If 


+ ( J.  -I  d^oo.  6)(^ 


2a 


la 


=17’’ 


(A.  24) 


It  might  be  noted  that  the  limit  of  T^  as  0 goes  to  zero  is 


P0-*  0 


eW  , , , 

-2-^ 

2 c 


1 c 


(A.  25) 
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Equations  (A.  10)  can  be  substituted  into  Equations  (A.  9).  (A.  16),  (A.  20), 

and  (A.  24)  to  define  the  forces  and  torques  acting  on  the  projectile  due  to  its 
interaction  with  the  rail  system.  These  equations  are  coupled  with  the 
equations  for  other  forces  and  torques  acting  on  the  p rojectile  to  describe 
the  motion  of  the  projectile. 
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APPENDIX  B 

A DESCRIPTION  OF  THE 
PROJECTILE  DYNAMICS  COMPUTER  PROGRAM 
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INPUT  VARIABLE  DEFINITIONS  FOR  THE  DYNAMIC 
RAIL/ PROJECTILE  SIMULATION  PROGRAM 


Variable 

Definition 

Units 

TEMPI 

Temperature  or  range  gas. 

°K 

PRESSI 

Pressure  of  range  gas. 

atmospheres 

ANGLEI 

Cone  half- angle. 

degrees 

FLAG  1 

Upper  limit  on  number  of  reflected 
waves  between  rail  and  projectile. 

Integer 

FLAG  2 

Controls  selection  of  weak  or  strong 
reflected  shocks.  If  set  equal  to 
zero,  weak  shocks  are  called  for 
unless  program  decides  that  this  is 
impossible.  If  set  equal  to  1,  then 
the  first  (and  only)  reflected  shock 
is  considered  to  be  a strong  shock. 

Integer 

DROC 

Diameter  of  projectile  body. 

m 

FLEN 

Distance  from  the  center- of- gravity 
of  the  projectile  to  the  front  sholder 
of  the  projectile  cylindrical  body. 

m 

ALEN 

Distance  from  the  center-of-gravity  of 
the  projectile  to  the  aft  end  of  the 
projectile  body. 

m 

MASS 

Mass  of  the  projectile. 

kg 

MOMIN 

Moment  of  Inertia  of  the  projectile 
about  an  axis  perpendicular  to  its 
longitudinal  axis  through  center-of- 
gravity. 

, 2 
kg*m 

WR 

Rail  width. 

m 

ZETA 

Nose  bluntness  ratio,  r /r.  . 

n b 

dimensionless 

ALNSI 

Distance  from  base  of  cone  to 
where  bow  wave  hits  rail. 

m 

CONELH 

Distance  from  base  of  cone  to 
nose  tip. 

m 

TETAS 

Shock  wave  angle  where  bow 
wave  hits  rail. 
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Figure  B-1.  Sketch  Showing  Definitions  of  Rail  and  Projectile  Input 
Dimensions . 


Figure  B-2,  Sketch  Showing  Definitions  of  ALNSl,  CONEDH,  TEXAS,  and 
DCYLCF. 
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Variable 

Definition 

Units 

DCYL.CF 

Diameter  of  cylinder  on  a hemisphere- 
cylinder- flare  projectile 

m 

IBLUNT 

Code  number  identifying  the  type 
of  projectile 

0 designates  a sharp  tip  cone 

1 designates  a blunted  cone 

2 designates  a hemisphere- 
cylinder- flare- -read  in 

CONELH  and  DCYLCF, 

ALNSI  and  TETAS  are 
calculated  in  subroutine  SOKSAP 

3 designates  a hemisphere- 
cylLnder-flare-- read  in  ALNSI, 
CONELH,  TETAS,  and  DCYLCF 

integer 

NT 

Number  of  segments  which  constitute 
the  top  rail  of  the  system. 

dimensionless 

NB 

Number  of  segments  which  constitute 
the  bottom  rail  of  the  system. 

dimensionless 

YYT 

The  y- coordinate  of  the  beginning 
point  of  the  top  rail.  The  x-axis 
runs  parallel  to  the  longitudinal 
direction  of  the  rail  system.  The 
y-axis  is  perpendicular  to  the  x-axis. 

m 

YYB 

The  y-coordinate  of  the  beginning 
point  of  the  bottom  rail. 

m 

IRULE  {1,1) 

Code  number  identifying  the 
geometric  shape  of  the 
top  rail  segment  as  follows; 

1 designates  a linear  segment 

2 designates  a circular  segment 

3 designates  a parabolic  segment 

4 designates  a sinusoidal  segment 

dimensionless 
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Variable 


VAR(1, 1) 


DIS(I,  1) 


Definition 


For  IRULE  (I,  1)  = 1,  the 
equation  for  the  rail  seg- 
ment is 

y = mx  + b. 

For  IRULE  (1, 1)  = 2,  the 
equation  for  the  rail  seg- 
ment is  A ■ ■ 

y = k +.\/r^-(x-h)2  . 

For  IRULE  (I,  1)  = 3,  the 
equation  for  the  rail  seg- 
ment is 


y 


k + 


(x-h)^ 

4f 


For  IRULE  (1, 1)  = 4,  the 
equation  for  the  rail  seg- 
ment is 


y = 


-A  + 
u y 


cos  ( 


2 TT  X ^ 


The  x-axls  is  parallel  to  the 
longitudinal  direction  of  the 
rail  system.  The  y-axis  is 
perpendicular  to  the  x-axis. 

When  IRULE  (1,1)  is  1,  VAR(I,  1) 
is  the  value  of  m in  the  equation 
defining  the  I^k  top  rail  segment. 

When  IRULE  (I,  1)  is  2,  VAR(I,  1) 
is  the  value  of  r in  the  equation 
defining  the  I^^  top  rail  segment. 

A positive  value  of  r indicates  the 
curvature  is  upward;  a negative 
value  of  r indicates  the  curvature 
is  downward. 


When  IRULEd,  1)  is  3,  VAR(1, 1) 
is  the  value  of  f in  the  equation 
defining  the  I^k  top  rail  segment. 

When  IRULE  (1, 1)  is  4,  VAR(1, 1) 
is  the  value  of  y^  in  the  equation 
defining  the  I^h  top  rail  segment. 

The  X position  at  which  the  I^k 
top  rail  ends. 


Units 


dimensionless 


m 


m 


m 


m 
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Variable  Definition  Units 


CX(I,1) 

When  IRULE  (1, 1)  is  1, 

CX(1, 1)  is  the  value  of  b 
in  the  equation  defining  the 
top  rail  segment. 

m 

When  IRULE  (1,1)  is  2, 

CX(1, 1)  is  the  value  of  h 
in  the  equation  defining  the 

Ith  top  rail  segment. 

m 

When  IRULE  (1, 1)  is  3, 

CX(1, 1)  is  the  value  of  h 
in  the  equation  defining  the 
ith  top  rail  segment. 

When  IRULE  (1,1)  is  4, 

CX(1, 1)  is  the  value  of  A 
in  the  equation  defining  tKe 
ith  top  rail  segment. 

CY(I,1) 

When  IRULE  (1, 1)  is  1, 

CY(1, 1)  is  not  used; 
any  value  can  be  read-  in. 

When  IRULE  (1, 1)  is  2, 

CY(1, 1)  is  the  value  of  k 
in  the  equation  defining  the 
ith  top  rail  segment. 

m 

When  IRULE  (1, 1)  is  3, 

CY(I,  1)  is  the  value  of  k 
in  the  equation  defining  the 
ith  top  rail  segment. 

m 

When  IRULE  (1, 1)  is  4, 

CY(1, 1)  is  the  value  of  L 
in  the  equation  defining  the 
ith  J.QP  segment. 

m 

IRULE(I,2) 

VAR(I,2) 

DIS(I,2) 

CX(I,2) 

CY(I,2) 

These  variables  are  defined 
the  same  as  IRULE  (I,  1), 

VAR(I,1),  DIS(I,  1),  CX(I,  1), 
andCY(1, 1),  respectively, 
except  they  apply  to  the  ith 
bottom  rail  segment  rather 
than  the  Ith  top  rail  segment. 

X(l) 

The  Initial  x- coordinate  of  the 
projectile  center- of- gravity 
(usually  0. 0). 

m 
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Variable 

X(2) 


X(3) 

X(4) 

X(5) 

X(6) 

IPCHK 


IV  A R 


TO 

TF 


Definition 

The  Initial  /-coordinate  of  the 
projectile  center- of- gravity, 

(Must  be  between  the  rails  In 
such  a manner  that  the  projectile 
body  does  not  extend  outside  the 
rails). 

The  Initial  value  of  the  angle  ( 0 ) 
between  the  projectile  longitudinal 
axis  and  the  x-axls. 

The  Initial  value  of  x. 

The  Initial  value  of  y. 

(Usually  0,0), 

The  Initial  value  of  6. 

(Usually  0.  0) , 

A value  of  1 Indicates  that 
the  output  should  be  printed 
at  specified  time  Intervals. 

A value  of  2 Indicates  that 
the  output  should  be  printed 
at  specified  Intervals  of  x. 

Indicates  which  variable  should 
be  checked  to  determine  If  the 
numerical  Integration  step  should 
be  adjusted.  The  choices  are 

IVAR  = 1,  check  on  the  rate  of 
change  of  x. 

IVAR  = 2,  check  on  the  rate  of 
change  of  y, 

IVAR  = 3,  check  on  the  rate  of 
change  of  0. 

IVAR  = 4,  check  on  the  rate  of 
change  of  x. 

IVAR  = 5,  check  on  the  rate  of 
change  of  y. 

IVAR  = 6,  check  on  the  rate  of 
change  of  0. 

The  beginning  time  of  the  simulation. 
(Usually  0,  0), 

The  final  time  of  the  simulation. 


Units 

m 


radians 

m/sec 

m/se,c 

rad/ sec 

dimensionless 


dimensionless 


sec 


sec 
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Variable 

Definition 

XF 

The  final  X- distance  of  the 
simulation. 

TPRNT 

The  time  interval  at  which 
output  will  be  printed  if 

IPCHK  = 1. 

XPRNT 

The  distance  interval  at  which 
the  output  will  be  printed  if 

IPCHK  = 2. 

CP 

The  upper  limit  of  central  processor 
time  to  be  used. 

GEAR 

The  mean  aerodynamic  chord 
associated  with  the  projectile. 

CDO 

The  zero  angle  of  attack  drag 
coefficient  for  the  projectile. 

CDA 

The  change  of  projectile  drag 
coefficient  with  change  in  angle 
of  attack. 

CMO 

The  zero  angle  of  attack  pitching 
moment  coefficient  for  the  projectile. 

CMA 

The  change  in  projectile  pitching 
moment  coefficient  with  change  in 
angle  of  attack. 

CLO 

The  zero  angle  of  attack  lift 
coefficient  of  the  projectile. 

CLA 

The  change  In  projectile  lift 
coefficient  with  change  in  angle 
of  attack. 

CLAD 

The  change  in  projectile  lift 
coefficient  with  the  time  rate  of 
change  of  the  angle  of  attack. 

CMAD 

The  change  in  projectile  pitching 
moment  coefficient  with  the  time 
rate  of  change  of  the  angle  of  attack. 

CLQ 

The  change  in  projectile  lift 
coefficient  with  the  change  in 
pitching  rate  (6). 

CMQ 

The  change  in  the  projectile  pitching 
moment  coefficient  with  the  change 
in  pitching  rate  (6). 

Units 

m 

sec 

m 

sec 

m 

dimensionless 
per  radian 

dimensionless 
per  radian 

dimensionless 
per  radian 

per  radian 

per  radian 

per  radian 

per  radian 
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Variable 

G 


DELTA T 

HMIN 

HFAC 


VALl 


VA  L2 


Definition 

The  acceleration  due  to  gravity. 

If  the  two  rails  being  simulated 
lie  in  a horizontal  plane  or  if 
gravity  is  to  be  ignored,  G should 
be  read-in  as  0.0. 

The  initial  time  step  of  the 
numerical  integration. 

The  minimum  time  step  to  be 
used  in  the  numerical  integration. 

The  factor  by  which  the  time  step 
is  adjusted.  If  an  increase  is 
needed,  the  present  time  step  will 
be  multiplied  by  HFAC.  If  a decrease 
is  needed,  the  present  time  step  will 
be  divided  by  HFAC. 

The  value  used  to  determine 

if  the  time  step  should  be  increased. 

The  test  is  made  on  VAL  where 


VAL 


RRKOl  - RRK(2) 

RRK(2)  - RRK(l) 

RRK(  i)  is  the  estimate  of  the 
rate  of  change  in  the  variable  designed 
by  IVAR,  Four  estimates  of  the  rate 
of  change  of  each  variable  are  made 
at  each  time  step  by  the  fourth  order 
Runge-Kutta  numerical  integration 
scheme.  If  VAL  is  less  than  VALl, 
the  time  step  is  increased  by  multiplying 
the  present  time  step  by  HFAC.  A nominal 
value  of  VALl  might  be  0.  008,  but  this  may 
depend  upon  the  value  of  IVAR  and  the 
particular  system  being  simulated. 


The  value  used  to  determine  if  the  time 
step  should  be  decreased.  If  VAL  is 
greater  than  VAL2,  the  time  step  is 
decreased  by  dividing  the  present  time 
step  by  HFAC.  A nominal  value  of 
VAL2  might  be  0.08,  but  this  may 
depend  upon  the  value  of  IVAR  and  on 
the  particular  system  being  simulated. 


Units 

m/sec^ 

sec 

sec 

dimensionless 

dimensionless 


dimen  s ionle  s s 
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Variable 

HABLAT 

HABL 

TWALL 

RHOLEX 

AMUSL 

STIFFC 


Definition 

Total  heat  of  ablation,  L. 

Heat  of  ablation  due  to  heat 
of  vaporization  and  heat  of 
polymerization . 

Ablation  temperature. 

Density  of  projectile  material. 

Coefficient  of  sliding  friction, 

ij. 

Elastic  deformation  coefficient, 

k . 

e 


Units 
Joules /kg 
Joules /kg 

°K 

kg/m^ 

N/m^ 
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A BRIEF  DESCRIPTION  OF  THE  FUNCTION  OF  THE  PROGRAM 
SUBROUTINES 

Section  V provides  a brief  description  of  the  Projectile  Dynamics 
Computer  Program  (PDCP).  Figure  15  shows  the  basic  design  of  the  pro- 
gram. The  user  should  note  that  the  basic  system  of  units  used  in  the  PDCP 
is  the  MKS  system.  All  of  the  gasdynamic,  thermodynamic,  and  transport 
property  subroutines;  however,  use  the  CGS  system  of  units.  The  user 
should  also  note  that  specific  values  of  extensive  thermodynamic  quantities 
are  used  in  two  different  ways  in  the  program.  When  quantities  such  as 
enthalpy,  entropy  and  specific  heat  at  constant  pressure  appear  in  gasdynamic 
subroutines,  they  are  expressed  on  a per  gram  basis, whereas  when  they  appear 
in  strictly  thermodynamic  subroutines  these  quantities  appear  on  a per  gram- 
mole  basis.  A brief  description  of  each  of  the  37  individual  subroutines  follows. 
Where  applicable,  references  to  sources  of  information  used  to  develop  a 
particular  subprogram  are  given.  The  listing  of  the  program  is  replete 
with  comment  cards  which  attempt  to  describe  many  of  the  details  of  the 
methods  and  procedures  used  in  developing  the  various  subroutines. 

RAIL 


This  is  the  main  program.  The  first  function  performed  is  to  call 
INPUT  which  reads  in  and  initializes  quantities.  With  all  quantities  initialized, 
RAIL  then  calls  RK  which  carries  out  the  numerical  integration  of  the 
equations  of  motion.  After  RK  has  completed  each  integration  step,  all 
information  on  the  time  elapsed,  the  position  of  the  center-of-mass  of  the 
projectile,  the  first  and  second  derivatives  of  x,  z,  and  fl  and  the  wear  ex- 
perienced by  the  top,  bottom  and  side  rails  during  the  integration  time  step 
are  returned  to  RAIL.  RAIL  then  stores  this  information  along  with  the 
results  of  the  previous  two  integration  steps.  RAIL  then  checks  to  see  if 
the  output  should  be  printed.  If  necessary  the  Lagrangian  interpolation 
subroutine  is  called  if  the  output  is  desired  at  a time  or  distance  inter- 
mediate to  the  present  value  of  time  or  distance  travelled.  The  reason  for 
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always  having  three  sets  of  output  stored  is  so  that  the  LagrangLan  inter- 
polation can  be  carried  out,  RAIL  terminates  the  run  when  time  TF  is 
reached  or  when  the  distance  XF  is  reached.  If  the  velocity  of  the  pro- 
jectile decreases  by  more  than  Z0%  from  its  origiani  value,  RAIL  reinitializes 
the  gasdynamic  conditions  and  and  to  the  current  conditions. 

It  is,  of  course,  possible  to  update  these  quantities  at  the  end  of  each  in- 
tegration step,  but  since  these  quantities  do  not  change  drastically  or 
strongly  influence  sensitive  quantities  such  as  interface  normal  forces  and 
moments,  it  is  not  necessary  to  perform  these  computations  at  the  end  of 
each  integration  step. 

INPUT 

This  subroutine  reads  in  all  input  data.  It  computes  the  Mach 
number  of  the  projectile  from  the  input  data  and  then  computes  the  shock 
angle  of  an  attached  conical  shock.  This  computation  of  the  shock  angle  is 
performed  by  means  of  a function  which  is  an  empirical  formula  developed 
by  the  authors  based  on  a curve  fit  to  real  gas  solutions  to  the  conical 
shock  equations.  SOKSAP  is  then  called  to  compute  the  thermodynamic 
and  gasdynamic  conditions  at  key  locations  on  the  projectile.  The  initial 
angle-of-attack  is  computed  and  then  DRAGCO  is  called  to  determine  initial 
values  of  Cn  , and  Ct  , and  Cm  • HEATRS  is  called  to  determine  an 
average  aerodynamic  heating  heat  transfer  coefficient  for  the  interface. 

RK 

This  subroutine  is  in  essence  a standard  Runge-Kutta  fourth-order 
numerical  integration  subroutine  having  variable  time  step  capability.  The 
Runge-Kutta  method  basically  computes  four  separate  predictions  of  the 
unknowns  to  arrive  at  new  values  of  the  unknowns.  During  each  prediction, 

RK  first  computes  the  angle-of-attack.  Then  it  calls  ORIENT  to  determine 
the  orientation  of  the  projectile  relative  to  the  rails.  If  any  of  the  separation 
distances  are  negative,  ELASTC  is  called  to  determine  the  elastic  deformation 
forces  and  moments.  With  the  gap  height  established  for  the  interfaces,  INTER 
is  called  to  determine  the  net  interface  force  components  and  the  net  moment 
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due  to  the  Interface  forces  and  rail  shock  Interaction  forces.  Subroutine 
DRAG  is  then  called  to  determine  the  total  aerodynamic  lift,  drag,  and  moment 
coefficients  at  the  current  angle-of-attack.  Subroutine  F is  then  called  which 
determines  the  force  components  in  the  x and  z directions.  After  the  four 
predictions  on  the  unknowns  have  been  computed  , RK  checks  to  see  if  the  time 
step  used  was  too  large.  If  it  was,  the  time  step  is  reduced  and  the  entire 
procedure  indicated  above  is  repeated.  If  the  time  step  becomes  less  than  a 
minimum  time  step  read  in  as  input,  the  program  is  terminated  in  RK.  If 
the  time  step  is  adequate,  a weighted  average  of  the  four  predictions  for  each 
unknown  is  computed  which  defines  the  new  values  of  the  unknowns.  The  new 
values  of  the  unknowns  are  then  returned  to  RAIL. 

F 

This  subroutine  simply  determines  the  total  net  force  on  the  pro- 
jectile in  the  x and  z directions  and  the  net  torque. 

DRAG CO 

The  basic  function  of  this  subroutine  is  to  determine 
and  determined  for  both  the  conical  and  cylindrical 

portions  of  the  projectile  using  formulas  taken  from  reference  28.  The  method 
by  which  CDQ'*iB  determined  is  fully  described  in  Section  4.1. 

DRAG 

Given  the  angle-of-attack  and  pitching  rate,  this  subroutine  de- 
termines the  total  drag  coefficient,  C , the  total  lift  coefficient  C , and 

the  total  moment  coefficient  about  the  center  of  gravity,  C . 

m 

INTER 

The  complete  version  of  this  subroutine  is  given  in  Section  4.2  and 
the  simplified  version  is  described  in  Appendix  A. 

TURBFP 

In  this  subroutine,  the  Karman-Schoenherr  turbulent  skin  friction 
coefficient  equation  is  solved  for  a given  value  of  the  Reynolds  number.  The 
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equation  being  solved  Is  Equation  53  which  is  an  Implicit  equation.  Newton's 
iterative  method  for  solving  for  the  roots  of  a transcendental  equation  is  used 
to  determine  C^. 

TURBCE 

This  subroutine  provides  the  skin  friction  coefficient  for  turbulent 

Couette  flow  given  the  Mach  number  and  Reynolds  number.  Re,  . The  Kor- 

h 

kegi- Briggs  implicit  skin  friction  equation  for  equal  heat  transfer  rates 
(Equation  51)  is  solved  by  use  of  Newton's  iterative  procedure. 

ORIENT 

This  subroutine  determines  the  interface  gap  heights  at  four  different 
locations  on  both  the  top  rail  interface  and  the  bottom  rail  interface.  The  four 
locations  are  the  front  of  the  interface,  and  the  three  positions  on  the 
cylindrLcal  portion  of  the  projectile  (see  Figure  17).  Given  the  current 
position  of  the  projectile  (x,z,  and  0 of  the  center-of-gravity)  and  geometric 
quantities  related  to  the  dimensions  of  the  projectile  this  subroutine  determines 
the  space  coordinates  of  the  four  positions  indicated  above  for  both  the  top  and 
bottom  of  the  projectile.  Subroutine  CONFIG  is  then  called  to  determine  the 
space  coordinates  of  the  top  and  bottom  rails  at  these  four  locations.  From 
all  of  this  information,  ORIENT  then  computes  the  gap  heights.  Note  that 
gap  heights  are  not  computed  for  the  side  rail  interfaces  since  the  program 
does  not  account  for  lateral  motion  of  the  projectile. 

If  a negative  gap  height  is  obtained,  ORIENT  signals  RK  through  the 
value  of  ICOLIN  that  this  has  occurred.  RK  then  calls  ELASTC  . ORIENT 
prints  out  a message  that  an  apparent  collision  with  the  wall  has  occurred 
and  tells  where  on  the  projectile  the  apparent  collision  has  occurred. 

ELASTC 

This  subroutine  is  called  from  subroutine  RK  whenever  ORIENT 
determines  that  the  distance  between  a bearing  surface  and  the  rail  has  be- 
come negative,  based  upon  rigid  body  behavior,  ELASTC  then  computes  by 
means  of  the  equations  given  in  Appendix  C,  the  normal  forces  required  to  de- 
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form  the  projectile  bearing  surfaces  to  zero  separation  distance.  The 
moments  associated  with  these  normal  forces  are  also  computed  in 
ELASTC  . 

LGRNG 

This  is  a standard  3 point  Lagrangian  interpolation  subroutine. 

CONFIG 

This  subroutine  gives  the  space  coordinates  of  either  the  top  or 
bottom  rail  at  any  position  x along  the  rail  system.  The  shape  of  the  rails 
is  defined  by  input  information  which  is  fully  described  at  the  beginning  of 
this  Appendix. 

GASDYN 

This  subroutine  solves  normal  and  oblique  shock  problems  for  air 

and  serves  as  an  executive  subroutine  for  carrying  out  other  gasdynamic 

computations.  This  subroutine  along  with  the  equilibrium  air  subroutines 

and  the  undissociated  air  subroutines  can  be  used  over  a temperature  range 
o o 

of  50  K to  about  8500  K.  At  the  low  end  of  the  temperature  range  the 
user  has  to  use  some  judgment  since  the  program  cannot  tell  if  or 

will  liquify.  At  the  upper  end  of  the  temperature  range  some  judgment 
must  also  be  used  since  the  ionization  of  atomic  species  is  not  included  in 
the  equilibrium  air  subroutines. 

Some  of  the  details  of  this  subroutine  have  been  presented  in  Section 
3.4.1  and  will  not  be  repeated  here.  The  first  computation  performed  in 
the  subroutine  is  to  determine  the  conditions  immediately  downstream  of  the 
bow  shock  where  it  hits  the  rail.  Given  the  shock  angle  at  the  rail  and  initial 
conditions  GASDYN  determines  the  turning  angle  of  the  flow  and  the  downstream  flow 
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conditions.  The  second  computation  is  to  determine  the  conditions  which 
exist  on  the  cylindrical  portion  of  the  projectile.  This  is  accomplished 
by  a call  to  subroutine  PMEXP  which  performs  a Prandtl -Meyer 
expansion  from  conditions  existing  behind  the  conical  shock  to  conditions 
which  exist  behind  a flow  which  has  turned  through  an  angle  equal  to  the  cone 
half -angle.  The  third  computation  is  to  determine  the  conditions  down- 
stream of  the  reflection  of  the  bow  shock  from  the  rail.  The  conditions  up- 
stream of  this  first  reflected  shock  are  those  just  previously  computed 
and  the  turning  angle  is  the  same  as  the  turning  angle  of  the  flow  as  it 
passed  through  the  bow  shock.  The  shock  angle  for  this  first  reflected 
shock  is  unknown.  An  estimate  of  the  shock  angle  based  on  a solution  to 
the  perfect  gas  with  constant  specific  heat  oblique  shock  equations  is  ob- 
tained from  subroutine  CUBIC.  The  oblique  shock  equations  are  then  solved 
(three  equations  in  three  unknowns)  as  described  in  Section  3.4.  1.  The 
program  will  ordinarily  assume  that  weak  shock  solutions  are  desired  unless 
a strong  shock  solution  is  specifically  requested  by  setting  FLAG  2 equal 
1 (see  the  first  part  of  this  Appendix  for  a full  description  of  the  input 
quantities).  If  the  turning  angle  is  too  large  for  a given  Mach  number,  a 
weak  shock  solution  is  not  physically  possible  in  which  case  the  program 
will  do  a strong  shock  solution.  The  user  specifies  the  upper  limit  on  the 
number  of  weak  reflected  shocks  that  he  wants  to  consider  by  FLAG  1. 
GASDYN  will  determine  the  conditions  behind  each  of  these  reflected  waves; 
or  if  regular  reflection  is  not  possible,  for  a given  shock,  it  will  do  a strong 
shock  solution  and  tell  the  user  the  number  of  the  shock  which  has  been  con- 
sidered to  be  a strong  shock  which,of  course,  must  be  the  last  shock  then 
regardless  of  what  the  value  of  FLAG  1 is.  The  initial  conditions  plus  con- 
ditions behind  each  subsequent  shock  are  stored  in  COMMON.  After  con- 
ditions behind  all  shocks  are  determined,  conditions  at  other  points  of  in- 
terest are  determined.  First,  sonic  conditions  for  the  flow  behind  the  last 
shock  are  established  by  a call  to  CONSTS.  Next,  conditions  at  state  4 
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(conditions  after  expansion  from  conditions  behind  the  last  shock  in  the 
system  to  the  pressure  existing  behind  the  attached  conical  shock,  P^) 
are  determined.  If  the  flow  is  still  supersonic  after  the  last  shock,  then 
PMEXP  is  called  to  do  this  expansion.  If  the  flow  behind  the  last  shock  is 
subsonic,  then  PMEXP  is  called  to  expand  the  flow  from  sonic  conditions 
to  P^  if  the  pressure  at  sonic  conditions  is  greater  than  P^  • 
pressure  at  sonic  conditions  is  less  than  P^,  CON  STS  is  called  to  isentro- 
pically  compress  the  air  from  sonic  conditions  to  P^.  Finally,  stagnation 
conditions  corresponding  to  the  flow  behind  the  last  shock  are  determined 
by  another  call  to  CONSTS, 

The  thermodynamic  and  flow  properties  at  each  of  the  conditions 
described  above  are  stored  in  COMMON,  Table  5 shows  how  the  various 
points  are  labelled, 

PMEXP 

This  subroutine  performs  Prandtl-Meyer  expansions  and  com- 
pressions, Three  options  are  provided  for  handling  three  different  kinds 
of  processes.  The  calling  agrument  of  this  subroutine  includes  a variable 
called  IPM  and  another  variable  called  ISTATE,  If  IPM=1,  the  turning  angle 
of  the  flow  is  specified  (through  DELNU),  If  IPM=2,  the  final  pressure  is 
specified  (through  PF),  If  IPM=3,  the  final  Mach  number  is  specified  (through 
AMF),  The  initial  conditions  are  specified  through  ISTATE  which  designates 
the  number  of  the  initial  state.  Initial  state  conditions  are  presumed  to  be 
available  in  COMMON  labelled  STATE.  Conditions  at  the  end  of  the  Prandtl- 
Meyer  process  are  also  stored  in  STATE  at  the  location  specified  by  FSTATE 
which  is  also  a part  of  the  calling  argument  for  PMEXP.  Values  of  the 
Prandtl-Meyer  function,  the  thermodynamic  properties,  and  the  flow  velocity 
and  Mach  number  are  stored  at  each  integration  step  so  they  can  be  used 
in  other  subroutines  to  help  define  flow  conditions  at  the  entrance  to  the 
interface  region. 

■ The  first  part  of  the  subroutine  uses  ideal  gas  relations  to  solve  Prandtl- 
Meyer  flows.  The  second  part  of  the  subroutine,  which  is  used  for  high 
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TABLE  1.  ISTATE  DESIGNATIONS 


The  gasdynamlc  conditions  at  each  state  in  the  system  are  stored  in 
labeled  COMMON  designated  by  STATE.  The  following  numbering  system 
for  labeling  the  states  is  used  in  the  program. 


State  No. 
1 
2 

3 

4 

5 

6 

7 

8 


9 


10 


Gasdynamic  Description  of  State 

Free- stream  conditions. 

Behind  conical  (or  first)  shock. 

Behind  first  reflected  shock  (first  reflection 
off  rail). 

Behind  second  reflected  shock  (first  reflection 
off  cone). 

Behind  third  reflected  shock  (second  reflection 
off  rail). 

Behind  fourth  reflected  shock  (second  reflection 
off  cone). 

Conditions  after  Prandtl-Meyer  expansion  around 
cone -cylinder  shoulder  (conditions  along  side  of 
rail-projectile  interface). 

Dummy  - stores  intermediate  state  results.  If 
conditions  behind  last  shock  are  subsonic,  then 
conditions  after  expansion  to  M=1  are  stored 
here.  Also  used  to  store  intermediate  results 
when  conditions  behind  the  last  shock  are  highly 
supersonic.  In  this  case,  if  CONSTS  subroutine 
does  not  converge  when  attempting  to  locate  sonic 
conditions  the  results  after  Prandtl-Meyer 
expansion  to  M=4  are  stored  here. 

Conditions  after  expansion  from 
pressure  existing  behind  last  shock 
to  pressure  existing  behind  the 
first  shock. 

Stagnation  conditions  behind  the  last 
shock  in  the  system. 
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temperature  conditions,  uses  numerical  integration  to  solve  the  Prandtl- 
Meyer  flow  pressure  as  the  integration  variable.  CONSTS  is  used  in  the 
numerical  integration  to  determine  the  change  in  thermodynamic  and  flow 
properties  with  each  incremental  change  in  pressure.  A Modified  Simpson's 
Rule  is  used  for  the  numerical  integration.  The  equation  which  is  integrated 
is  given  by 


d ( Jtn  P) 
oV^/P 


where  dw  is  the  angle  through  which  the  flow  turns  as  the  result  of  an  in- 
cremental change  in  pressure  d(  AnP). 

This  subroutine  actually  possesses  much  more  capability  than 
is  used  in  the  present  version  of  the  Projectile  Computer  Program.  For 
example,  this  subroutine  Together  with  expressions  for  determining  the 
displacement  thickness  of  the  three-dimensional  boundary  layers  on  the 
rail  and  the  projectile  could  be  used  to  define  much  more  precisely  the 
pressure  distribution  in  the  interface  for  cases  where  a significant  inviscid 
core  region  exists  between  the  two  boundary  layers.  This  subroutine  could 
be  used  to  perform  detailed  computations  of  the  circumferential  expansion 
of  the  inviscid  core  flow  in  the  interface. 


CUBIC 

This  subroutine  is  used  to  solve  the  cubic  equation  which  relates  the 
shock  wave  angle  to  Mach  number  and  specific  heat  ratio  for  two-dimensional 
flow  over  a wedge  of  a perfect  gas  with  constant  specific  heat.  This  sub- 
routine is  called  to  provide  an  estimate  of  the  shock  angle  when  the  Mach 
number  and  turning  angle  are  known.  The  equation  solved  in  this  sub- 
routine is  Equation  150a  in  Reference  28.  This  subroutine  also  includes 
a hypersonic  similarity  law  empirical  formula  that  gives  the  shock  angle 
in  terms  of  the  Mach  number  and  wedge  half-angle.  This  formula  is  used 
when  the  wedge  angle  is  large  and  the  Mach  number  is  also  large. 
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EQAIR 

This  subroutine  solves  for  the  equilibrium  composition  of  dissociated 
(but  not  ionized)  air  given  the  temperature  and  pressure  of  the  air.  The  com- 
position of  the  air  in  a non -dissociated  state  is  taken  to  be  dry-  air  having 
a molar  composition  of  20.  979%  O^i  78.  087%  N^,  and  0.  934%  Argon. 

The  formula  weight  of  one  mole  of  this  air  is  28,  963  grams/gram-mole. 

Some  of  the  details  of  this  subroutine  are  discussed  in  Section  3.4.  1,  The 
basic  formulation  of  the  coupled  set  of  equilibrium  thermod-ynamic  equations 
and  the  method  of  solution  is  too  lengthy  to  present  here.  The  details  are 
given  in  Reference  30,  Suffice  it  to  point  out  that  the  general  development 
follows  that  originally  presented  in  Reference  31.  A great  amount  of 
thought  and  effort  has  gone  into  EQAIR  in  order  to  produce  a computationally 
efficient  program.  The  procedures  used  in  this  subroutine  are  not  straight- 
forward or  easy  to  understand.  The  Newton-Raphson  procedure  is  used 
to  solve  a system  of  six  coupled  non-linear  algebraic  equations.  The 
equations  consist  of  three  chemical  equilibrium  equations  for  the  following 
reactions  '• 

20 

« 2N 

1/2  0^+  1/2  :3=S=  NO, 

two  mass  balance  equations:  one  for  conservation  of  oxygen  atoms  and  one 
for  conservation  of  nitrogen  atoms,  and  Dalton's  law  of  partial  pressure. 

IGUESS 

This  subroutine  determines  initial  estimates  of  the  composition  of 
the  equilibrium  air  given  the  temperature  and  pressure  of  the  air.  These 
initial  estimates  are  used  in  EQAIR  to  start  the  Newton-Raphson  iterative 
procedure.  The  equations  used  to  determine  the  initial  estimates  of  the 
unknowns  were  developed  along  lines  similar  to  those  discussed  in  Reference 
30. 
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SONIC 

This  subroutine  determines  the  speed  of  sound  of  equilibrium 
air  given  the  temperature,  density,  and  composition  of  the  air  as  deter- 
mined in  EQAIR.  The  method  used  is  based  on  that  presented  in  Reference 


31. 

DERIV 


This  subroutine  determines  the  thermodynamic  derivatives 

which  are  used  in  GASDYN  to 


carry  out  the  Newton -Raphson  procedure  applied  to  solution  of  the  normal 
shock  and  oblique  shock  equations.  This  subroutine  calls  CSOLVR  to  solve 
a system  of  linear  algebraic  equations.  The  unknowns  in  one  set  of  linear 


equations  are  the  thermodynamic  derivatives  a ( 

known  in  another  set  of  linear  algebraic  equations 
derivative 


ni)~i 

""  J 

are  th 


and  the  un- 


e thermodynamic 


ivativesT  3(  Ai  njM  . These  thermodynamic  derivatives  together  with 

L d An  T)  Jp 

temperature,  pressure,  density,  composition,  and  individual  species 
thermodynamic  properties  allow  one  to  determine  the  four  thermodynamic 
derivatives  defined  at  the  beginning  of  this  paragraph.  The  general  method 
used  here  is  described  in  Reference  31. 


CONSTS 

The  basic  function  of  this  subroutine  is  to  perform  thermodynamic 
process  computations  for  both  dissociated  air  and  non-dis sociated  air. 

This  subroutine  can  also  be  used  to  determine  thermodynamic  properties 
when  conditions  other  than  temperature  and  pressure  are  known.  For 
instance,  this  subroutine  is  used  in  INTER  to  determine  the  temperature 
corresponding  to  a state  where  the  pressure  and  enthalpy  are  known  (that 
is  to  determine  the  temperature  corresponding  to  the  Sommer -Short  ref- 
erence enthalpy  at  a pressure  equal  to  the  pressure  at  the  edge  of  a boundary 
layer).  The  primary  function  of  this  subroutine;  however,  is  to  determine 
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conditions  at  the  end  state  of  an  isentropic  expansion  or  connpression.  The 
subroutine  contains  5 basic  options  defined  byr  IFLG.  If  IFLG=1,  the  final 
entropy  and  pressure  are  specified,  if  IFLG  =2,  the  final  entropy  and  tem- 
perature are  specified,  if  IFLG=3,  the  final  entropy  and  enthalpy  are  specified, 
if  IFLG=4,  the  final  entropy  and  Mach  number  are  specified,  and  if  1FLG=5, 
the  final  pressure  and  enthalpy  are  specified.  In  each  option,  given  the 
specified  quantities  and  initial  velocity  or  Mach  number,  the  other  end  point 
thermodynamic  properties  and  velocity  and/or  Mach  number  will  be  deter- 
mined. The  development  of  this  subroutine  is  based  on  the  method  outlined 
in  Reference  3 1.  One  non-obvious  procedure  used  to  implement  option  five 
needs  to  be  explained.  The  initial  enthalpy  must  be  read  in  as  HT  and 
the  specified  final  enthalpy  must  be  read  in  as  MWT. 

CSOLVR 

This  is  a linear  equation  solver  subroutine.  The  method  used  is 
Gaussian  elimination  with  pivotal  condensation. 

TRANSP 

This  subroutine  determines  viscosity  and  thermal  conductivity 
of  high  temperature  air  given  the  temperature  and  mole  fractions  of  the 
constituents  of  the  air.  The  methods  used  in  this  subroutine  are  taken 
directly  from  Reference  22. 

TCHEM 

This  subroutine  determines  the  individual  specie  ideal  gas  thermo- 
dynamic properties  at  the  input  temperature  and  a pressure  of  one  atmos- 
phere. The  change  in  the  Gibbs  function;  change  in  entropy,  and  change 
in  enthalpy  for  each  of  the  three  reactions  discussed  in  the  description  of 
EQAIR,  evaluated  at  one  atmosphere,  are  also  computed  in  this  subroutine 
as  well  as  the  equilibrium  constants  for  the  three  reactions.  This  sub- 
routine serves  as  an  executive  subroutine  with  most  of  the  computations 
actually  performed  in  the  following  four  subroutines  which  represent  a 

statistical  thermodynamic  package  for  the  determination  of  thermodynamic 

*•  32 

properties 
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SPECF 

This  subroutine  calculates  thermodynamic  properties  of  N^.  O^, 

O,  N,  NO  and  Ar.  Specific  heats,  enthalpy,  and  entropy  at  one  atmosphere  are 
computed  from  partition  functions  for  translation,  rotation,  vibration,  and 
electronic  excitation  contributions.  At  high  temperatures,  anharmonicity 
corrections  are  applied. 

TRANS 

Evaluates  the  translational  partition  function  for  each  species. 

VIBR 

Evaluates  the  vibrational  partition  function  for  each  species. 

CORREC 

Evaluates  anharmonicity  corrections  for  each  species. 

MU STAR 

A function  subprogram  which  returns  the  viscosity  of  undissociated 
air  given  the  temperature  of  the  air, 

KSTAR 

A function  subprogram  which  gives  the  thermal  conductivity  of 
undissociated  air. 

SIG 

A function  subprogram  which  gives  the  Prandtl  number  of  un- 
dissociated air. 

CPSTAR 

A function  subprogram  which  determines  the  specific  heat  of  un- 
dissociated air. 

GAMMA 

A function  subprogram  which  gives  the  specific  heat  ratio  of  un- 
dissociated  air. 
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ENT 

This  function  determines  by  means  of  table  lookup  the  enthalpy 
of  undissociated  air  considered  to  be  a perfect  gas.  This  function  returns 
an  enthalpy  divided  by  a constant  given  by  . 0686132  x 273.  16.  The  data 
in  the  table  were  calculated  with  SPECF  of  this  program  and  are  nearly 
identical  to  those  given  in  Reference  33.  The  temperature  range  is  50  to  3000®K. 
The  table  is  in  increments  of  10"K  up  to  800®K  and  incrementsof  SO^K  from 
800 ’K  to  3000'’K. 

KELVIN 

This  is  the  inverse  function  of  ENT.  It  returns  the  temperature  of 
air  given  the  enthalpy  divided  by  0,  0686132  x 273.  16. 

SOKSAP 

The  primary  function  of  this  subroutine  is  to  determine  the  shape  of 

the  bow  shock  and  compute  both  the  axial  location  where  the  bow  shock  would 

intersect  the  rail  if  the  gap  height  were  zero  and  the  angle  of  the  bow  shock  at 

this  location.  The  method  used  to  estimate  the  shape  of  the  bow  shock  is  that 

given  by  Inouye.^"*  The  shock  standoff  distance  is  computed  by  Inouye's 

method  after  GASDYN  has  been  called  to  determine  the  density  ratio  across 

the  normal  shock.  After  the  position  and  angle  of  the  bow  shock  at  the 

radial  distance  r have  been  determined,  the  quantities  ALNSI  and  TETAS 
b 

defined  in  Figure  B-2,  GASDYN  is  called  to  determine  the  gasdynamic  con- 
ditions at  the  states  defined  in  Table  1 of  this  Appendix. 

HEATHS 

This  subroutine  is  used  to  determine  boundary  layer  heat  transfer 
coefficients.  The  methods  used  are  taken  directly  from  Appendix  C of 
Reference  27. 
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APPENDIX  C 

LINEAR  ELASTIC  DEFORMATION  ANALYSIS 
OF  PROJECTILE  BEARING  SURFACES* 


*This  analysis  was  developed  by  Dr.  J.B.  Hicks  of  the  Aeroballistic 
Branch  of  VKF  at  AEDC. 
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LINEAR  ELASTIC  DEFORMATION  ANALYSIS 
OF  PROJECTILE  BEARING  SURFACES"!' 


The  basic  concept  used  to  determine  the  elastic  deformation  of  the 
projectile  bearing  surface  is  to  allow  the  surface  to  deform  whatever  amount 
is  necessary  to  prevent  a negative  separation  distance  from  occurring 
between  the  projectile  and  the  rail.  The  elastic  deformation  analysis  of 
the  projectile  is  not  coupled  to  the  gas  bearing  analysis,  that  is,  no  elastic 
deformation  is  considered  when  separation  distances  are  positive. 


Subroutine  ORIENT  determines  the  separation  distances  at  the  front 

and  rear  of  the  top  and  bottom  bearing  surfaces.  These  distances  are 

defined  in  Appendix  A as  h , h , h , and  h where  the  subscripts  1 and  2 

XI  XcL  ul 

refer  to  the  top  and  bottom  bearing  surfaces,  respectively  while  the  sub- 
scripts f and  a refer  to  the  front  and  rear  (aft)  ends  of  the  bearing  surfaces, 
respectively.  Whenever  any  of  these  separation  distances  are  negative 
subroutine  ELASTC  is  called  to  determine  the  normal  force  required  to 
deform  the  bearing  surface  to  zero  separation  distance. 


Three  simplifying  assumptions  are  made  in  the  development  of  the 

elastic  deformation  analysis.  First,  the  force  per  unit  area  P required 

to  deform  the  surface  by  an  amount  6 is  assumed  to  follow  a linear  re- 
ef 

lationship  of  the  form  given  by 


P 

ef 


= k 6 . 

e e£ 


(C-1) 


where  k is  an  elastic  coefficient.  Second,  the  deformation  over  the  width 
e 

of  bearing  surface  is  assumed  to  be  constant.  Third,  the  deformation  over  the 
length  of  the  bearing  surface  is  assumed  to  be  linear.  The  total  normal  force 
is  then  obtained  by  integrating  Equation  (C-1)  over  the  length  of  the  bearing  sur- 
face where  negative  separation  distances  occur. 


'i'This  analysis  was  developed  by  Dr.  J.B.  Hicks  of  the  Aeroballistic  Branch 
of  VKF  at  AEDC. 
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The  methods  used  to  arrive  at  the  working  equations  used  in 
subroutine  ELASTC  based  upon  integrating  Equation  C- 1 are  next  described. 
Three  separate  cases  are  defined  for  determining  the  working  equations 
depending  on  the  signs  of  the  separation  distances.  These  three  cases  are 
defined  for  the  top  surface  as  follows: 


Case  1 : 

^a 

£ 0 and  h^^ 

£0 

Case  2 : 

h. 

& 0 and  h. , 

< 0 

la 

If 

Case  3 : 

h. 

< 0 and  h , 

aO 

la 

If 

Similarly,  three  cases  are  defined  for  the  bottom  surfaces. 

Case  1 for  Top  Bearin0  Surface  (h^^^  and  h^^  both  negative) 

The  following  loading  diagram  sketch  for  Case  1 (top  surface)  is 
used  to  define  certain  quantities  and  illustrate  the  method  used  to  determine 
the  total  normal  force  and  moment. 
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For  Case  1 this  sketch  Illustrates  that  the  total  loading  can  be  considered 

as  the  sum  of  a constant  loading  and  a variable  (linear)  loading.  The 

normal  force  due  to  the  uniform  loading,  F is  simply 

c 


F = - W Bk  h, 
c r e la 


while  the  normal  force  due  to  the  linear  loading,  F^  is 


F 

r 


-- W Bk  (h,,-  h,  ) 
2 r e If  la 


(C-2) 


(C-3) 


The  total  normal  load,  F is  then 

n 


F = - F - F 
n c r 


The  total  normal  load  on  the  top  bearing  surface  F is  then  (with  the  same  sign 

nl 

convention  as  used  in  Section  V and  Appendix  A). 


F 


nl 


■^W  Bk 
2 r e 


(3  h 


la 


- h^^) 


(C-4) 


The  moment  about  the  center- of- gravity  of  F is  the  sum  of  the  moment  due 

n 

to  F^  and  F^  and  can  be  written  by  inspection  of  the  loading  diagram  as 

T.  = -F  (f-«a)-F  f|B-jia)  (C-5) 

i C V 

where  again  the  sign  convention  is  that  used  in  Section  V and  Appendix  A. 

Case  2 for  Top  Bearing  Surface  (h^^.  negative , h^^  positive) 

The  loading  diagram  for  Case  2 is  as  follows: 
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The  normal  force  and  moment  can  very  simply  be  expressed  in  terms  of 

the  distance  from  the  rear  of  the  projectile  to  the  point  where  the  loading 

is  zero,  x , From  the  law  of  similar  triangles,  x can  be  written  as 
o 


X 

o 


h 


la 


B 


{C-6) 

or 


X 

o 


B 


and  from  inspection  of  the  loading  diagram 

W k h 


_ r e If  I 

F = - F = (B  - X ) 

nl  1 2 o 


(C-7) 


and  the  moment  about  the  center-of-gravity  is 


T = - F I X + - (L  - X ) - A 
Ij  1 I o 3 o a 


] 
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which  can  be  simplified  to 


T = - F (—  B 
13^ 


) 


(C-8) 


Case  3 for  Top  Bearing  Surface  negative,  h^^  positive 
The  loading  diagram  for  Case  3 is  as  follows 


The  law  of  similar  triangles  applied  to  this  loading  diagr2un  leads  to 

O o 


•h 

la  If  la 


or 


(C-9) 


X = 

o 


(1 


) 


and  from  inspection  of  the  loading  diagram 


F = W k h X 
1 2 r e la  o 


(C-10) 
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and  the  moment  is 


\ « a- 


(C-11) 


The  normal  loads  and  moments  due  to  elastic  deformation  of  the 
bottom  bearing  surface  for  each  of  the  three  cases  can  be  written  by  in- 
spection of  the  loading  diagrams  for  each  corresponding  top  surface  case. 
The  same  sign  convention  is  used  for  the  bottom  surface  as  for  the  top, 
that  is,  forces  are  positive  in  the  upward  direction  and  counterclockwise 
moments  are  positive. 


Case  1 for  Bottom  Surface  (h^^  and  h^^  both  negative) 


F = - W Bk  h, 
c r e 2a 

(C-12) 

F = - W Bk  (h, , - h,  ) 

V 2 r e 2f  2a. 

(C-13) 

' = - 4 W Bk  (3h-  - h,,) 

nZ  Z r e 2a  Zf 

(C-14) 

B ZB 

= F^  i a)  + Fv  (Y 

(C-15) 

Case  2 for  Bottom  Surface  (h^^  negative  and  h_^  positive) 


X 

o 


B 


(1 


2f 

'Za 


(C-16) 


F , = W k (B  - x^)  h 
nZ  2 r e o 2f 


(C-17) 


T = F , + Y - -I  ) 

n2  3 3 a 


(C-18) 
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Case  3 for  Bottom  Surface  negative  and  positive) 


B 

~^Zf 


(C.19) 


F , = --  Wxkh, 
nZ  Z r o e Za 


(C-ZO) 


T = - F U ) 

nZ  a 3 


(C-Zl) 


The  elastic  deformation  analysis  developed  here  is  incorporated 
in  the  computer  program  in  subroutine  ELASTC.  Note  that  this  analysis 
is  based  upon  assuming  that  the  interface  region  can  be  represented  by  one 
straight  line  segment.  This  procedure  is  thus  inconsistent  with  the  method 
used  to  track  the  wear  of  the  top  and  bottom  interface  regions  as  described 
in  Section  5.4.  Two  straight  line  segments  are  used  to  describe  the 
geometry  of  the  top  and  bottom  wear  surfaces.  Therefore,  an  obvious 
improvement  in  the  accuracy  of  the  program  could  be  made  by  using  two 
straight  line  segments  to  define  the  shape  of  the  elastically  deformed  interface. 
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Area  of  interface 

Ablating  contact  area 

Interface  area  affected  by  lubrication 

Sound  speed 

Length  of  projectile  bearing  surface 

Increase  in  B caused  by  wear  into  the  conical  surface 

Correction  coefficient  to  account  for  mass  added  to 
boundary  layer 

Proportionality  constant  for  shear  stresses 
Distance  from  center-of-gravity  to  center-of-pressure 
Aerodynamic  drag  coefficient 
Aerodynamic  lift  coefficient 

Skin  friction  coefficient  for  laminar  Couette  flow 

Average  skin  friction  coefficient 

Aerodynamic  pitching  moment  coefficient 

Numerical  factors  involving  pressure  drag  defined 
after  Equation  35 

Coefficient  defining  the  change  in  the  drag  coefficient 
with  a change  in  angle  of  attack 

Base  drag  coefficient 

Drag  coefficient  from  pressure  on  the  forebody 

Drag  coefficient  arising  from  the  interference  of 
the  shock  wave  by  the  rails 
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iYi/i 


C and 

”‘q 


m 


SI 


p max 


Total  projectile  drag  coefficient 

Viscous  drag  coefficient  on  entire  projectile 

Drag  coefficient  at  zero  angle  of  attack 
Lift  coefficient 

Coefficient  defining  the  change  in  the  lift  coefficient 
with  a change  in  angle  of  attack 

C_  Damping  coefficients  associated  with  the  lift 

Lift  coefficient  arising  from  the  interference  of 
the  shock  wave  by  the  rails 

Lift  coefficient  at  zero  angle  of  attack 

Moment  coefficient 

Lift  coefficient  at  zero  angle  of  attack 

Coefficient  defining  the  change  in  the  pitching 
moment  coefficient  with  a change  in  angle  of  attack 

C Damping  coefficients  associated  with  the  pitching 
"*01  moment 

Pitching  moment  coefficient  arising  from  the 
interference  of  the  shock  wave  by  the  rails 

Term  defined  especially  for  Equation  36 

Aerodynamic  drag  force 

Frictional  force  between  projectile  and  rail 
Projectile  diameter 

Proportionality  constant  involving  transverse 
pressure 

Viscous  force 

Average  viscous  force 

Normal  force  due  to  elastic  deformation 
Sliding  friction  force 
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LIST  OF  SYMBOLS  (Continued) 

Lubrication  force 

Net  tangential  force  due  to  interface  phenomena 

Proportionality  constant 

Acceleration  of  gravity 

Gas  enthalpy 

Reference  enthalpy 

Gap  heighth 

Minimum  gap  heights 

Gap  heighth  to  top  rail 

Gap  heighth  to  bottom  rail 

Distance  between  upper  rail  and  rearward  end  of 
projectile 

Distance  between  upper  rail  and  center  of  projectile 

Distance  between  upper  rail  and  forward  end  of 
projectile 

Distance  between  lower  rail  and  rearward  end  of 
projectile 

Distance  between  lower  rail  and  center  of  pro'jectile 

Distance  between  lower  rail  and  forward  end  of 
projectile 

Minimum  gap  heighth 

Pitching  moment  of  inertia  of  projectile 
Boltzmann  constant 

Thermal  conductivity  of  interface  gas 
Elastic  coefficient 
Heat  of  ablation 
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LIST  OF  SYMBOLS  (Continued) 

Aerodynamic  lift  force 

Lift  force  from  projectile/rail  interface 
Distance  along  projectile  body  from  center  of  gravity 
Distance  from  rear  of  projectile  to  center  of  gravity 
Distance  from  front  of  projectile  to  center  of  gravity 
Increment  along  projectile  position  parameter,  1 
Mach  Number 
Projectile  mass 
Molecular  weight 

Individual  molecular  mass,  or  projectile  mass 
Ablative  mass  loss 

Ablation  mass  loss  rate  from  sliding  friction 
Number  density  of  atoms 
Static  pressure 

■•ig. 

Average  static  pressure 

Transverse  gas  pressure  in  interface 

Average  pressure  on  cone  surface 

Interface  pressure  caused  by  ablation  products 

Local  elastic  deformation  force 

Pressure  on  projectile  from  upper  rail  interface 
Gas  pressure  in  region  2 defined  in  Figure  6 
Energy  dissipation  from  sliding  friction 
Pitch  rate  (0) 

Nondimensional  heat  flux  at  the  rail 
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LIST  OF  SYMBOLS  (Continued) 

Heat  flux  to  ablative  wall 

Gas  constant  for  ablation  products 

Reynolds  Number 

Lubrication  Reynolds  Number 

Reynolds  Number  based  on  length  B 

Reynolds  Number  based  on  gap  heighth,  h 

Reynolds  Number  based  on  free  stream  conditions  and 
cone  slant  length 

Reynolds  Number  based  on  Ax 

Reynolds  Number  based  on  x 

Cone  base  radius 

Nose  radius 

Cross  sectional  area  of  projectile 
Entropy 

Gas  temperature 
Aerodynamic  torque 

Torque  from  projectile/rail  interface 

Torque  on  projectile  from  pressure  couple  in  interface 

Torque  on  projectile  from  shear  couple  in  interface 

Ablation  temperature 

Rail  temperature 

Range  gas  temperature 

Sliding  velocity 

Local  gas  velocity  component  along  rail  axis 
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LIST  OF  SYMBOLS  (Continued) 

V Projectile  velocity 

V Local  gas  velocity  component  along  circumference 
of  projectile 

W Projectile  velocity  component  normal  to  rail 

W Width  of  rails 

r 

w Local  gas  velocity  component  normal  to  rail 

W'  Average  bearing  pressure  above  inlet  pressure 

Ax  Distance  along  rail  from  strong  -shock  to  gap  entrance 

Axc-  Distance  from  cone  base  to  reflected  shock  im- 

pingement point 

X Denotes  direction  of  projectile  travel 

y^  Laminar  boundary  layer  thickness 

z Vertical  coordinate  of  projectile  center  of  gravity 

z^  Vertical  coordinate  of  upper  rail 

Zj  Vertical  coordinate  of  lower  rail 

GREEK  SYMBOLS 

A Bearing  number 

a Angle  of  attack 

6 Define  parameter  for  Equation  51 

V Angle  between  velocity  vector  and  rail  axis 

Local  elastic  deformation 

a Cone  half-angle 

6 Angle  between  longitudinal'  axis  of  projectile 

and  rail  axis 

Shock  angle  from  rail  axis 
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ii 
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n 


C 

P 

T 

T. 


W 


xz 


yz 


Shock  angle 

Shock  reflection  angle 
Gas  viscosity 

Coefficient  of  sliding  friction 

Kinematic  viscosity  using  pressure  from  region  n 
of  Figure  6 and  temperature  of  rail 

Kinematic  viscosity  (n  refers  to  gas  region  from 
Figure  6 ) 

Bluntness  ratio 

Gas  density 

Equivalent  shear  stress  of  interface  region 

Stress  defined  in  Equation  45 

Viscous  shear  stress  in  x direction 

Viscous  shear  stress  in  y direction 

T evaluated  on  upper  rail 

T evaluated  on  lower  rail 

Non-dimensional  ablative  mass  loss 

Prime  specifies  turbulent  fluctuations  of  quantity 


SUBSCRIPTS 


aw 

b 

e 

number 


w 


adiabatic  wall  conditions 

Projectile  base  conditions 

Conditions  at  edge  of  boundary  layer 

Gas  dynamic  conditions  for  regions  defined  in 
Figure  6 and  Section  3.3.2 

Free  stream  conditions 

Conditions  evaluated  at  a surface 
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